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ERRATUM 


The Scattering of Deuterons and Protons by Alpha-Particles, by R. G. 
FREEMANTLE, T. GRoTDAL, W. M. Grsson, R. McKnacue, D. J. PROWSE, 
and J. Rorsnat, 1954, Phil. Mag., 45, 1093. 


An incorrect Q-value for the reaction *He(d, p)>He was included 
inadvertently ; the value actually obtained was 3-3-+ 0-1 Mey. 
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CXXII. A Search for Natural Radioactivity in Neodymium, 
Rhenium and Osmium 


By D. Dixon and A. McNair 
Department of Natural Philosophy, University of Glasgow* 


[Received May 6, 1954] 


ABSTRACT 


The element neodymium and the naturally occurring isobaric pair 
Rhenium 187 and Osmium 187 are investigated, as internally mounted 
sources in a well-shielded proportional counter, for evidence of natural 
radioactivity. The B-decay probability of Neodymium 150 corresponds 
to a half-life of not less than 101° yr. No beta instability of rhenium is 
observed and 1015 yr. is suggested as the minimum possible half-life for 
the isotope '8*Re. The emission of L x-rays from osmium is observed, 
but it is found that dense materials show similar emission, and it is 
attributed to the photoelectric absorption of y-rays (residual background) 
in the materials. It is concluded that there is no evidence of decay 
in osmium and the estimated minimum possible half-life by electron 
capture exceeds 1015 yr. 


§1. InTRODUCTORY 


THERE have been various conflicting reports concerning the radioactivity 
of naturally occurring neodymium. In 1934 Libby attributed a 
B--activity of half-life 1-510! yr. and maximum energy 11 kev to 
neodymium and in later years evidence, reviewed by Mulholland and 
Kohman (1952), began to accumulate in favour of 6--activity. More 
recently Curran et al. (1952), in a preliminary survey of neodymium 
failed to detect any instability within the limits of the sensitivity of 
their apparatus, and a similar investigation by Mulholland and Kohman 
(1952), who laid stress on the need for chemical purification in the 
preparation of sources, increased the minimum possible half-life of °°Nd 
to 2x10!5 yr. It was suggested that the B--decay of 1°°Nd, which 
might be expected from a study of nuclear f-stability limits (Kohman 
1948), was highly forbidden, due to a large spin change involved in the 
transition. Interest in the possible instability of naturally occurring 
1445Nd was heightened on account of its position in the nuclear table 
with regard to the possible ‘missing’ stable isotopes of promethium. 
Considerations of the systematics of stable nuclei of odd Z (>7) restrict 
the occurrence of f-stable isotopes of promethium (Z=61) to odd A values 
eee ese ek te ee a ee er ee ee 
* Communicated by S. C. Curran, F.R.S. 


4C2 


1100 D. Dixon and A. McNair on a Search for 


and further limitations are expressed by the position of known artificially 
produced radioactive isotopes of Pm. The occurrence of a f--stable 
isotope, 14°Pm, was considered, and in such circumstances the B--decay 
of its neighbouring naturally occurring isobar 145Nd should have been 
energetically possible (Feather 1952). The discovery, however, by 
Butement in 1951 of the captive-activity of 14°Pm, established indirectly 
the B--stability of “°Nd. The isotopes “°Nd and 1478m thus remain as 
one of three instances of pairs of B-stable nuclei having the same odd 
number of neutrons (Z>7), stability being achieved seemingly by a 
single even number of protons in all other nuclei of odd N values. 
Evidence of alpha-instability in isotopes in the region N=84, is provided 
in several instances among artificially produced radio-isotopes (‘Thompson 
et al. 1949, Rasmussen et al. 1950) and the possibility of «-emission from 
naturally occurring neodymium isotopes (as suggested by Kohman 1949), 
was investigated by Mulholland and Kohman (1952) and an activity of 
less than 0-02 «’s per sec per g was obtained, compared with 0-3 «’s 
per sec per g obtained by Bestenreiner and Broda in 1949. Recently 
Waldron, Schultz and Kohman (1953) using impregnated nuclear plates, 
have found a natural «-activity in neodymium and have ascribed it to 
the isotope “4Nd (N=84). This report on the investigation of the 
possible radioactivity of neodymium was completed in March 1953, 
with the exception of a renewed search for «-activity, which was under- 
taken following a private communication from Kohman reporting an 
a-activity of “4Nd of 1-9+0-1 Mev energy and of half-life 1-5 1015 yr. 
In our earlier experiments a very weak «-activity of 1-6 e.p.m. from 
1-18 g of Nd was observed and was thought to arise mainly from 
slight radioactive contamination. 


§2. EXPERIMENTAL ARRANGEMENT 


Previous investigations of the radioactivity of some natural radio- 
elements (Curran e¢ al. 1952) showed that the maximum possible efficiency 
of detection was required, and our experimental equipment has been 
improved towards this end. The system now in operation (and to be 
described in detail in a subsequent paper) consists essentially of a large 
specially designed proportional tube spectrometer completely surrounded 
by two rings of G.M. counters (72 cm in length) in anticoincidence, and 
shielded by 8 in. of lead and 1 in. of iron on all sides. This system is 
situated in the most shielded part of the building, some 10 ft. below 
ground level. The proportional counter is of the end-corrected type 
(Cockroft and Curran 1951) with 10in. of its cathode fully operative 
and with a diameter of 5-65 in. giving a maximum internal source area 
of 1100cm?. The inner walls of the counter are lined with stainless 
steel and, together with other refinements, this helped to reduce the total 
background count rate to less than 16 c.p.m., a value of the same order 
(for comparative counting volumes) as that observed by Libby e¢ al. 
(1951) and Kulp and Tryon (1952) in their experiments on !@C. This 
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has made possible the more successful investigation of certain radio- 
elements (Dixon et al. 1954) and represents a reduction in the background 
by a factor of 20 relative to that observed in most of our earlier experi- 
ments. The pulses from any possible radioactive source are displayed 
oscillographically and a photographic record is analysed in amplitude 
to yield the pulse energy spectrum. The proportional counter, used with 
the anticoincidence system described above, is insensitive to the detection 
of radiations for only 1% of the time operated. The K x-rays of silver 
were used for energy calibrating purposes and the resolution, at high gas 
pressures, corresponded to a full width at half-height equal to 12% of 
the pulse amplitude. This resolution was maintained throughout the 
period of operation of the counter and energy calibrations made at 
intervals during the course of experimental runs up to 12 h duration, 
invariably remained constant. 


Fig. 1 


ENERGY INTERVAL 


COUNTS PER 


ENERGY INTERVAL 


PER 


so 


COUNTS 


° 


Spectra of background (indicated by broken histogram) and background 
+1-18 g of neodymium (full histogram). 


§ 3. NEODYMIUM 


The source of neodymium was mounted on the steel inner lining of 
the counter in the form of NdCl,, of thickness 2-1 mg/cm?, obtained from 
the pure metal. The results of the spectrum analysis for the energy 
ranges 6-60 kev and 20-300 kev are shown in fig. 1. The counting 
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rates of the total background and background +-source at lower energies 
suggest that if the material is radioactive the half-life Is not less than 
1015 yr. In view of the complete absence of any appreciable difference 
in the pulse spectra of background and background source it is likely 
to be much longer. We conclude that 15°Nd is effectively B-stable, a 
conclusion which is not in contradiction with any known nuclear stability 
rules. Following the communication from Kohman concerning the 
suggested «-activity of “4Nd, an attempt was made to measure the 
energies of the a-rays previously observed in the experiments described 
above. Due to the unfavourable nature of the source and its thickness 
(2:1 mg/em?) our results, which did however indicate the presence 
of a-rays of energy <2-0Mev (~1-5 c.p.m.) must, in the main, be 
regarded as inconclusive. Earlier investigations of low energy «-rays 
following the natural decay of !47Sm (Dixon 1953) illustrated the success 
that can be achieved in energy measurements of «-activities of long 
life-time, using thin uniform sources of finely divided powder internally 
mounted in proportional spectrometers. Energy resolution of 5% 
(width at half maximum height) was obtained for the 2-22 Mev «-rays 
from 17Sm of half-life 1011 yr. It is hoped to re-examine Nd for 
z-activity in more favourable conditions. 


§4. RHENIUM AND OsMIUM 
Introductory 

The isotopes 18’Re and 18’Os form one of the few instances of naturally 
occurring neighbouring isobars, so that one should be radioactive. Some 
of the earlier studies by Scherrer and Zingg (1939), Lougher and Rowlands 
(1944), suggested that 18’Os was radioactive by K- and L-capture, but 
later results did not confirm this reported activity. The examination 
of both rhenium and osmium by Naldrett and Libby (1948) indicated 
that 18’Re was B--active with a half-life, 7, 4-1 x 1012 yr., whilst Sugerman 
estimated the half-life to be almost twice as great. In their examination 
of osmium, Seliger and Bradt (1948) failed to find any significant 
activity. More recently Gauthé and Blum (1953) used nuclear plates to 
confirm the activity of 187Re and work done in Glasgow seemed to do 
likewise, but in our case there was a notable difference in the value 
obtained for the maximum f-ray energy. Naldrett and Libby proposed 
43 kev on the basis of absorption measurements, while a preliminary 
report by Curran (1952) on our own investigation of the form of the 
spectrum suggested that the f-rays had a limiting energy of about 
400 kev. Since that time we have considerably reduced our ‘ background ’ 
as indicated above, and in the process found that most of the f-particles 
previously thought by us to be emitted from 187Re were due entirely 
to the differing radioactive contamination of the aluminium sheets on 
which the source had been mounted, and that which was used as the inner 
lining of the proportional counter. This result indicates the need for 


very careful and constant checking of all materials used in the mounting 
of sources. 
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Hxperimental Arrangement 


The same apparatus was used as in the earlier neodymium experiments. 
Two sources of thickness 0-686 mg/cm? and 0-423 mg/cm? were prepared 
using spectrographically pure rhenium metal and both were mounted in 
turn on the steel lining of the counter in the form of finely divided metallic 
powder held in suspension in a solution of amyl acetate, in which was 
dissolved a trace of nylon to act as an adhesive agent. The source areas 
were 900 cm?. 


COUNTS PER ENERGY INTERVAL PER I2HRS. 


fe) 190 200 300 400 KEV 


Spectra of background (indicated by broken histogram) and background -++source 
(full histogram) for rhenium (0-62 g in main curve and 0-38 g in inset 


figure). 


Spectrum Analysis 


The results of this more thorough search for activity in 18’Re are 
shown in fig. 2, which covers two energy ranges. The spectra show no 
evidence of activity from the element, as hitherto assumed, and in the 
two energy ranges the counting rates of background and background+- 


source were : 


Wt. of | Range in Background Background Increase 
source kev counts/12 h +source/12 h per 12h 


5 2853-£11: —-185 +227 
0-617 30-350 13038 + 115 128538 +113 + 
0-3808 2 3-3-20 2204 +47 2292 +48 + 88+95 
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The corresponding minimum half-lives calculated from the above 
results for 187Re are 1-3X 1016 yr. and 1-5 x 10° yr. respectively, but the 
complete absence of any trend in the curves leads us to conclude that 
1016 yr. is a conservative estimate of the minimum possible half-life 
unless the energy of the B-rays is exceptionally Jow. Counting rates at 
lower energies confirmed the absence of B-particles down to about 
1-0 key. The occurrence of L x-rays of Re was noted on the lower energy 
curve and reference to this is made later in the paper. If allowance for 
the presence of these L x-rays is made, the above calculated half-lives 
of 187Re will both be approximately 1-4 101° yr. for the two energy 


ranges covered. 


Fig. 3 
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Spectra for osmium (1) Background marked © 2 g 
2 1KG2. ° . ’ ° Background + 1-5 g Os 
indicated b broken histo am a ( . : 
full hist y gr nd (3) background-++-7 5 g Os, shown 


§5. Osmium 


The negative results for rhenium given above led us to examine osmium 
in the form of a chemically stable oxide of purity greater than 99% 

The sources employed contained 1:5 g and 7-5 g of osmium and - 
particularly close search was made for both K and L x-radiation, as it 
seemed not unlikely that decay by electron capture would take lace 
It had previously been verified that no y-radiation, includin annihilate ; 
quanta, was emitted so positron decay was eliminated. ; J 
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Nylon in amyl acetate was used as a holding material for the sources 


and the results of the spectrum analysis in the region of the K x-rays of Re 
are tabulated below. 


Background Background 
Background 8 eae 
l +1-5 g Os per +7-5 g Os per 
per 12 hours 12 hours 12 hours 
1987 +44 2127 +46 2264 +47 


—— 


The spectrum for this K x-ray region, fig. 3, shows no evidence of a 
peak so that we conclude that if K-capture occurs it is not clearly set 
in evidence by these experiments. The gas absorption for the Re K x-rays 
was not high at the gas pressures employed (4 atm of argon) so the 
question of K-capture remains uncertain. In the lower energy region 
positive results were seemingly obtained and fig. 4 illustrates clearly 
the occurrence of x-rays in the region of Re L x-ray energies. Two 
separately observed difference spectra using different sources of the 
same mass gave similar peaks in the energy range 7 to 11 kev and the 
counting rates between the limits 6-0 and 13-3 kev were in good agreement 
over a period of 12 h. 

The origin of the L x-rays was studied further. It was noted that 
a slight general rise in the background was obtained on the insertion 
of osmium into the counter and the possibility of L x-rays arising from 
external conversion of y-rays in the osmium source was investigated. 
The effect of introducing a 20 mg/cm? lining of tungsten powder and a 
sheet of platinum 54:6 mg/cm? are shown in fig. 5 where L x-rays of 
tungsten and platinum are clearly observed. We must, therefore, presume 
that most of L x-rays observed, using the osmium source, are produced 
in the photoelectric interaction in the K and L shells, of y-rays passing 
through the proportional counter, undetected by the G.M. anticoincidence 
system. Accurate quantitative data regarding the comparative number 
of L x-rays produced by osmium and tungsten as sources is lacking 
owing to an unavoidable change of site for the equipment. It would 
seem, however, from the following table, giving the number of L x-rays 
produced by three different elements, that by far the greater fraction 
of L x-rays arising from the osmium must be attributed to this photo-effect. 


Pee ee 


Klement Source of thickness No. of L x-rays 


mg/em? observed per hour 
Os 8:3 38 
W 10-2 43 


Pt 54:6 80 
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It is hoped to re-examine osmium in the near future, comparing the 
emission of L x-rays with those arising from a similar mass of tungsten 


under identical conditions. Je , 

A search for K x-rays from Os was made with a scintillation counter 
heavily shielded from background radiation. The Nal crystal was 
surrounded first with tungsten powder to obtain the background spectrum 


12 HR. 


PER ENERGY INTERVAL PER 


COUNTS 


50 


nd ass Lo 5 20. KEV 


Low energy region of osmium spectra. Histogram 2 represents the background, 
I the background+source of 7-5 g Os and 3 the difference of 1 and 2. 
Histogram 4 is the difference spectrum of another source of 7:5 g 
osmium. The peak due to L x-radiation is clearly shown between 


~11 kev. 


COUNTS PER ENERGY INTERVAL PER 12 HRS 


9 10 20 KEV 


Low energy region of platinum and tungsten spectra. Histogram 1 represents 
the background, 2 the background+ platinum source 54-6 mg/cm? and 
3 the background -+tungsten source 10-2 mg/cm?. 


Fig. 6 
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60 


INTERVAL 


40 
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° 50 100 exe) KEV 


Spectra obtained with scintillation counter surrounded by tungsten and 
osmium. (1) Background represented by full histogram, (2) back- 
ground-+tungsten source marked x, and (3) background--osmium 
source marked ©. 
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and subsequently with a similar mass of osmium >10 g. The results 
of this investigation are shown in fig. 6, from which we conclude that 
if K-capture does occur in 1*’Os its decay probability is very small. 


$6. DiscussIon 


The search for radioactivity in rhenium and osmium shows that the 
minimum value of the half-life, for whichever is unstable, exceeds 1015 yr. 

The very small decay probability that must be associated with either 
suggests that the energy release is small, a conclusion in agreement 
with the energetic studies of Jensen and Suess (1951). It is unlikely 
that there exists a naturally occurring stable isobaric pair and it is 
assumed that the instability associated with either one of these isotopes 
is beyond the limits of detection of our apparatus. 
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CXXIII. Time Variations of Extensive Air Showers and the Origin of 
Cosmic Rays 


By T. E. CransHaw and W. GALBRAITH 
Atomic Energy Research Establishment, Harwell* 


[Received May 28, 1954] 


ABSTRACT 


An experiment is described in which sixteen Geiger counters each of 
area 200 cm? arranged on a square lattice of 54 m pitch were used in 
coincidence to count showers of energies about 101%, 2x 101°, 5x 1018 and 
1017 ev. From the shower rates, no evidence of a deviation from the 
power law spectrum, R(>Z) cc E-1’, was found. The results are analysed 
in solar and sidereal time. In solar time no variations greater than 1% 
for the 1016 ev showers or 10% for the 1017 ev showers were found, and 
these amplitudes could easily arise as a result of statistical fluctuations. 
In sidereal time, no variation greater than 0-6°% was observed at an energy 
of 1016 ev. At 2x10! and 5x101 ev, small variations were observed 
which had about 1% chance of being due to random fluctuations. The 
bearing of these results on theories of the origin of cosmic rays is discussed. 


§ 1. IyTRODUCTION 


EXPERIMENTS have been reported in which the rates of recording cosmic 
ray events as a function of time were studied, and when variations were 
observed, attempts have been made to assign the cause. Two types of 
variation can be distinguished, and these have been called primary and 
secondary variations. Primary variations are variations which have to be 
ascribed to a property of the primary radiation at some distance from 
the earth, whereas secondary variations can be ascribed to local meteoro- 
logical or geomagnetic conditions. At low primary energies (~101° ey) 
such as can be studied by counting neutrons or the total ionizing radiation, 
both types of variation have been observed (Elliot 1952). It appears 
that at these low energies, most of the primary variation can be described 
as a distortion by electromagnetic fields in the solar system of an external 
isotropic distribution, together with the emission of particles by the 
sun. Large increases in this emission have been shown to occur during 
solar flares (Adams 1950, Simpson 1953). The data for these energies 
have been analysed in sidereal time by Elliot and Dolbear (1951), who 
conclude that a variation of about 0-:02% with a maximum at 0500 h 
L.S.T. is probable. 

he eee ee ee ce 

*Communicated by W. J. Whitehouse. 
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At somewhat higher energies (10! ev), effects due to electromagnetic 
fields in the solar system may be expected to be much reduced, and 
primary variations would then be observed in sidereal time, and would 
relate to the radiation outside the solar system. Primary particles of 
energy about 10!*-10!4ev have been studied in attempts to’ observe — 
such effects by counting .-mesons underground. In an. experiment of 
this type Cocconi (1951) found no correlation with either solar or sidereal 
time to within a statistical accuracy of 3% in the rate of arrival of p- 
mesons. Barrett and Eisenberg (1952) and Sherman (1953) also found 
no evidence for anisotropy greater than 2% and 0-5% respectively. 
In a similar experiment, Sekido et al. (1951) reported an effect of 10% 
with a 2%, significance level, the maximum being at 0500 h L.S.T. Other 
observations at this and higher energies have been made by counting 
the extensive air showers produced by interactions of the primaries with air 
nuclei. Hodson (1951) recording showers of estimated energy 5 x 1014 ev 
found a variation in sidereal time of (1:15-+0-61)°% with the maximum 
at 2330 h L.S.T. Daudin and Daudin (1953) found a sidereal variation 
of about 2% with the maximum between 20 and 23 h for showers of 
estimated energy about 10!ev. Citron (1952) however found no 
significant sidereal effect, but a small solar effect of about 1%. Recently, 
Farley and Storey (1954) have reported a small but apparently significant 
sidereal variation for showers of energy probably between 101°—1014 ev. 
When corrected for a seasonal effect, the amplitude of the first harmonic 
of the sidereal variation is 1-5+0-239% with a maximum at 1700 h L.S.T. 

It seems that the only conclusion which can be drawn from these 
results is that the radiation of energy less than 101° ev is isotropic to 
within 1% or 2%. The dissentient results can hardly be considered as 
statistically significant when considered as part of the total number of 
experiments. 

In the present experiment we have extended this result to an estimated 
energy of 101° ev by counting extensive air showers. At an estimated 
energy of 5 x 1016 ev some indication of a sidereal variation is reported, and 
a possible interpretation suggested. The secondary variations observed 
during the course of this work appear to be due almost entirely to changes 
in barometric pressure, and will be described in a later report. ‘ 


§ 2. DESCRIPTION OF THE APPARATUS 


The counter arrangement used in the present experiment is shown in 
fig. 1. ABC... P represent stations, separated by 54 metres. Each 
station contains a single G.M. counter of area 200 cm?, filled with an 
ethyl formate—argon mixture, which is found to eliminate the need for 
temperature control. To obtain the greatest reliability over long 
periods, each counter fed into a quench unit, which, at each count, reduced 
the potential on the wire below the Geiger threshold for a period of 
0-5 msec. The quench pulse also triggered a blocking oscillator circuit, 
which generated a 30 v, lysec pulse in a coaxial cable, terminated 
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by a 70 ohm resistor in the central trailer at Q. Delay lines were inserted 
in the cables from F, G, J and K to make the total delay from all the 
stations the same. 

In the trailer, the pulses from the stations are arranged to cut off the 
current in one of sixteen triodes, whose anodes, connected together, are 
supplied from the cathode of another triode whose grid is maintained at 
+150 v (fig. 2). The current through each of the sixteen triodes is 


Fig. 1 


G.M. Counter arrangement. 


stabilized at 2-5 ma by a large cathode resistor connected to a negative 
supply. Thus each pulse changes the total current by 2-5 ma, and makes 
a pulse of 10 v in a 4K resistor. These pulses are fed into three discrimi- 
nators, with thresholds at 22, 32, and 42 volts. Thus the first is triggered 
by a 3-fold coincidence, and the second and third by 4- and 5-fold 
coincidences. 

The events were recorded by five pens on a moving chart. Two pens 
were used to make a mark each solar and sidereal half hour, and the other 
three pens recorded 3-, 4- and 5-fold coincidences. 

The resolving time of the coincidence circuit was 1 psec, and the total 
rate of all sixteen counters was about 100/sec. The calculated random 
counting rate of 3-fold events was about 0:5% of the true rate and 
was neglected. 

An additional piece of apparatus was a neon hodoscope. This con- 
sisted of twenty-four coincidence circuits, eight of which do not concern 
the present experiments, each connected to a double triode circuit arranged : 
to switch off a small neon lamp when triggered. This hodoscope was 
more reliable in operation than one in which the neon lamps were switched 
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on when a shower was detected. The coincidence circuits were fed from 
the sixteen stations, and also from a master pulse produced whenever the 
3-fold discriminator was triggered. A photograph was then taken 
automatically. In this way, a record of which counters were involved 
in each coincidence was obtained. These records were used to make the 
analysis of 8-fold coincidences referred to later. A report on the study 
of the structure of air showers from these records is also in preparation. 


Fig. 2 


+450 V 


4Ka 


f+ TO DISCRIMINATORS 


Pulse adding circuit. 


A ratemeter and pen recorder was used to provide automatic checking 
of the operation of the counter stations. The input of the ratemeter 
was connected to the signal lead from each of the stations in succession 
for a period of half an hour. Thus each counter was checked three times 
aday. Faults and breakdowns could easily be detected, and the suspect 
part of the record was excluded from the analysis. 

In order to eliminate spurious pulses from the mains, the apparatus 
was run from a motor alternator. At intervals, the apparatus was run 
with the H.T. on the Geiger counters below their threshold. No spurious 
coincidences were observed during these periods. 


§ 3. ANALYSIS OF THE RESULTS 
Calculations based upon different assumptions about the lateral 
structure of air showers, including the Moliére structure function 
(Moliére 1946), were made to estimate the energy of the primary particles 
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producing the showers. All the calculations broadly agree that the 
most probable number of particles in showers which cause 3-fold 
coincidences is about 7105, and the mean number about 10°. This 
number is then used in conjunction with cascade theory (Rossi and 
Greisen 1941) to arrive at the energy value quoted in the table, col. 1. 


SSS.» '' 


‘<) c . - 
Shower Solar Time Sidereal Time 
Events 


Fundamental 2nd harmonic Fundamental 2nd harmonic 


Energy |Rate/|Prob.|Amp.|Phase|Prob.|Amp.|Phase|Prob. Amp.|Phase|Prob.|/Amp.|Phase 
0 0 0 
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f= LO 10-1; 51 0-61) — 30 | 0-82) — | 49 0-63 | 200°} 12 hell |) Tal 
1>2X 101% 3-18] 87 | 0-50) — 10 | 2-0 —SaeOoneo OO) TASS O27 alo:Om PTTL 
Eo lO 1-27) ean, || 2:2 | — 35 | 2-2 == || Ore || ale) || days | took |) aiess )) ee 
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The other energy values have been derived in the same way. The table, 
col. 2, shows the rate at which showers were observed. If we could 
assume that the ‘ sensitive area ’, i.e. the area inside which the cores of 
the showers must fall in order that a coincidence should be recorded, was 
independent of energy, we could deduce directly that the energy spectrum 
in this region was given by R(>H) o H-1:’, R(>E) being the rate of 
incoming primaries with energy greater than H. In fact, the sensitive 
area is likely to be smaller for large showers than for small ones. The 
value 1-7 for the exponent is thus an upper limit. We therefore conclude 
that the usual power law spectrum (Hilberry 1941, Barrett et al. 1952) 
continues at least to 1017 ev. 

The results were also analysed in the following way. The chart from 
the 5 pen recorder was examined, and the number of events of each type 
occurring between the half-hour marks in solar and sidereal time counted. 
Forty-eight columns were then drawn up containing the number of 
counts obtained between 0000 h and 0030 h, between 0030 and 0100h 
and so on through the day. A record of the number of days examined 
for each half hour interval was also taken so that when all the charts 
were analysed, the counts in each half hour interval could be normalized 
+o the same number of running days. A similar analysis was undertaken 
from the hodoscope films for 8-fold events. The results obtained by 
taking the counts in 24 hourly intervals are shown in fig. 3. The heavy 
line is the mean rate, and the dotted lines represent the standard deviation. 
By visual inspection one can see that there are no large variations with 
either solar or sidereal time. In order to put this result more quanti- 
tatively, we have calculated the amplitude of the best sine wave through 
the points and. also the amplitude of the second harmonic. These 
values are given in the table, columns 4, 7, 10, 13 and in columns 3, 6, 
9,12 are given the probabilities that such an amplitude should arise as a 
result of random fluctuations (Chapman and Bartels 1940). 
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It is immediately clear that there are no significant variations in solar 
time. Out of 8 analyses, only one has a significance level greater than 
10%, a result which might easily be obtained if the counts were distributed 
at random. We conclude that the showers producing 3-fold coinci- 
dences in this apparatus are isotropic to better than 1% in solar time, 
the 5-fold coincidences to about 2%, and the 8-fold coincidences to 
about 6%. 


Fig. 3 


SOLAR TIME 


[2 4 6 8 10 12 14 16 18 20 2 24 


Shower Energies: A>10!%ev; B>2x10%ev; C>5x 1016 ev. 


SIDEREAL TIME 


2 4 6 8 10 12 14 16 18 20 22 24 


In the sidereal analysis, we again have no very compelling result. 
However, two out of the eight analyses reach a significance level of about 
1%, and this must be regarded as suggestive. The phases and amplitudes 
of the sine waves have been plotted on a harmonic dial in fig. 4, where 
the dotted circles represent the standard deviation, and the solid circles 
the 10% probability level. On this diagram, the probability of the ‘ true’ 
point being outside the standard deviation circle is 61%, and therefore 
the standard deviation circles are quite likely not to intersect, even if 
the true phase and amplitude are the same in each case. However, it 
can be seen that a constant phase and amplitude of about 1000 h and 
3% would be consistent with the results, but it also seems reasonable 
to suggest a constant phase and an amplitude which increases with the 
energy of the showers. 

The relatively large value of the second harmonic indicates that the 
distribution shows some rather rapid change. This appears to be a 
decrease in the region 20-24 h L.S.T. However, we wish to emphasize, 
that although the statistical uncertainty looks quite small, experience 
shows that it is unwise to rely on any single result. 
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§ 4. Discussion 


We wish now to consider these results in the light of some theories 
of the origin of the cosmic radiation which have been put forward in 
recent years. 

It was suggested (Richtmeyer and Teller 1949, Alfvén 1949) that the 
origin of all the cosmic radiation was the sun and that the particles 
emitted were kept circulating about in the neighbourhood of the solar 
system by the action of an extended magnetic field of strength ~10-5 
gauss. By circulating in this field for period 10°-108 years, the cosmic 
radiation would become isotropic. The upper limit in energy at which 
isotropy is observed is that at which the magnetic field is no longer able 


Fig. 4 


(a) Fundamental. (6) 2nd Harmonic. 
Dotted circles—standard deviation. Full circles— 10% probability. 


to retain the particles in orbits through the solar system. On this theory 
it was expected that no protons in excess of 1014 ev and no other particles 
in excess of 1016 ev would be observed. 

The result of the present experiment shows that there is a remarkable 
degree of isotropy in the primary radiation even up to 10%-10% ev and 
that the power law energy spectrum still holds in this energy region. 
This contradicts the conclusions of the Richtmeyer-Teller theory and. 
supports the hypothesis that the highest energy primary particles are 
not solar in origin. The radius of curvature for a particle of 1017 ev ina 
field of 10-5 gauss is 30 light years (1 lyr =101* cm) which shows that 
the orbits of such particles would be in interstellar space and thus the 
particles would not be retained in the neighbourhood of the sun. In 
other discussions (Cocconi 1951, Barrett et al., loc. cit.) it has also been 
pointed out that, although a solar origin may be ascribed to the lowest 
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energy primaries, the highest energy particles must come from the 
interstellar space of the galaxy as a whole. 

Another theory developed at the same time as Richtmeyer’s and 
Teller’s was that of Fermi (1949). In this, he suggested a mechanism 
whereby the cosmic ray particles could be accelerated as a result of 
collisions with the magnetized gas clouds which wander through inter- 
stellar space. The particle’s energy increases in opposing collisions and 
decreases in overtaking collisions, but since the former are more probable, 
a nett gain of energy results. Provided that the nett gain of energy 
exceeds the energy lost by ionization, the particle energy increases until 
it is finally lost in some catastrophic collision. The sources of the highest 
energy particles (i.e. the sites of the last accelerating collision) would be 
more or less evenly distributed throughout interstellar space. Fermi 
showed that these assumptions led to a power law energy spectrum of 
the primary radiation, but with the rate of increase of energy assumed 
by him, it was not possible to accelerate the heavy ions, owing to their 
large ionization loss. 

The Fermi theory has recently been modified and extended by Morrison, 
Olbert and Rossi (1954). These authors assume a more rapid energy 
gain and a loss of particles to intergalactic space. This modification 
allows the acceleration of heavy nuclei. By applying diffusion theory 
to the scattering by the gas clouds and determining the age velocity 
distribution of the particles, it is shown that a power law spectrum 
is again obtained for both protons and heavy nuclei. To explain the 
isotropy of the low energy primaries it has to be assumed that the average 
distance between scattering centres is only 1 light year, and the theory 
predicts that the power law spectrum should fail above about 10 ev, 
at which energy, anisotropy should become apparent. This again is 
contradicted by the results of the present experiment. 

More information about the magnetic field in interstellar space has 
been given in recent papers (Fermi and Chandrasekhar 1953). It is now 
believed that this field (still of average strength ~10-5 gauss) is directed 
along the spiral arms of the galaxy, and it is supposed that magneto- 
hydrodynamic waves, of wave length ~104 light years and amplitude 
~10 gauss, travel along these arms with a velocity of the order of 10° — 
cm sec”. It has been suggested to us by W. B. Thompson that such waves 
may play an important part in accelerating the cosmic ray particles by 
means of a rather complicated betatron effect, which occurs during the 
time when the magnetic field is increasing, and which accelerates the 
particles between collisions of the type suggested by Fermi. The particles 
would in general be confined within the spiral arm. If this suggestion 
proves to be correct, the time variations of the primary radiation will 
have to be related to our position in the spiral arm, and to the motions 
of the particles within it. Speculations on this point would lead us to 
look f: or anisotropy of the high energy radiation connected with (1) a drift 
of particles along the arm, (2) a drift outwards along the radii, (3) motion 
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round the arm. Such speculations will have little weight until both 
the theory and the observations are more definite. Even before this 
suggestion had been made, it had, however, occurred to us that the 
minimum indicated at 20-24h L.S.T. takes place at a time when the 
apparatus is facing inwards along the spiral arm in which the solar system 
is situated. 

It is evident from this short discussion that the results of the experi- 
ment have a considerable bearing on our knowledge of the origin of cosmic 
rays, since each of the current groups of theories (i.e. those based on a 
solar origin and those based on acceleration in interstellar space) would 
have to be modified to explain the apparent isotropy and power law energy 
spectrum in the region 1016-10!’ ey. 

It is unfortunate that observations of showers at sea level cannot, as 
far as we know, give any information about the nuclear charge of the 
primary particle. This can only be studied in nuclear emulsions at 
great altitude and the area of the apparatus is then so small that statis- 
tically valid measurements are virtually impossible for energies above 
1015 ev. However, further knowledge of the distribution in energy 
(i.e. the power law) of the highest energy primary particles and of their 
distribution in direction, independent of the nuclear charge, will be 
valuable. This can only be obtained by increasing, very greatly, the 
area in which extensive showers can be observed. We are therefore 
planning a larger experiment of the type described here, to test the power 
law and search for time variations in the region 1017-101% ev. In this 
connection it is interesting to note that the radius of curvature of a 
proton with an energy of ~1018 ev in a field of 10-5 gauss, is 3 10° cm 
which is nearly equal to the supposed radius of the spiral arm. We shall 
therefore be approaching a region of energy in which we might expect 
some change in the behaviour of the radiation, quite independent of the 
details of any theory of the acceleration of the particles. 
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ABSTRACT 

This paper is a sequel to a previous one (Read 1954) dealing with the 
acceptor centres associated with dislocations in germanium. The previous 
paper developed two approximations for the statistics of occupation of 
dislocation acceptors. This paper presents some further approximations 
that take account of the non-uniform spacing of accepted electrons along 
the dislocation. It also treats the case where the dislocation acceptors 
have a range of energy levels (rather than a single level). The several 
approximations are compared by plotting the fraction of acceptors that 
are full vs. temperature for a particular example. 


§ 1. INTRODUCTION 

Tus is the second in a series of papers on the electrical effects of dis- 
locations in germanium. The first paper (Read 1954) discussed the 
observed acceptor-type behaviour of dislocations in germanium in terms of 
the dangling bonds on the extra plane of a dislocation having some edge 
component. In the present paper it is sufficient to regard a dislocation in 
germanium as a row of acceptor centres, or sites, where extra electrons can 
be added. For dislocations in or near the edge orientation, the spacing 
between sites is the order of a few Angstroms. The energy level of these 
dislocation acceptors appears to be slightly above the middle of the gap— 
see Pearson, Read and Morin (1954). 

In n-type germanium with the Fermi level above the dislocation acceptor 
level, the dislocation accepts electrons from the conduction band and 
becomes negatively charged. Surrounding the negatively charged dis- 
location is a cylindrical region of fixed positive space charge. ‘The first 
paper derived the electrostatic energy as a function of the fraction f of 
sites that have accepted electrons and acquired a negative charge. Two 
statistical approximations were obtained for the temperature variation of 
f, and f vs. T plots were given for a particular example. 

The present paper forms a sequel to the first one and is a further study 
of the statistics of occupation of dislocation acceptors in n-type germanium. 
The problem is peculiar in that the added electrons on a dislocation are 
so closely spaced. Because of the strong coulomb forces between neigh- 
bouring electrons, the electrons in dislocation energy levels form a one- 
dimensional lattice. 
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In the first paper the statistics of occupation of dislocation acceptors 
was treated on the basis of two approximations. In one, the minimum 
energy approximation, the accepted electrons were assumed to be 
uniformly spaced along the dislocation. In the other, which was based on 
Fermi statistics, all ways of arranging the electrons in sites along the 
dislocation were assumed to have the same energy. The present paper 
gives some further approximations that take account of the effect on 
energy of the non-uniform spacing of the accepted electrons on the 
dislocation. 


§ 2. SmncLe Enercy Levet—Non-UnIFORM SPACING 

As in the first paper we consider the case where the line of acceptor 
states have a single electronic energy level &,. (Later in this paper we 
consider the case where the dislocation acceptors have a range of energy.) 
Let &, be measured upward from the top of the valence band at the dis- 
location. (The energy level &, of the top of the valence band will vary 
in the space charge region around the dislocation.) 

At the absolute zero of temperature the accepted electrons, or full sites, 
are uniformly spaced along the dislocation at intervals of a=c/f where f is 
the fraction of sites that are full and c is the spacing between sites. As — 
the temperature is raised above 7'=0°K some of the electrons on the 
dislocation acquire enough thermal energy to jump into neighbouring 
sites; so the spacing between electrons becomes non-uniform. Let 
5&, be the added electrostatic energy of non-uniformity and let S be the 
entropy. As in the first paper we let &, be the electrostatic energy for 
uniformly spaced electrons. Again measure the Fermi level &; upward 
from the valence band in the normal n-type material. (See fig. 2 of the 
first paper.) The free energy F per electron measured relative to the free 
energy of an electron in the normal n-type material is 

F=6,—6;+ €,+86,—TS. . . .. . (2.1) 

We shall define 

6F=0€,—TS<0 2 tee 


as the added free energy of non-uniformity. To find the f vs. 7’ relation 
we minimize the free energy fF per site with respect tof. This gives 
é*(f)= 6y(T)— €,—85F* (f,T) . . . . . (2.3) 
ses é*(f)=(d/df) fé , (f) was found in the first paper. In this paper we 
n 
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§ 3. ENERGY or Non-Untrormiry 8 6, 
Let x; be the displacement of the ith electron from the position it has 
when the electrons are uniformly spaced, 6&, is a function of all the 2;,’s. 


Since the electrons are in acceptor centres 2; must be an integral multiple 
of the spacing c between sites. 
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In this section we find 6é, using an Einstein approximation, in which 
each electron is assumed to move in the average field of the other electrons. 
Then the energy of each electron is a function only of its own displacement. 
In § 5 we discuss in some detail the conditions under which this approxi- 
mation is good. r 

Since the field of the fixed positive space charge vanishes along the 
dislocation, the energy of non-uniformity comes entirely from the forces 
between electrons. An electron that has a displacement 2 will be acted on 
by a field H(~) due to the other electrons and will have an, added energy 


7 
q| B(x) de | 
, Isa 
where q is the electronic charge and 
n= 6 > 
ta ne | (hoa) ae 
2 2 
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where x is the dielectric constant. The linear stiffness coefficient K is a 
function of f 4g2_ 1 4808 
d pee) Se Me ae 0 £3 
k= s—-—_a ee F (3.2) 


&,=@/kc is the energy of interaction of electrons in neighbouring sites. 

We shall deal with the case where the displacement x is small compared 
to a, so that the linear term is a good approximation. 

The added electrostatic energy is }Ka?. Thus if the displacement is 
nc (where n must be an integer), the added electrostatic energy is 
86,=4K c'n?=(2:40 &)f%) n?. 36, is the average of 8&,,; so 

SOW CARI Be Ee oe 
where n2 is determined from statistics as will be discussed in the following 
section. 

§ 4, STATISTICS 


In this section we consider the statistics of the system of electrons 
spaced non-uniformly in sites along the dislocation. For a given f the 
statistics of the system is specified by the distribution of electrons among 
the values of »; that is, we divide the electrons into groups where all 
those in the nth group have added energy ¢,,=2-40 6, f*n?. Thus the 
system is one to which Boltzmann statistics applies. The number of 
electrons having displacement nc is proportional to exp (— &,,/kT) where 
kis Boltzmann’s constant. It will be convenient to define a parameter 

2-40 & of? 
The probability of finding an electron with displacement nc is proportional 
to exp (—zn?). We shall express the final results as a function of the 
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parameter exp (—z), which is a measure of the non ne It is bee 
ratio of the number of electrons having displacement t= +c to the num x3 
having displacement c=0. For the uniform distribution exp (—z)=0 an 
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Variation of decrease —6F in free energy and increase 5&, in electrostatic energy 
with non-uniformity, which increases with exp (—z), where z=2°40 & 9 f?/kT. 


for a completely random distribution exp (—z)=1. However, the present 
approximation breaks down before we reach a completely random 
distribution. (We discuss this in detail in § 5.) 

The added free energy SF is 


6F=—kT In J exp(—zn?) . . . . . « (42) 


where the sum is taken from — oo to oo. (See J. C. Slater, 1939, Intro- 


duction to Chemical Physics, p. 50.) Figure 1 shows a plot of —SF/kT vs. z. 
For exp (—z)>0-4, SF =—(kT/2) n (77/2). 
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From eqns. (3.3) and (4.1) 


ee a ee ete | ee (4.8) 
where 

= d 

n= — zn 2'exp (—2n’). esebesh. cake » (454) 
Thus 


d In &' exp (—zn?) 
dlnz . 
Figure 1 also shows a plot of 86,/kT'vs.z. For exp (—z)>0-4, 
06 =kT/2. 
From eqns. (4.1), (4.2) and (4.5) 
O5F d06F dlInz 
Gili fameo line, on 7, 


We now have for 6/™* in eqn. (2.4) 
Steer eo Cpe eee. sf  .  n(4-7) 


So 5F* can be found from fig. 1 for any values of 7 and f. Thus eqn. 
(2.3) can be solved for fvs. 7. The first paper gave &*(f) and &,(7’) for 
the case of a particular dislocation in 1:7 ohm cm n-type germanium where 
c—44A4. The donor and acceptor concentrations were N g=2N ,=2 x 10"° 
cm-3,- The dislocation acceptors were taken as lying 0-225 ev below the 
conduction band. Figure 2 shows the two approximate f vs. 7’ relations 
derived for this example in the first paper. The approximations are the 
curves marked M.E. (Minimum Energy) and F. (Fermi). The solid 
curve marked B. (Boltzmann) in fig. 2 is the f vs. 7 relation as given by 
eqns. (2.3) and (4.6) of this paper. We shall call this the Boltzmann 
approximation. As expected this more accurate approximation lies 
between the other two, which were based on upper and lower limits for the 
free energy. 


Secs Me 
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§ 5, LIMITATIONS AND REFINEMENTS 


This section discusses how the assumptions made in deriving the 
Boltzmann approximation limit the range of its validity. Some refine- 
ments are considered which extend the range of validity and indicate the 
accuracy. 

Correlated Displacements 

The Einstein model gives the correct average energy of non-uniformity 
per electron 5&, when the displacements are completely uncorrelated ; 
that is, when the probability that an electron will have a given displacement 
is independent of the displacements of the other electrons. At low 
temperatures and high values of f the energy per electron 56&, is small 
compared to the energy of interaction of neighbouring electrons ; so the 
strong coulomb forces between electrons cause the displacements to be 
highly correlated. That is, if the displacements were analyzed into a 
Fourier series, the long wavelengths would account for a major part of the 
deviation from uniformity. Thus the Einstein model overestimates the 
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energy 5é, for a given 72. However, this error does not seriously affect 
the f vs. 7’ curve as we now show by using another approximation—one 
that underestimates 5 & ,. 

Suppose we consider only the interactions between nearest neighbours, 
which are spaced at intervals a in the uniform distribution. Now 
define a, as the change in distance between the ith and i+ 1st electrons. 


Fig. 2 


Variation of the fraction f of dislocation acceptors that are full with 7'=absolute 
Uae ee The approximations are: B.=—Boltzmann, F.—Fermi 
-E.=Minimum Energy, and H.O.=Harmonic Oscillator. 


Again x,=nc where n is an integer. If we require that x, be small com- 
pared to a and take account of the fact that positive and negative dis- 
placements are equally probable, then we find that §¢,= Ef? n2 which 
differs from eqn. (3.3) by a factor of 2-4. This gives a lower limite to 
8é » because it neglects all interactions except those between nearest 
neighbours. The curve marked B’. in fig. 2 was obtained in the same way 
as the one marked B. except that the relation 8¢.= & f3n2 was used 
instead of eqn. (3.3). The two curves differ only sient | 
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Non-linearity 

As the root mean square displacement increases, the non-linear terms in 
the electric field become important. These are not as serious as might be 
thought, however. For example, from eqn. (3.1), we see that the 
cubic term is only about 10% of the linear term even for v=a/4. Neglect 
of the non-linear terms leads to underestimating 5 &, and, therefore. tends 
to compensate for the effect of assuming the displacements are uncorre- 
lated. 

Lower Limit for Free Energy 


As discussed in the first paper Fermi statistics gives a lower limit for the 
free energy. This corresponds to taking 5 & ,=0 and giving S its maximum 
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Decrease —8F in free energy vs. f at T=300°K. The two dashed curves 
represent Fermi and Boltzmann statistics, respectively, which are 
limiting forms. 


value, which is found from the number of possible ways of filling a fraction 
f of the acceptors. However, if eqn. (4.2) is extended to sufficiently 
low f it gives values of 8F lower than the lower limit given by Fermi 
statistics. Physically this means that the x; become so large that the 
electrons are no longer confined to discrete intervals on the dislocation ; 
consequently some of the electron distributions that are counted as different 
in Boltzmann statistics are, in fact, indistinguishable. This error will be 
most serious at higher temperatures. In the remainder of the section we 
find f at 7=300°K by a further approximation that combines Boltzmann 
and Fermi statistics. 

Figure 3 is a plot of —dF vs. f at T=300°K. The two dashed curves are 
from Boltzmann and Fermi statistics respectively. In Fermi statistics 


1126 W. T. Read, Jr. on 


—$F=+7S where S is uniquely determined by f. The Boltzmann curve 
depends on the parameter & in eqn. (4.1). For the particular example 
we have considered (where c=4A and k=16), & =0-225 ev. At 
high values of f the root mean square displacement is small compared. 
to a, and the Boltzmann curve becomes an increasingly good approxi- 
mation. At low, SF approaches the lower limit given by Fermi statistics. 
The solid curve in fig. 3 approaches the Boltzmann curve at high f and the 
Fermi curve at low f and makes a smooth transition between the two 
limiting forms. The point B.-F. in fig. 2 was found by taking 6#'* from 
the solid curve in fig. 3. This corrected point lies only slightly below the 
uncorrected Boltzmann curve in fig. 2. The reason for the small change 
is that errors in SF and 50F/d Inf tend to compensate. The final result 
is relatively independent of the exact form of the solid, or transition, curve 
in fig. 3. 
Uncertainty in &* at low f 

The above uncertainties are peculiar to the Einstein model and resulting 
Boltzmann statistics. At low values of f all the approximations to the 
fvs. T curve are subject to the uncertainty in &* at low f. As explained in 
§ 7 of the first paper, there is no definite space charge unless the spacing a 
between accepted electrons is smaller than both & and the mean spacing 
L=(N,—N ,)-1? between excess donors. An even more stringent require- 
ment comes from the assumption, made in deriving &,(f) and &*(f), that 
the positive space charge is uniformly distributed. This will become 
exact as the radius F of the space charge cylinder becomes large compared 
to the mean spacing between excess donors. As shown in the first paper 
eqn. (7.5), 

R fL 

Lo Nae’ 

In the example of fig. 2, L=108 A. Thus even at the maximum f, R is 
only about 3L. At f=0-0126,R=L. This uncertainty in the formula for 
é*(f) may be as serious as any other single source of error in the final 
f vs. T curve. 


(5.1) 


§ 6. THe Harmonic OSCILLATOR APPROXIMATION 


In the previous sections we have considered the case where the dis- 
location acceptors have a single energy level &. This will certainly be 
the case when c is large enough, that is, when the dislocation is near 
enough to the screw orientation. As the angle between the dislocation and 
its slip vector increases, and the dislocation moves toward the edge 
orientation, c decreases and approaches a minimum of 0-886 b where b is 
the interatomic spacing in a (110) direction (see § 2 of the first paper). 
In germanium b=44. As c decreases, the wave functions of adjacent 
dangling bonds overlap more. Thus what, at large c, was a single energy 
level may become a range of levels, or energy band. 

Let &, be the minimum energy of an electron in a dislocation acceptor 
and let &, be the band width. Then an accepted electron can be acceler- 
ated and can acquire a maximum kinetic energy &, while still remaining 
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in the dislocation acceptors. The energy & of an accepted electron will be 
a function of its momentum P along the dislocation and will be a minimum 
&é, at P=0 and a maximum at P=h/2c. We shall use the approximation 
&=6,+ &,[1—cos (27cP/h)]. The equivalent mass m* for an electron 
whose kinetic energy is small compared to &, is given by 


1 Zé 47°? € , 
=—(i),_,= ip Pe eter 8 (6.1) 


where / is Plank’s constant (=4-14 x 10-1 ev sec). 

Again we use an Einstein model, where each electron moves in the 
average field of the other electrons and the maximum displacement is 
small enough that the restoring force is approximately linear. Then each 
electron oscillates with frequency v where (27v)?=K/m*. From eqns. (3.2) 
and (6.1) we have for the quantum hy of vibrational energy 


hv=f2/24/(480 &, &,). pen ee oe ee (G:2) 
Each electron is thus a harmonic oscillator having one degree of 


freedom. The average number 7 of energy quanta per electron above the 
zero point energy 4 hv is 


1 
(See any treatment of the Einstein theory of specific heat of a solid.) 

The average vibrational energy per electron (%-+-4) hv is half kinetic and 
half electrostatic (in contrast to the previous, static, approximations, in 
which the energy of non-uniformity was assumed to be entirely static). 
The. entropy per electron is S=k [(1-+-7) In (1+7%)—% In 7] and the free 
energy of non-uniformity is 


SF =(i+4) by—T 8 (i) 


hy hy 
=" 407 n| 1—exp(— zp) |. foe aes wr (Ord) 


It is convenient to regard 5F as an explicit function of # and v. 
Equation (6.3) comes from (0/d”) dF (nm, v)=0. If we use this and 
remember that v is a function of f, we have from (2.4) 


0 dF av 
6F*=8F +f ay Of 
06F dlnv 
ee 6.5 
pe wee Een ee?) 
From (6.4) (0 5F/dv)=(a-+4) h and from (6.2), (dInv/aInf)=3. Thus 
spre (i+ 3) SST aln (1s eee ae 28 (6.8) 


where (6.3) gives 7 vs. (hv/kT') and (6.2) gives hv as a function of f. From 
(6.6) and (2.3) we can find the f vs. 7’ curve. 

The harmonic oscillator approximation involves an additional parameter, 
m* or &,, not present in the previous, static, approximations. There is 
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also another limitation on the range of validity : the kinetic energy 
8€,=1 hv (Xn +4) must be small compared to &, otherwise the electron 
will not stay on the parabolic part of the & vs. P curve and the average 
kinetic energy will not be equal to the average potential energy. If 
8 &, comes out larger than the band width, then we can conclude that the 
oscillator approximation is a poor one and that the energy of non- 
uniformity is mainly static. When 5&,<@&, is small compared to 6 6 55 
then a static approximation will be valid. 

We now return to the example which we used to illustrate the f vs. T 
curves (fig. 2) for the static approximations. The curve marked H.O. in 
fig. 2 is the oscillator approximation for the same example with &,—0-10ev 
The equivalent mass m* is then 4:75 times the mass of a free electron. 
Equation, (6.2) becomes hv=0-33 f?/2._ The harmonic oscillator curves for 
f vs. T are subject to the same uncertainties discussed in §5. Again the 
assumption of uncorrelated displacements introduces errors at low 7’, 
which are relatively unimportant since the energy of non-uniformity has a 
small effect at low 7’. We could get another limiting approximation by 
considering only nearest neighbour interactions. We would then have a 
Debye model rather than an Einstein model. A still more accurate 
approach would be to treat the electrons as a one-dimensional plasma and 
investigate the oscillations. At higher temperatures errors due to non- 
linear restoring force come in and tend to cancel the errors due to correlated 
displacements. At still higher temperatures the displacements become 
large enough that violations of the exclusion principle are significant. 
However, for the H.O. curve of fig. 2, the free energy drop —8F* is well 
below the limiting value given by Fermi statistics even at 7=300°. A con- 
struction like fig. 2 had little effect on the 300° point. On the other 
hand the uncertainty in &*(f) is greater than for the Boltzmann approxi- 
mation because f is lower; this may well be the most important single 
source of error at the high temperature end of the curve. 

The requirement that 5 &, be small compared to &, is satisfied through- 
out the temperature range. 
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ABSTRACT 

A restricted proof is given of the third law of thermodynamics. The 
restriction resides in the fact that the proof applies only in the limit of 
an infinitely large system (the ‘limit L’). An additional restriction is that 
it applies only to non-interacting and structureless quantum mechanical 
particles whose quantum states satisfy certain requirements (the * class O ’ 
of systems). The essence of the proof depends on showing that the free 
energy of the system, as averaged over a grand canonical ensemble and 
expressed per unit volume of the system, is of the form a—bT"(b>0, 
d>1) near the absolute zero. 


§ 1. INTRODUCTION 


‘An important task of statistical mechanics is to furnish an explanation 
of thermodynamics. In particular, the problem of establishing the second 
law, by proving the so-called H-theorem, has been attacked by several 
writers. The first important achievement to be based on quantum 
mechanics was due to Pauli (1928), who utilized the method of transition 
probabilities. His result is therefore valid only for restricted time 
intervals. An assumption of random a priori phases, to replace the 
classical assumption of molecular chaos was sufficient to yield the H- 
theorem in the form dH/dt<0 for these restricted time intervals, H being 
a certain quantity characteristic of the ensemble under consideration. 
An extension to arbitrary time intervals is due to Klein (1931) and 
Tolman (1938). It was shown later that, subject to the above and 
certain additional assumptions, the ensemble average of a quantity H 
at a certain time, and the long time average of a quantity H for the 
system of interest, are almost always very nearly the same (von Neumann 
1929, Pauli and Fierz 1937). 

No quantum statistical proof is, however, available of the third law. 
The problem was recently taken up by Simon (1951). He considered 
under what conditions it is possible to infer Nernst’s theorem from the 
experimental fact that the specific heat of physical systems vanishes as 


ewes eee SSS 
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the absolute zero is approached. He points out that a system for which 
this inference is not valid must have rather curious properties, and asks 
if one needs a supplementary postulate to express the fact that they do 
not occur in nature. ‘The answer seems to us to be in the affirmative, 
since even a system whose free energy satisfie: the simple relation 
F=a(z...)—b(x...)T near T=0, 21, 2... being external parameters 
of the system, has a specific heat which vanishes near 7'=0, but does not 
necessarily satisfy Nernst’s theorem. Simon pointed out that the 
condition pV=sU where s is a constant and U the internal energy of 
the gas, is a sufficient supplementary condition. This supplementary 
condition has the merit that it follows from the quantum statistics of 
ideal gases, but it still leaves the vanishing of the specific heat to be 
explained. 

We shall show here that the position is somewhat more favourable. 
For a careful formulation of the quantum statistics of the ideal gas 
furnishes in fact an expression for the free energy from which both the 
third law and the vanishing of the specific heat can be deduced. By 
restricting attention to non-interacting particles without internal degrees 
of freedom, we are, of course, dealing only with the theoretically simplest 
case. A generalization of the present argument to imperfect gases would 
be desirable. 

The most important difficulty which has to be overcome is that one 
cannot exclude on quantum mechanical grounds the possibility of 
degeneracy of the lowest energy level of a system (Ludloff 1931, Schottky | 
1943). This suggests that the third law in Planck’s form (for our purposes 
this may be stated in the form ‘ the entropy of a system tends to zero as 
the absolute temperature tends to zero ’) may not be rigorously true for 
finite systems, but that it may be a law which attains validity only in 
certain limiting cases. 

The case with which we shall deal here is that of a system of non- 
interacting particles in the ‘limit ZL’. The limit L is the limit in which 
the total number of particles NV, as averaged over a grand canonical 
ensemble, and the volume V of each system of the ensemble, are each 
allowed to tend to infinity so as to keep the ratio N/V a finite and non-zero 
constant, which will be denoted by p. Under these conditions the 
distinction between the coarse-grained and the fine-grained probabilities, 
which plays such an important part in the proofs of the H-theorem, can 
be neglected, since it cannot affect the results. The following definition 
will be required. A system of class O is one which consists of structureless 
non-interacting quantum mechanical particles, such that each satisfies 
the following requirements : (i) there is an infinity of particle quantum 
states whose particle energies are H,, H,...; the labelling of the states 
is arranged so that H;<E,,, for all j ; (ii) EH, is finite if 7 is finite ; (iii) 
all energy levels with finite j behave like E,;=E,+6,V~° for large volumes 
V, where 6; is independent of the volume, and E, is the value to which 
all those H,’s tend in the limit L for which j is finite; (iv) for large 
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volumes there exists a finite positive integer M, such that a continuous 
spectrum approximation is valid for j>M, and the jth quantum state is 
given by £,;=H,+cV~°%j" for j>M. M must exceed the degeneracy 
w, of the lowest particle energy H, of the system of interest when its 
volume is finite. Under (iv) one might equivalently specify the density 


of states function 
7 1 Ve ae 1/7-1 
NO Nts aa (emg) et, 


In all cases of interest 6=7 (Landsberg 1954 a), and this will be assumed. 

We need some assumption which will enable us to extrapolate physical 
laws to temperatures which are infinitely close to the absolute zero. 
The assumption which we shall adopt, and which will be referred to as 
assumption (A), is “the relations for ensemble averages near T=0 and in 
the limit L, are assumed to remain valid as T tends to zero”’. 

All our conclusions concerning conditions at the absolute zero will 
be deduced from the following result: For a system of class O the free 
energy per unit volume (as averaged over a grand canonical, or a canonical, 
ensemble) is, at temperatures near the absolute zero and in the limit L, 
given by 

f=a—bT* (6>0, d>1). eae Aer ie ( 1.) 


A proof of this result is given in §2. Basically, it depends on certain 
simple theorems which hold for systems of class O in the limit L 
(Landsberg 1954c). 

The relation between the third law and Carathéodory’s (macroscopic) 
formulation of thermodynamics is an additional problem. It does not 
seem to have been discussed so far, and will be the subject of a sub- 
sequent communication. 


§ 2. PROoF OF (1) 
We shall average over a grand canonical ensemble, and define the 
ensemble average of the entropy per system by 
ee ee) log P(N, esq es); 


Ny Ne 
P(t, ta, - - -) = eXP {Q4+2(u—BE,)n,}, B= VT. 
a 
Q and p are independent of the suffix 7, and n, is the number of particles 
in the ith quantum state. The n,-summation is over all permitted 
occupation numbers and the i-summation is over all particle quantum 
states of a system. Let U be the mean energy per system. Then the 
mean free energy per system is given by 


Fa 0—TS=HpkTNFEPE log (1+ exp (u—n)], m= Hk, (2) 


where 


where the top signs apply to Fermi—Dirac systems, and the bottom signs 
to Bose -Einstein systems. We must distinguish four cases. 


4E2 
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(a) FermiDirac systems. For any finite integer in 
1 /kT\ Ut @ eWt-1 dx a 
OMe aa 
where ~ denotes equality in the limit Z. The first relation is due to 
the fact that 7;/V ~ 0 when j is finite, where n, is the mean occupation 
number of the jth quantum state. The last expression is obtained by 
taking m=WM and using statement (iii) in the definition of a system of 
class O. The average free energy per unit volume of a system of the 
ensemble, f= F/V, satisfies similarly by (2) 
kT 1/7 oO gilt dx ( 
CL Pree (aoe ————— re - 
matte (GF) | ope ee 
In virtue of assumption (A) the integral in (3) must diverge like Ta 
asT +0. It follows that ~—y) > « as T' > 0, so that the approximation 
fore) s +1 2 1 
eg [Os]. s>—l, (5) 
9 l+exp(#—a) s+1 a 
can be used near 7’=0. By combining (3) and (5) one can obtain an 
approximate expression for » near 7=0. This expression and eqn. (5) 
are then used to eliminate both the integral and yx from (4). The result is 


T 


1 2 6 
p~pe fexp (u—7;) +1} 74 ~ 
j=m 


cpt = ep*t* /axkT\2 F 
f ~~ pHot 7+1 = om. Cp? ; ( ) 
(b) Bose-Einstein system with 7>1. Our basic equations are now 
1D /(dePN Sieh gil Ske 
(i 
tLe) Jo exp (@—p-+%)—1 
kKTNis f° eu da: 
~ pkT p—kT | — | —__—__———, 8 
sagt yet ( = 0 Exp (T—n+%)—1 (8) 


since such a system does not show the phenomenon of condensation 
(Landsberg 1954 a). We need the results 


Eee 

{, axp wal Base alee 8220, us ey 
xs da 
if Seal =I\(s+1)F'(—s)(1—e*)*,, 0<—a<1, —l<s<0. (10) 


I(x), ¢(v) are the usual Gamma and Riemann Zeta functions respectively. | 
(10) is due to de Groot et al. (1950). In virtue of assumption (A) the 
integral in (7) must diverge like 7-1/7 as 7'— 0, so that, by (10), > Np 
as T’+>0. By combining (7) and (10) one can obtain an approximate 
expression for « near 7=0. This expression and eqn. (9) are then used 
to eliminate both the integral and yz from (8). The result is 


Fe pHy—P(U+ 1/7) C+ 1/r)e-V(eT 1H, | 
(c) Bose-Einstein systems with t=1. The relations (7) and (8) hold 
again, but the integral in (7) can in this case be evaluated in closed 


form, and » is now found to depend exponentially on the temperature. 
The final result is again given by (11) (with 7=1) 
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(d) Bose-Einstein systems with r<1. Systems of this type are below 
their condensation temperature near 7'=0, so that eqn. (7), which 
determines p,, is replaced by 

[eet eae RPP ea poets (LZ) 
in this case. The equation for f is again given by (8), but is now subject 
to (12). It therefore reduces immediately to (11). 

This establishes the result (1) for the case when one averages over a 
grand canonical ensemble of systems of class O for which the mean total 
number of particles per unit volume of the system of interest is p= N/V. 
By averaging over a canonical ensemble of systems of class O for which 
the fixed number of particles per unit volume of the system of interest, 
N/V, is also p, one obtains the same result in the limit L, since the distri- 
bution functions are then the same (Landsberg 1954b). This proves 
assertion (1). 

For completeness we give below the values of » near 7’=0: 


| Zo Cp? 


LT zi LP (F.D.) 

E P(ijr) Pa TAN) (Ee D\ ae) 

ia a] (=) (B.E., 7>1) 

f= \ 2 ; (13) 

ap exP (- “7 (B.E., 7=1) 

Ey 

a7 (BE, 7< 1) 
aE 


§ 3. Discussion 


It follows from (1) that for a system of class O in the limit L the mean 
total energy per unit volume (wu), the mean entropy per unit volume (s), 
and the mean heat capacity at constant volume expressed per unit volume 
(cy) are given by 


u=a+(d—1)bT", s=bdT@-1, ¢,=d(d—1bTt. . (14) 


Thus both entropy and specific heat tend to zero, and, in this sense, the 
third law is established for these cases. 

The principle of the unattainability of the absolute zero may be formu- 
lated as follows. Let 21, 7,... be the external parameters of a physical 
system, and let 


[G2F/OT Axles ayo 


es ee 5 (le 
sf Obs 41> Vo ee j= fozr/eT\. ae, se 0 ( ) 
Then the absolute zero is unattainable if 
lim B(T, 21,...)=0 Ria chou a et ee 


T—0 
for all physical systems, and whichever external parameter is chosen to be 
v,. It is easy to verify the accuracy of this statement. or instance, 
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Bennewitz (1926)* showed that the quantity (T/c,)(Op/eT), must tend 
to zero as the temperature tends to the absolute zero, and our statement 
(16) is equivalent to a generalization of this. We can infer from (1), as an 
approximation, that for a large but finite system of classO F=aV—bVT4, 
where V is the volume of the system. It follows that (one may take 2, to 
represent the volume V) 


T [do log (bv) 
Joyidl: U1, X%9,-- Faas Ber) eel Ae . ee ee et" (17) 


Tf all variables are continuous functions of the external parameters, it 
follows that the systems of class O satisfy (16), and therefore the un- 
attainability principle. (The general result that Planck’s form of the 
third law implies the unattainability principle could not be used in the 
present discussion, since the results (14) hold in the limit L, while a finite 
system is contemplated in (16).) 

It might be thought that a possible formulation of the third law is to 
stipulate the unattainability principle and the condition that 
(02F/0T?),,,... OF ¢,/7' remain finite as the absolute zero is approached. 
We see from our results that this form of the third law is violated by those 
Bose-Einstein systems of class O whose density of states for high quantum 
numbers is proportional to an inverse power of the energy (7>1). Such 
systems, which include the perfect Bose—Einstein gas in a one-dimensional 
box, are of no practical importance at the present time, but a satisfactory 
formulation of the third law should presumably apply to them also, 
particularly so, since they do satisfy the unattainability principle as well 
as Planck’s form of the third law. A satisfactory formulation of the third 
law should, therefore, not lay down that 0?F/dT? or c,/T be finite at the 
absolute zero of temperature. To cover the case of the perfect gases it is 
sufficient to require that (1) hold. Nernst’s theorem in Planck’s form, 
the unattainability principle and the vanishing of the specific heat at low 
temperatures all follow from this result. 


APPENDIX 


A thorough investigation of the relationships between the various 
formulations of the third law is due to Bennewitz (1926). We merely 
wish to put on record here that three results which he states (p. 155) may 
be disproved by counter examples. The results are : 

(i) If U=F at T=0, then U and F are finite at T=0. 
(iil) If(OU/eT),.. .=0 at T=0, then (OF /OT),,... is finite at T=0. 

(ii) If (OF /0T),,...=(0U/0T),,... at T=0, and F is finite at fe) 
then (OF /0T),, ...=0 at T=0. 


eee esate Se 
* Some remarks on Bennewitz’s article are made in the Appendix. 
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The counter examples (numbered to correspond) are : 
(i) If F=a+0b/T, then U=F—T(OF/0T),, ... =a+-20/T. 

(ii) If (OF /0T),, ...=—a log (—log dT), then 

(QU /OT),,, .. .=—a/log bT(b>0). 

(iii) If F=a+6T” log T, then U=a--bT? [(1—p) log T—1], 0<p<l. 
It follows that 0F/0T and dU/0T tend to —co as the absolute zero is 
approached. 

In these examples a, 6 are non-zero, finite and continuous functions of 
the external parameters only, and all equations are supposed to represent 
conditions near the absolute zero of temperature. Although Bennewitz 
does not deal with this point explicitly, we believe we follow him in making 
the convention that if two quantities diverge, both to + 00 or both to — 0, 


at a certain temperature 7’, the equation A=B is regarded as applicable 
at the temperature 7’. 
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SUMMARY 


Attempts have been made to interpret the observed galactic radio 
emission in terms of the integrated radiation from sources having a 
spatial distribution similar to that of the common optical stars. It 
is found, however, that such models fail to explain the large intensity 
observed towards the direction of the galactic poles. The origin of 
this additional component is discussed and the two most probable 
explanations seem to be (a) the incidence of an isotropic component from 
extra-galactic sources, and (b) the existence of a new class of object having 
a wide spatial distribution relative to the galactic plane. 

The integrated radiation from normal extra-galactic nebulae is shown 
to be inadequate to account for the observations. An estimate has been 
made of the intensity likely to be produced by sources of the type 
associated with the collision of two galaxies, and it is shown that such 
sources will increase significantly the emission per unit volume. An 
independent figure can be derived from data on the large scale irregularities 
in the distribution of galaxies. 

The total intensity to be expected from extra-galactic sources is found 
to depend strongly on the more distant regions of the Universe, and the 
value obtained, therefore, depends on the Cosmological theory adopted. 
Accurate measurements of the extra-galactic component might thus be 
important in distinguishing between different theories. 

It is possible to account for the observed polar intensity with extra- 
galactic sources using Relativistic Cosmology, but not on the Steady- 
State theory. Recent measurements of the Andromeda nebula have, 
however, suggested that part of the polar emission may be due to galactic¢ 
sources and a conclusive distinction cannot therefore yet be made. 


§ 1. INTRODUCTION 


WESTERHOUT AND OorT (1951) have shown that the survey of galactic 
radiation at a frequency of 100 Mc/sec made by Bolton and Westfold 
(1950) may be explained by assuming that most of the radiation comes 
from radio stars. In their model these are distributed through the galaxy 
in the same way as the common Population II stars of types G and K. 
They pointed out, however, that in order to account for the observed 
distribution, it was necessary to superimpose on this model an isotropic 


component of the radiation. At 100 Mc/see this additional intensity 
ee een ea 
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corresponds to a brightness temperature 7', of 600°K, where 7’, is defined 
by the Rayleigh—Jeans law :— 
Biv) c? 


Qkv? ? 


T= 


where v is the frequency, B(v) is the brightness in units of flux density per 
steradian, / is Boltzmann’s constant and c is the velocity of light. 

In §2 a number of possible causes of this isotropic component are 
discussed, one of which is related to the emission from extra-galactic 
sources ; this possibility is examined in detail in §§ 3 and 4. 

In addition to the emission from normal extra-galactic nebulae, recent. 
observations have revealed the existence of a different type of radio source 
whose emission exceeds by a factor of some 10° that of the normal nebulae. 
Following the accurate determination of the position of the intense radio 
star in Cygnus by Smith (1951), Baade and Minkowski (1954 a) were able 
to identify it with a very peculiar galaxy whose spectrum shows the 
presence of high excitation in the interstellar material ; they concluded 
that the observations could be explained if the object were in fact two 
galaxies in collision. Subsequently two similar events have been located 
at the positions of other radio sources. 

From data on the spatial density and size of galaxies an estimate can be 
made of the spatial density of galaxies in collision at any time. If the 
relationship between the type of collision and the radio emission were also 
known, it would then be possible to derive the contribution of such objects 
to the radio emission per unit volume of space. 

In order to deduce the total intensity produced, it is necessary to consider 
the effect of the red-shift. It is shown in § 4 that owing to the difference 
between the spectral law of the radio sources and the spectral law of 
visible stars, the integrated radio intensity will be affected to a much 
greater degree by the more distant parts of the Universe and the total 
thus depends critically on the Cosmological theory adopted. 


§ 2. PosstBLE SOURCES OF THE ISOTROPIC COMPONENT 


Westerhout and Oort suggested three possible explanations for the 
radiation not accounted for by their model :— 

(a) Radiation from the ionosphere. 

(b) Galactic sources with a very unusual spatial distribution. 

(c) The integrated radiation from all the normal extra-galactic nebulae. 

In addition the following will be considered :— 

(d) Radiation from ionized. hydrogen clouds surrounding the sun. 

(e) The integrated radiation from colliding galaxies. 


These will now be dealt with in turn :— 

(a) Westerhout and Oort point out that if this mechanism were of 
importance considerable diurnal variations in the observed intensity would 
be expected. Apart from this the optical depth of the ionosphere 
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(r=fx(v)ds where x(v) is the absorption coefficient and ds is an element of 
length) may be calculated. It may be shown that for a reasonable 
distribution of electron density and temperature with height the radiation 
at 100 Mc/sec corresponds to an additional temperature of only about 
5°K. 

(b) If the large polar intensity were to be explained in terms of galactic 
sources, it would be necessary to suppose that they were distributed to 
large distances from the galactic plane, and also that they extended to 
radial distances larger than that of the solar system. No such population 
of optical stars appears to have been observed. 

(c) The early observations of the nearer extra-galactic nebulae suggested 
that there would be little hope of accounting for the isotropic component 
in terms of the emission from normal galaxies. Hanbury Brown and 
Hazard (1952) however, later suggested, on the basis of observations of the 
Perseus cluster, that there was a discrepancy between the radio and optical 
emission from faint galaxies; on this supposition the integrated radio 
emission would then be considerably greater than that previously 
estimated. Recently Baldwin and Elsmore (1954) have made some de- 
tailed measurements of the Perseus cluster and have shown that 75% of 
the radiation is due to NGC 1275, one of the members of the cluster. The 
remaining emission is in agreement with the value expected if the other 
members of the cluster have a radio emission comparable with the nearby 
galaxies which have been observed individually. It is, therefore, con- 
cluded that the ratio of optical to radio magnitudes is constant, and in 
later sections a value of the radio emission from a galaxy of optical 
magnitude +10 will be taken as 5-3 x 10-2 w m-? (c/s)! at 100 Mc/see 
(cf. Hanbury Brown and Hazard 1953a). As will be seen in § 4 this 
figure is inadequate to account for the brightness temperature of 600°K 
observed in polar directions. 

(d) Siedentopf (1953) has shown from measurements of the polarization 
of the zodiacal light that the electron density in the solar system is much 
higher than was previously thought to be the case. Strémgren in 1939 
gave an expression which relates the diameter of the region within which 
the interstellar gas is ionized, to the temperature and radius of the 
exciting star and the density of the gas. Using the values of electron 
density determined by Siedentopf and assuming an average electron 
temperature of 10 000°K, it is possible to compute the optical depth of the 
entire region as viewed from the earth. The results obtained show that 
the contribution to the brightness temperature at 100 Me/sec is unlikely 
to exceed 1°K,. 

(¢) The contribution due to colliding galaxies is considered in more 
detail in the following paragraphs. 


§ 3. THE OBSERVATIONAL Data 


In this section an attempt will be made to derive from the very limited 
observational results enough numerical data to enable an estimate to be 
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made of the emission per unit volume due to colliding galaxies. It is 
first necessary to find their spatial density and the mean radio emission 
from such sources. 


3.1. The Spatial Density of Colliding Galaaies 


Suppose that space were filled with a uniform distribution of galaxies, 
each having a radius a parsecs and a random velocity V parsecs per second, 
the average spatial density being p galaxies per cubic parsec. If it is 
assumed that an average collision occurs when the centres of two galaxies 
pass within a distance a of each other, the collisional cross section may be 
written 7a? and each galaxy will make 1/27a?Vp collisions per second. 
The number of collisions per cubic parsec is then given by 1/27a?V p? and if 
each collision is assumed to last for a time a/V, the number of collisions 
occurring at any time in each cubic parsec is n= 1/27a'p?. 

The greatest uncertainty in deriving n is due to the difficulty of assessing 
the mean density of galaxies. As has been shown by many workers, 
galaxies exist in clusters, the density in these being anything up to one 
thousand times that of the field nebulae. Since the collisional probability 
is proportional to p?, it might be thought that clusters only need be con- 
sidered as contributing to collisions. The fraction of space occupied by 
dense clusters, however, is probably small and also the proportion of 
elliptical galaxies in clusters is much greater than among the field nebulae. 
Indeed Spitzer and Baade (1951) have suggested that this characteristic 
is due to the removal of the interstellar gas by collisions. The develop- 
ment of spiral nebulae is thus prevented and the radio emission, which on 
present evidence should be related to the gaseous material, may be smaller 
for collisions involving elliptical nebulae. 

Since there are no reliable data on the fractions of space having different 
galactic densities, the present calculations are based on the mean density, 
4x 10-18 galaxies per cubic parsec (Neyman, Scott and Shane 1953), 
taking in this and further calculations, the new distance scale to be twice 
the old. Using this figure and taking a as 10 000 parsecs, the collisional 
density n=7 x 10-8 per cubic parsec. Since the effect of clustering will be 
to increase the number of collisions, this represents a minimum value. 
There will also be a contribution from peculiar galaxies such as M 87 in 
Virgo, which has an abnormally intense radio emission, but this cannot at 
present be estimated. 


3.2. The Known Colliding Galaxies 


3.2.1. Cygnus—A, 19-01 (Ryle, Smith and Elsmore 1950), 19-+4 (Mills 1952) 

From the photographs obtained by Baade and Minkowski (1954 a), 
two late-type spiral galaxies appear to be in broadside collision, since 
emission lines are prominent over the whole disc. Using the observed 
red-shift, which corresponds to a recessional velocity of 16 830 km/sec the 
relationship between red-shift and distance (corrected for the new 
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distance scale) gives the distance as 66x 10° parsecs, and the measured 
intensity at the earth is :— 

1350 x 10-25 w m~ (c/s)! at 81-5 Me/sec 
and 1300 x 10-25 w m-2 (¢/s)~1 at 100 Mc/sec. 


3.2.2, NGC 1275, 03-02 (Ryle, Smith and Elsmore), 03+-4 (Jills) 


This object is interpreted (Baade and Minkowski 1954 b) as a tightly- 
wound spiral penetrating a late-type spiral galaxy. From the red-shift 
the distance is found to be 20 x 10° parsecs. 

The measured intensity at the earth is :— 

12x 10-25 w m=? (c/s)! at 81-5 Me/sec (Baldwin and Elsmore 1954). 


3.2.3. NGC 5128, Centaurus—A, 13—4 (Mills) 

Baade and Minkowski believe this to be a late-type spiral in collision 
with an elliptical galaxy. De Vaucouleurs (1953), however, is of the 
opinion that NGC 5128 and several other similar objects he has found in 
the southern sky near radio sources, resemble the nuclei of peculiar barred 
spiral galaxies. 

The distance is greater than 10° parsecs, and the measured intensity 
is 16 x 10-25 w m- (c/s)~1 at 101 Me/sec. 


3.3. The Radio Emission Associated with Colliding Galaxies 


Only three examples of colliding galaxies are known at present and 
these differ in the type of galaxy involved and in the geometry of the 
collision. They also show marked differences in the intensity of the radio 
emission, as can be seen from the following table which gives the power 
radiated per unit bandwidth per unit solid angle ; these figures have been 
corrected to a frequency of 100 Mc/sec assuming that the spectral law is 
dt ceirecie: 

Cygnus—A : P=5 x 10?° w steradian—! (¢/s)-1. 
NGC 1275: P=4x 103 w steradian—! (¢/s)-}. 
NGC 5128: P> 10? w steradian—! (c/s)—}. 


There is as yet no satisfactory mechanism to account for the emission of 
radio waves from colliding galaxies, but it seems probable that the 
radiation is due to the interaction of moving gas clouds, since all the radio 
sources which have been identified with visible objects show the presence 
of gaseous material in violent motion. Spiral galaxies such as those in the 
Cygnus source contain widely distributed interstellar material, whilst 
elliptical galaxies contain practically none. It has been suggested by 
Baade and Minkowski that the radio emission from the Centaurus source 
may be due to the tidal forces causing differential motion of the gas. This 
might be expected to produce relatively less intensity. 

Without a detailed theory for the emission mechanism it is clearly 
impossible to say how the emission might vary with the density and 
relative velocity of the gas clouds, or how long the radiation might 
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last after the galaxies have penetrated one another. These questions 
may be answered if a larger number of identifications can be made, but 
it is necessary to see what estimate of the total emission can be made from 
the present limited data. 

The observations suggest that a collision between two late-type spirals 
is very much more effective than other collisions, but the spectroscopic 
evidence of excitation over the whole disc indicates that the Cygnus 
source itself may be a rather special example having a rare geometrical 
configuration ; the emission may therefore exceed that expected from a 
typical collision. This conclusion may easily be verified from radio data, 
for if the Cygnus-type sources had a density of 7 10~*3 per cubic parsec, 
previously estimated as a minimum value for colliding galaxies, then over 
the whole sky there would be some 100 000 radio stars with an intensity 
greater than 1074 wm (c/s)-!. Using preliminary data from the new 
Cambridge survey of radio stars, it may be estimated that there are some 
60 stars with a flux larger than this value. This then fixes a maximum 
permissible density for Cygnus-type sources ; and their radio emission per 
unit volume is less than 1/10 of that of normal galaxies. 

Having dismissed Cygnus—A as a special case, the consequences of 
assuming a mean power for collisions of the order of that observed from 
the Perseus and Centaurus sources must be examined. The loss of 
energy by radio and optical radiation in the time required for the stars in 
two galaxies to pass by each other, a/V, is very small compared with the 
collisional energy of the gas clouds, and since the total emission per unit 
volume depends on the time during which the radiation persists, it may be 
that the effective source density is considerably greater than that already 
estimated. An upper limit may again be placed on this increase by the 
number of sources actually detected. For example, if 1/3 of the 60 
sources with flux larger than 1024 w m- (c/s)~1 were colliding galaxies then 
the density could not be increased by more than a factor of 10 to 
n=7 10-22 per cubic parsec. The emission per unit volume due to 
colliding galaxies would then become nP=7 x 10-22 x 3x 107—210 w 
steradian—! (c/s)-! parsec-? which may be compared with the figure due to 
normal galaxies of 216 w steradian~t (c/s)~* parsec”* (n=4 X 10718 per cubic 
parsec and P=5-4x 10" w steradian* (c/s)-?). 

Thus on this basis we should expect the radio radiation from a large 
assembly of galaxies to be about twice that calculated for the normal 
galaxies alone. 

Independent evidence for such a conclusion has been provided by the 
observations of Kraus (1954), and Hanbury Brown and Hazard (1953 b) 
who have related certain irregularities in the distribution of radio emission 
to large scale irregularities in the distribution of extra-galactic nebulae. 
The latter authors conclude that the emission per unit volume exceeds by 
1-7 magnitudes (a factor of 4-8) that expected from a consideration of the 
emission of normal galaxies. More accurate observations of such large- 
scale structure have recently been made by Baldwin (in preparation) with 
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the new Cambridge radio telescope, and he finds that the factor is probably 
at least 5. (The optical data on galactic counts are insufficient to get a 
very precise figure.) It therefore seems likely that the emission per cubic 
parsec lies in the range 400-1200 w steradian-! (c/s)~1. In the next 
section the integrated intensity associated with such an emission will be 
derived. 


§4. Tue InreGRATED RADIATION FROM ExTRa-GALACTIC SOURCES 


Before the integrated intensity due to distributed sources of radiation 
can be derived, it is important to examine the effects caused by the 
apparent expansion of the Universe. Because of the red-shift, the 
spectrum and intensity from radiating objects will be dependent on their 
distances. At optical wavelengths the spectrum from nearby galaxies 
approximates to that from a black body at a temperature of about 
10 000°K, so that galaxies situated at distances where the apparent: 
recession velocity is of the order of half the velocity of light will produce 
relatively little radiation in the visible spectrum, owing to the rapid fall of 
intensity on the short wavelength side of the black body curve. It 
therefore becomes extremely difficult to observe optically at distances. 
greater than this, and the integrated optical radiation is highly convergent. 

The radio spectrum both of our own galaxy and of the intense extra- 
galactic sources, however, shows a comparatively gradual decrease towards 
the shorter wavelengths. It should, therefore, be feasible to detect indivi- 
dual sources at very much greater distances, whilst the integrated radiation 
will depend more on the distant regions of the Universe. It is of interest 
to note that although the Cygnus source is at a distance about 1/10 of 
that to which optical telescopes have penetrated, its intensity is of the 
order of 10% times that of the weaker radio sources now being 
detected ; it therefore seems probable that a number of the latter are 
beyond the range of optical instruments. 

It is now possible to derive the integrated radiation at a given wave- 
length in terms of a given Cosmological theory, if the spectrum of the 
source is known. 

Let the total power per unit volume radiated in the frequency range v to 
oe Av be 47nP(v) dv w parsec~* then the total flux observed at the earth 
is 
°nP(v’) Av’ = rdr 
0 M(1+2)? * (1-2) 


where z is the red-shift, r is the distance, v’ is the frequency at which the 
radiation was emitted, that has been shifted to the frequency of reception v, 
Av’ is the bandwidth of emission, Av is the bandwidth of reception, and « 
is a constant factor=0 for Relativistic Cosmolo é 
=3 for the Steady-State ho A Bondi Pera 
The circumstances under which this expression holds are considered later. 
The red-shift z is observed to be proportional to distance and z=r/T 


Biv) Av= | 
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where 1/7’ is Hubble’s constant (corrected for the new scale), also 
v'/v= Av'/dv=(1+2)=(1+7/T). Substituting in the expression above, 
© nP(v') dr 
Bo | pay 
Sciama (personal communication) has shown that it is unnecessary to take 
into account the optical energy range when this has been shifted to radio 
frequencies. 
Assuming that P(v) ocv-%5, then P(v’)/P(v)=(v'/v)5=(1+7/T)°? and 
nP(v) dr nP(v)T 
Bw)= [os = Ra 
o (taer/T)eet= (05+) 
Boje? __ nP(v)T A? 
2ky? ~~ 2k(0-5-+a) * 
It should be noted that this result is correct for the Steady-State theory, 
but holds for Relativistic Cosmology only if it is possible to neglect 

(a) Curvature, 

(b) the variation of P(v) with time (since the Universe is evolving and 
the galaxies observed at greater distances will be at an earlier stage ot 
development), 

(c) a possible variation of the red-shift vs. distance relationship. 

The importance of these corrections will depend on the convergence of 
the expression for the integrated radiation. To test this convergence it 
can be integrated to a distance R at which the apparent velocity is half 
that of light. This distance can be taken as approximately the ‘ limit ’ of 
the optically observable Universe and the distance beyond which the 
effects mentioned will play an increasing part. In this case 

vie=r/T «. R=1/2T 
nP(v)T 
and [B(v)ye=4. hoeey : 
Hence only 1/3 of the total radiation comes from a distance less than 
1/27’, and it is possible that the actual radiation expected using Relativistic 
Cosmology may be rather different from that calculated here. 
However, taking A=3 m=10~1° parsecs, 
k=1-4 x 10-*3 joules °K~, 
T' =1-1 x 10° parsecs 
and substituting nP=400—1200 w steradian-} (c/s)~! parsec-* in the 
expression for 7',, then 7’, lies in the range 
314°-942°K when «=0, i.e. for Relativistic Cosmology, 
45°-135°K when «=3, ie. for Steady-State theory. 


Therefore T ».= 


§ 5. CONCLUSIONS 
The observed intensity plot of radio emission shows a greater value 
towards the galactic poles than can be accounted for by sources distributed 
like the common optical stars. The integrated emission from extra- 
galactic sources, and the existence of a population of galactic sources 
having a wide dispersion have been suggested as possible explanations. 
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The determination of the total emission from extra-galactic sources 
depends on a knowledge of the average emission per unit volume of extra- 
galactic space, and on an investigation of the effects of the red-shift. 

Estimates are made of the emission per unit volume due to normal 
galaxies, as well as that to be expected from the rare but intense colliding 
galaxies. Other evidence has been obtained of the total emission per unit 
volume by observations of the enhanced radiation from large scale 
irregularities in the galactic distribution. 

It is found that the integrated intensity depends strongly on the Cosmo- 
logical theory assumed. Using the data derived here it is seen to be 
possible to account for the polar intensity with the extra-galactic sources 
alone on the Relativistic theory but not the Steady-State theory. The 
present results cannot yet, however, be used to distinguish conclusively 
between the two theories, both because of the possible errors in the 
assumptions made in the calculations, and also because recent observations 
of the Andromeda nebula by Baldwin (1954) have suggested that the 
radio emission extends over an appreciably larger angle than the optical 
emission. There is, therefore, some evidence in favour of a population of 
radio sources in the Andromeda nebula which has a more extended distri- 
bution than the optical stars; if such a population exists in our own 
galaxy also, it might be sufficient to account for a considerable fraction of 
the observed polar intensity. More accurate observations may allow 
the contribution of extra-galactic sources to be derived well enough to 
enable a definite distinction to be made between the two theories. 


This work was carried out at the Cavendish Laboratory as part of a 
programme of radio research supported by the Department of Scientific 
and Industrial Research, to whom the author is indebted for a main- 
tenance grant. He also wishes to thank Professor H. Bondi for a useful 
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CXXVII. Double Pion Production and the Nature of Pseudoscalar 
Coupling 


By P. Bocconterr* and G. FritpMan 
Department of Mathematical Physics, The University, Birmingham? 


[Received May 26, 1954] 


ABSTRACT 


A perturbation calculation for single and double pion production is 
carried out for pseudoscalar, symmetric meson theory at energies in the 
Birmingham range (about 1000 Mev). It is found that the y, nature of 
the coupling does not favour double over single production and that, in 
fact, a large coupling constant would be required to explain a large double 
to single ratio. 


§ 1. INTRODUCTION 


PRELIMINARY experimental results (Fowler, Shutt, Thorndike and 
Whittemore 1953) indicate that double pion production in nucleon— 
nucleon collisions is rather large at nucleon energies of about 2-3 Bev in 
the laboratory system ; and that at these energies double pion production 
is probably larger than single production. This result would, at first 
sight, indicate a large coupling constant between nucleons and mesons. 
However, one might imagine that if the coupling between nucleons and 
mesons is pseudoscalar (y;), this would favour the occurrence of virtual 
pairs (in terms of perturbation theory). In that case, those reactions 
containing an even number of vertices (e.g. double pion production) 
would be favoured compared to those containing an odd number of 
vertices (e.g. single pion production). This presumed effect of y; should 
make itself more apparent if the coupling constant were indeed not too 
large. 

It is the purpose of this paper to carry out a perturbation calculation for 
double and single pion production in order to see whether a large double to 
single ratio might be explained in terms of y,; coupling rather than in 
terms of a large coupling constant. We find, in fact, that the y; nature of 
the coupling alone does not explain a large double to single ratio. In spite 
of the unsatisfactory nature of a perturbation calculation, we feel that it is 
worth putting such a calculation on record since, like similar calculations 
for other reactions, it may be useful for orientation. 

In § 2 we evaluate the densities of final states arising in single and double 


production and make a crude estimate of the matrix elements using 
ee ee 
* On leave from the Istituto Nazionale di Fisica Nucleare, Sez-di-Milano. 
+ Communicated by Professor R. E. Peierls, F.R.S. 


SER. 7, VOL. 45, NO. 370.—NOV. 1954 4F 


1146 P. Bocchieri and G. Feldman on Double Pion 


arguments similar to those of Fermi (1951). In § 3 we show how to write 
down, simply, the terms occurring in the cross section and we evaluate 


them. The results are discussed in § 4. 


§ 2 
8 3 
Throughout the paper we shall be considering the following two 
reactions : 


NitN,g>NgtNatm, . . 2 ss ss ee 

NiAN,>NatNatmtm, ...... . (2.2) 
where NV, are nucleons and 7, are mesons. In the laboratory system the 
threshold for (2.1) is 290 Mev and for (2.2) it is 600 Mev. We shall be 
calculating cross sections for an incoming nucleon of kinetic energy 
940 Mev (i.e. equal to M the nucleon rest mass). In the C.M. system of the | 
two initial nucleons the kinetic energy available is 0-45 M. The mass of 
the 7-meson is 0-15 M, and thus in reaction (2.1) the kinetic energy to be 
shared between the final particles is 0-30 W/, whereas in (2.2) it is 0-15 M. 
In evaluating the integrals from the density of final states below, we shall 
assume these energies to be non-relativistic. In view of the errors of other 
approximations to be made, this would seem adequate. 

The cross section for either (2.1) or (2.2) is 


Qa 
o— = || HP p(B) dQ, +e ate eke 


where H is the matrix element, p(Z) is the density of final states and dQ 
the suitable differential factors. v is the relative velocity of the incoming 
nucleons. We use units in which h=c=1. 

For single production 


ap, d3 
px(B) 4Q,=(ps-+ pata) (PAPE 4 Fg 1) dp, Bet. (2 


where ps, p, and q are the momenta of re outgoing nucleons and meson 
respectively, » is the meson mass and M the nucleon mass and 7’, is the 
kinetic energy available in the C.M. 7,=0-30 M for an incoming nucleon 
of 940 Mev kinetic energy. 

If we assume the matrix element is energy independent (or varies 
slowly with energy) we can integrate over the variables dQ, in (2.3) easily 
by making the transformation 


: Ps+Pa=a, P3—py=b. 
This gives 


M8 Ree 3/2 /T.\2 
Joxz) )dQ,= qs ¥ QM uti T 2x3 —— ~ Tap (4) (=) M3 p2. (2.5) 
For double Ree 
2 2 
{ p2(H) dQo= | (Ps +Pat+4i+q2)d (Pop + oH bes —1,) 


fe 
dq, a 
x dp, dp, re (2.6) 
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We can carry out this integral by making the substitution 


Ps+Pa=a, Ps—Pa=b | 
(2.7) 
and Gitqde=%, q3—q.=B. 
This gives 
1 1 128 /M\3/2 Te 72 
ee a Se eg ae -3 [2 8 
iy 1 ZS fea N Sia TION Hy ; 


It is of interest at this stage to make an estimate of the matrix element 
as carried out by Fermi (1951). 

Each meson vertex in fig. 1 contributes f/1/2w where f is the coupling 
constant (essentially scalar coupling) and w is the energy of the meson. 
This gives for the matrix element near threshold 


AJ 2. fee. aac eee gw PD) 
Fig. 1 


Single meson production. 

(The result is slightly different from Fermi since we have allowed for 
various permutations of the vertices.) For an incoming nucleon of 
940 mev we get for the matrix element 

fee et eae. es 3 (2.10) 
We shall crudely fix the coupling constant by fitting the results to near 
threshold single pion production, that is, at about 340 Mev. A 340 Mev~ 
incoming nucleon corresponds to 7';=§ and we get 


f2\3 
o,, (340) =1-5 x 125 (5) TA ar eer e es . ( abl) 
ct 4qr 


Assuming the cross section is 0-6 mb this gives f?/47~1/7. Foranenergy of 


940 mev 7',=2p and we get 
ALOE” SF ie Ss eee CR), 


We compare this with the result for total p-—p cross section for 1 Bev 
protons. Brookhaven (Leavitt and Chen 1954) find a total p—p cross 
section of 49 mb. The elastic cross section at 400 Mev is about 25 mb. 


4¥F2 
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We now estimate the double meson production cross section at 940 Mev 


and we get 


AoA y oie . + 4 oe ee 
At 940 mev, 7',=p and we are led to the following ratio : 

ml 980) Se 4 Oe 

o,(940) 8 


It is of interest to extrapolate these results to Brookhaven energies 
(2:3 Bev incoming nucleon). The validity of this extrapolation depends 
on (a) the ratio of the matrix elements for double and single production 
remaining fairly constant with energy and (6) the non-relativistic approxi- 
mation to the density of states (2.5) and (2.8) being valid at this energy. 
Although we could evaluate the density of states more accurately we feel 
it is not worth the trouble in view of (a). Under these assumptions we 
obtain 


pe RR 


~~ 


02,(2300) 3-5 
g,(2300) 1’ 


$3 
In this section we evaluate by means of a graphical technique the squares 
of the matrix elements for single and double meson production. These 
have already been calculated for single production by Morette (1949). 


First let us consider single production. One can draw the Feynman 
graphs for this reaction. One of these is shown in fig. 2. 


Fig. 2 


Typical Feynman graph in single pion production. 


There are eight such graphs, four arising from the meson ‘ q’ coming off 
each leg and the same number when particles 1 and 2 are interchanged 

Y : . - . . . . . : 
Thus the matrix element H, will contain eight terms. Since we are going 
to sum over final and average over initial spins and isotopic spins it is 
easier to work with the square of the matrix elements | H 1 |? which con- 


“eee” 
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tains 64 terms. We can directly draw graphs* for the terms appearing in 
| H,|?. They fall into two classes which we call the ordinary terms and 
the exchange terms. An example of an ordinary term is given in fig. 3(a). 
This is just a graphical representation of the square of the term given in 
fig. 2. The prescription used in writing down the terms which arise is as 
follows. 


yptiM 
Each —>——. gives oD eae WER TET 
: M yp+iM 
Kach —+—0-<— gives A,(p)=— — > . 
Z a 2M (3.1) 
Each — -— =- 7 A So 
aC See gives F(P) pepe 


; 1 
Each —— —-0—<-- gives Spe 845: 


Each vertex gives gy,7;. (pis a four momentum and yp=y, p,.) One 
must multiply by 1/16 coming from the averaging over the initial spin and 
isotopic spin states. All of the 27 factors have been taken care of in the 
definitions (2.3) and (2.4). We have assumed the interaction Hamiltonian 
to be ‘ 
H (x)= g(x) y 53) (&) 7h ;(2). Mee | (3-2) 
There are 32 of these ordinary terms. They arise when the meson line 
q begins on some line in the lower half of the diagram and ends on some 
line in the upper half. Each closed loop gives rise to a trace over the Dirac 
matrices and isotropic spin matrices. One can easily see that 16 of these 
terms are zero. They are those diagrams in which the meson q begins in 
one loop and ends in the other. This gives rise in y, coupling to three y,’s 
in each loop and the trace of three y;’s with the other y’s is zero. An 
example of the exchange terms is given in fig. 3 (b). This is one continuous 
loop and so gives rise to one trace which will be more difficult to evaluate. 
There are 32 such terms. 
The expression for | H, |? corresponding to fig. 3 (a) is 
es 1 1 (po - Pat M")[2(p1 - W)(Ps - q)+u7(p1 - Ps t+M?)] 
(N330 BEE B st (Pi -g-+H FP: Pat IPP 
a WE pe (3.3) 
At the energy of double pion production threshold #,=H,=M-+p, and on 
the average Hz=H,=M-+ p/4 and qy=3/2u. Neglecting powers of u/M, 
this term contributes 
eh Wee) a 943.4) 
40 p?M>* 
And finally the total matrix element is (the exchange terms contributing 
less than 10% at this energy) 


aus (3.5) 
B pews wel eee 


* This suggestion was made by Professor R. E. Peierls. 
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We now come to the graphs in double production. A typical graph is 
shown in fig. 4 (a). There are altogether 40 such graphs in the matrix 
element H, and thus the square of the matrix element | H, |? will contain 
1600 terms, 800 of which will be ordinary terms of the form shown in 
fig. 3 (a) (containing an extra meson line, however) and 800 will be of the 
exchange type (fig. 3 (b)). Now, as in single production a certain number 
of the ordinary terms will be zero whenever there are 3 y,’s in one loop. 
There are 384 such terms. The others can be evaluated using the 


prescription of (3.1). 


(a) 


vnc IQs . PAG 1 
Graphical representation of the square of the matrix element in pion production. 
(a) An ordinary term. (b) An exchange term. 


In evaluating the various terms we may notice just what effect the 
coupling has. Naively, one might expect that the contribution front 
fig. 4 (a) would be greater than that from fig. 4 (b) ; since if we believe that 
in a time ordered event y; has large matrix elements only between positive 
and negative energies, then the largest part of 4 (a) is the time ordered 
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graph 4 (c) consisting entirely of ‘ large’ matrix elements, whereas the 
largest part of 4(b) is 4(d) which contains only two ‘large’ matrix 
elements. However, it turns out that both 4 (a) and 4 (b) contribute 
about the same amount to the total matrix element. 


Fig. 4 


Typical Feynman graphs in double pion production. 


One can easily see* why the graphs 4 (c) and 4 (d) are of the same order 
of magnitude. In the three-dimensional treatment, the matrix element 
between a positive and negative energy component is of order 1, whereas 
between two positive energy components it is of order | p |/ where | p| 
is the momentum of the nucleon. At threshold, each incoming nucleon 
must have a kinetic energy (in the C.M. system) equal to » and therefore 
its momentum is \/(2uM). Thus the vertices in 4 (c) contribute (1)4=1 
and those in 4(d) contribute (1)?[1/(2u/M)P=2u/M. The energy 
Pebieess tS es eS 


* This very simple argument is due to Dr. G. E. Brown. 
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denominator of 4(c) is (1/2M)#1/u and for 4(d) it is 1/2M(1/p)2. Thus 
both 4(c) and 4(d) are of the same order of magnitude. 


Fig. 4 


The corresponding time-ordered graphs. 


The square of the term in fig. 4 (a) gives (at threshold) 
ee) 
) | 36 Me > a WR vehi eh Lt (3.6) 
neglecting higher powers of u/M; and the square of term 4(b) gives 
one-quarter of this result. 


After a good deal of algebra we arrive at the final answer for | H, |? 


(which includes the sum and average o final initial i . : 
respectively) ge over final and initial isotopic spins 


ete) 
167 uiMs* oi 0 fod 2 ee ns (3.7) 


Lye 
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The results for the cross sections are 


BSN fe \e? 4 ; TNA T \2 of Ne TE 3 

on (Z) (47) 3 Ven (=) (=) eas ‘e i yeste, 
g?\4 (w\528 (T\72 (1\2 T.\7/2 : 

a (=) @ 3(-2) (a) ~0-006 (=) (2) mb, (3.9) 


where 7',=2yu and 7'.,= p for a 940 Mev incident proton. 


This gives 
Fon(940) 1 g? 
o,(940) ~ 50 (5) 
If we fit the coupling constant to o,(940)~%25 mb, we find g?/4n~4-5 
giving a result of the order of 1 to 10, for the ratio of double to single 
production at this energy. If we extrapolate the results to Brookhaven 
energies (where 7', ~ 6u, T', = 54) under the same crude assumptions 
as stated in § 1 we find 


(3.10) 


O29(2300) 7 g? 26 

o(2300) 1247-1 * 

We may compare these results with the still cruder ones (2.14) and 
(2.15). 


(3.11) 


§ 4. CONCLUSIONS 


The results of the above calculation indicate that one must be careful 
when using arguments about y,; coupling. One cannot say that a reaction 
with an even number of vertices will be favoured over one with an odd 
number. Nor can one say that a graph with the maximum number of 
pairs produced will be larger than another graph. It is important to 
examine the momenta involved. A reaction at threshold need not 
necessarily involve small momenta. Double pion production requires 
each incoming nucleon to have a kinetic energy of at least ». This is 
essentially why the matrix elements shown in figs. 4(a) and 4 (6) 
are of the same order of magnitude. f 

If the ratio of the cross sections for double to single pion production 
is indeed large at Birmingham and Brookhaven energies this would 
indicate a large coupling constant for pseudoscalar theory but would 
give little indication as to whether the coupling is pseudoscalar or not 
(cf. eqns. (3.10), (3.11) and (2.14), (2.15)). 


We should like to thank Professor R. E. Peierls and Dr. P. T. Matthews 
for many interesting and valuable discussions. 
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SUMMARY 


The paper (which makes no claim to originality) gives an elementary 
survey of the physical causes believed to underly the more energetic 
features of the general circulation. 

The geographical distribution of temperature implies hydrostatically 
a strong fall of the pressure from the tropics to the poles at any fixed 
altitude in the upper troposphere. This produces geostrophically 
(in the main) the strong westerly winds found at all latitudes in that 
region. Simple estimates indicate why a typical wind speed is 20 m/sec. 

The wind distribution near the ground, a secondary effect, is intelligible 
only if there is everywhere a poleward flow of west-east angular momentum 
(taken about the earth’s axis). The required flow reaches a maximum 
of around 3 x 1026 dyn cm, at latitude 30°. In the tropics and sub-tropics 
the angular momentum must flow against its own gradient. It can do 
so to the necessary extent, however, as a result of the poleward motions 
down pressure gradients in the upper troposphere which are needed to 
maintain the energy of the general circulation, which is dissipated at 
a rate of around 5 x 10% ergs/em?/sec. In the cyclone belt, the assump- 
tion behind this statement (that isobars in the upper troposphere lie 
east-west, on the whole) is constantly violated, but there the flow of 
angular momentum (down its own gradient) presents no difficulty, 
and is doubtless achieved by the cyclones themselves. 


$1. InrTROoDUCTION 


THIS paper is intended as a contribution to the relations between 
meteorologists and general physicists. It concerns the ‘explanation of 
the general circulation ’. 

A physicist, wishing to arrive at some understanding of the basic 
physical causes underlying the general circulation of the atmosphere, 
cannot easily do so from a survey of the best known accounts of 
meteorology. He discovers early that the schoolbook descriptions of 
hot air rising at the equator and flowing towards the poles and being 
deflected into a westerly direction by Coriolis force are quite inadequate 


ee ee ene aa ae ne a 


* Communicated by the Author. 
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to explain the strength of the Westerlies in the temperature zones. 
Apart from this, he is told, for example, by Brunt (1939), in his invaluable 
Physical and Dynamical Meteorology, that ‘it is not possible at present 
to put forward any satisfactory theory of the general circulation ” 
(p. 413), “ nor at this stage can we say that the cyclone, anticyclone and 
general circulation are not merely three aspects of one problem ” (p. 403), 
and that ‘“ the latest writings on the frontal analysis of depressions, on 
the anticyclones, and on the general circulation, all indicate a recognition 
of the complexity of the phenomena, and less readiness to attempt 
general theories ” (p. 418). 

All this is very damping, as the general circulation is the most energetic 
and permanent feature of the atmosphere, and a physicist is inclined 
to believe that the general laws of mechanics and thermodynamics 
should be sufficient to give an order-of-magnitude explanation of such 
a predominant feature of a relatively accessible part of the universe. 
The details of the instability of fronts, and of the production of periodic 
quasi-regular disturbances of the general circulation, are problems of 
prime importance to the meteorologist, but secondary from the point of 
view of the physicist, who would not expect a simple description of these 
elaborate perturbation effects in terms of basic laws. 

In fact, however, putting aside the finer details (with which the mete- 
orological experts are naturally deeply concerned), enough information 
exists already to quieten the physicist’s fears lest the rules of physical 
reasoning are inadequate to illuminate even the broadest general features 
of the situation. The following remarks are not intended to be in any 
way original, and their content will mostly be familiar to meteorologists. 
They are put forward, principally, as a sample of what can most appro- 
priately be said about the subject to physicists in general—for example, 
when the topic is mentioned in sixth-form or university physics teaching. 

The primary phenomenon is the tendency to strong westerly mean 
circulation, at all latitudes, in the upper part of the troposphere. ‘This is 
treated first, in §2. The distribution of mean winds near the ground is 
a kind of by-product, whose discussion is postponed to §§3 and 4. 


§2. THE CIRCULATION IN THE Uprer PARTS OF THE TROPOSPHERE 


In terms of energy content, the mean wind distribution at the earth’s 
surface is a secondary effect compared with the much stronger mean 
circulation at greater heights. In the upper half of the troposphere the 
prevailing wind is westerly at all latitudes, with mean velocities of the 
order of 20 m/sec. It is natural to try to explain this, the most energetic 
feature of the situation, first. 

The circulation is associated with a mean pressure distribution in the 
upper parts of the troposphere which falls fairly steadily (at any fixed 
height above sea level) from the tropics to the poles. The mean isobars 
have a clear tendency to lie in east—west directions, and they encircle 
the earth without a break except for certain * highs ’ over the continents 
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in summer. The mean pressure excess at the equator over that at the 
same altitude at the poles is of the order of 50 mb and is considerably 
more than the difference at sea level. 

One needs to explain only this mean pressure distribution, for this 
appears quite sufficient to maintain the observed mean circulation as, 
in the main, a ‘ geostrophic ’ wind. In other words, the wind is approxi- 
mately along the isobars and provides, on the average, as much Coriolis 
force per unit volume as will balance the poleward force due to pressure 
gradient. In symbols, it satisfies 


We - Op 
CORN IES eee eee 
where w is the earth’s angular velocity (7-3 x 10° s~+), ¢ is the latitude, p 
the density, w the velocity along the isobars, and dp/dy the pressure 
gradient normal to them. 

Now the pressure distribution seems clearly explicable as a direct 
result of the temperature distribution, which in turn is maintained by 
the distribution with latitude of incoming solar heat. The mean tempera- 
ture decreases, of course, from the tropics to the poles, and the pressure 
gradient in the upper part of the troposphere is an inevitable hydrostatic 
consequence. For the pressure must fall off with altitude more gradually 
in the tropics, where the density is reduced. Hence the pressure at higher 
altitudes must be much greater in the tropics than at the poles—unless 
there is a large difference the other way at the ground. The latter is 
impossible, however, because of the frictional resistance of the ground 
to a strong mean geostrophic winds (this will be considered further 
in § 3). 

A simple numerical estimate of the effect is obtainable from the well- 
known distribution of pressure in an atmosphere in which the temperature 
is independent of altitude. This is 


P=Ppo exp (—gz/RT), ee i 
where py and 7’) are the pressure and temperature at sea level, z is height 
above sea level, and R is the gas constant for air. If p, and 7’ vary with 


Y, which we may take as signifying distance from the nearest pole (to fix 
the ideas), then at fixed altitude 


Op Apo Pogz aT 


This pressure gradient takes its maximum value at the altitude 


Ae —( on Po * dpy/dy 
gy Ty) * aT /dy}’ 
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This result is interesting as showing that, whatever the pressure gradient 


z 


(4) 


where its value is 
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on the ground,* the largest pressure gradient will be at altitude, and 
proportional to the sea-level temperature gradient. The altitude (4) 
at which this maximum pressure gradient is achieved varies little if one 
is dealing with annual mean pressures; the bracketed factor remains 
near 1,} because at sea level the mean temperature varies with latitude 
by a greater fraction of its absolute value than does the mean pressure. 
Under these circumstances, the altitude (4) corresponds to a pressure of 
around 300 to 400mb, such as is found in the upper parts of the 
troposphere. 

The mean geostrophic wind velocity wu corresponding to the mean 
pressure gradient (3) is given by (1) as 


Pee -1.0p 1 (Fe Tap) 


(6) 


Qw sindpdy 2wsingd\p, dy ' Ty dy 
This gives a linear increase of geostrophic wind with height, at a rate 
g 1 dT 5 
Dy sin oT, dy ° . ° ° . . . . (7) 


per unit of altitude. 

The above results are especially valuable because, when the calculations 
are repeated with a realistic value of the lapse-rate (rate of decrease of 
temperature with height), the conclusions are negligibly changed, the 
0-4 (which in (5) represented e~) being retained to one significant figure. 
The change in (7) is even smaller. These results are hardly surprising 
since the temperature at 300-400 mb is less than 20% below its sea-level 
value. 

As order-of-magnitude explanations of what is observed, the results 
(5) and (7) seem quite satisfactory. Equation (5) suggests that the 
change of mean pressure between tropics and poles in the upper part of 
the troposphere (where the pressure is around 300 to 400 mb) should 
be about 40% of the relative change of sea-level temperature, multiplied 
by the sea-level pressure (1000 mb). Taking the relative difference in 
absolute temperature as one-fifth, we obtain 80 mb for the pressure 
change, a slight overestimate. The © sea-level temperatures’ used 
should really be chosen as somewhat less extreme values, however, 
since they ought to represent typical atmospheric temperatures outside 
the thermal boundary layer near the ground. 

An average value of (7), ‘ weighted ’ by a factor sin ¢ (which is reasonable 
since the mean temperature gradient is smaller in lower latitudest), is 


AT 
Oe: (8) 
a De i Se nS 
* Subject to the limitation that the altitude (4) is positive. 


+ This has been used also in (5), the exponential of minus the bracketed 


factor being replaced by 0-4 to one significant figure. 
{ Actually, the theory given (‘thermal wind ’ theory) is less accurate in the 
tropics than elsewhere, since Coriolis force is smaller for small ¢, and so 


dominates the situation far less. 
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where r=6-4 108 cm is the mean radius of the earth (so that 
g/2wr=10-*s—1), 

and AT',/T,, is the relative temperature difference between tropics and 

poles. If the latter is again taken as one-fifth, this estimate of average 

velocity gradient upwards represents an increase of 20 m/sec in 10 km, 

which is of the right order of magnitude. 

As a final remark, it is worth noting that the analysis suggests that 
the kinetic energy of the thermal wind per unit altitude will not reach 
its maximum until a very great height—very roughly, twice that indicated 
in (4). Note that the momentum density pu is greatest at the height (4), 
but that the energy density }pw? is still increasing there, since wu is 
increasing. 

§3. THe MEAN WINDS NEAR THE GROUND 

It was suggested in § 2 that the mean winds near the ground are largely 
prevented from attaining speeds as high as those found at greater altitudes, 
because a high wind near the ground is resisted by such a great frictional 
force. We must now consider how it is that the air near the ground at 
a particular latitude, in the temperate zones for example, can possess any 
mean westerly component of velocity, in spite of the resulting frictional 
couple which must be tending (on the average) to reduce the west-east 
angular momentum about the earth’s axis of the air at that latitude. 
The gravitational and pressure forces on the belt of air at the latitude 
in question cannot alter its angular momentum about the earth’s axis, 
since they all act through that axis. Hence the angular momentum 
loss can be made up only by a steady mean inflow of angular momentum 
from the lower latitudes, where the mean winds at the ground have an 
easterly component, so that in those latitudes the frictional couple 
constantly creates west-east angular momentum rather than destroying 
it. 

There is no conjecture in this part of the theory. We can be certain 
that in the temperate zones, where the mean geostrophic winds near the 
ground are westerly (veering towards 8.W.* in the frictional layer, since 
therein the Coriolis force is reduced but the pressure gradient is unaltered), 
the inflow of west-east angular momentum from lower latitudes must 
exceed the outflow towards higher latitudes. Conversely, in the trade- 
wind zones, where the mean geostrophic winds near the ground are 
easterly (veering towards N.E.+ in the frictional layer, for the same 
reason as before), the outflow of west-east angular momentum towards 
higher latitudes must exceed the inflow from lower latitudes. This is 
because the friction of the ground is in such a direction as to reduce 
easterly winds, and so increase west-east angular momentum. 

Thus, to explain the mean winds near the ground, it is necessary to 
suppose a continual mean flow of west-east angular momentum from the 
ee ee ee ee 


* N.W. in the southern hemisphere. 
+ S.E. in the southern hemisphere. 
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tropics to the poles, a flow which reaches its maximum at around 30° 
latitude. Then the air at latitudes higher than this is continually 
gaining west—east angular momentum, and so can afford to have a certain 
mean velocity in this direction even at the ground, with the resisting 
couple which this entails. The air at latitudes less than 30° is continually 
losing west-east angular momentum. Hence the westerly velocity 
component of this air is reduced, and takes such a negative value at the 
ground that the resulting frictional couple just makes up, on the average, 
for the angular momentum lost. 

To put it more graphically—if friction were far stronger than it is, 
there would be little mean circulation at the ground, and the general 
increase with height described in §2 would be an increase from zero. 
In reality, however, the flow of angular momentum polewards causes 
a general reduction of this circulation in the trade-wind zones, and 
a general increase in the temperate zones, both changes being limited 
by the frictional resistance of the ground. 


$4. THEORIES OF THE PoLEWARD FLow oF West—Hast 
ANGULAR MOMENTUM 


The considerations of §3 have made clear what is the last main feature 
needing explanation, if an ‘ elementary reasoned account ’ of the general 
circulation is to be achieved. It is a poleward mean flow of angular 
momentum, reaching a maximum at around 30° latitude, and falling to 
zero at the poles as well as somewhere near the equator. Knowing as 
we do the sea-level mean wind distribution, and a rough approximation 
to the law which relates the frictional resistance of the ground to the 
. wind speeds near it, we can calculate (Priestley 1951) the order of mag- 
nitude of angular momentum flow which is required. The maximum 
to which it rises is about 31026 dyn cm. This corresponds to a mean 
flow of westerly momentum, crossing each centimetre of the thirtieth 
parallel of latitude, of around 1-6x 108 dynes. Such a flow has been 
observed directly by Starr and White (1951, 1952) by analysing records 
at a large number of meteorological stations situated at around 30°N. 

It is this poleward flow of angular momentum concerning which the 
greatest uncertainty and controversy among meteorologists eXISts. 
It appears from their discussions that several mechanisms for it exist. 

One mechanism, which is certainly among those present, and has high 
suggestive value, was put forward in 1936 by the great aerodynamicist 
Ludwig Prandtl, who reproduces it ( 1952) in his book The Essentials of 
Fluid Dynamics. The small attention which it has received among 
meteorologists may perhaps be due to Prandtl’s undoubtedly too rigid 
ignoration of other mechanisms. ee 

This mechanism must be mentioned, however, because it is concerned 
with the process by which the general circulation in the upper troposphere, 
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considered in §2, is maintained against any energy losses that may 
occur. The mean circulation loses kinetic energy by a variety of causes, 
most of which can be summed up in the phrase ‘the tendency of dis- 
turbances to the mean circulation (such as turbulence, cyclones, and 
even annual changes like the great continental anticyclones) to equalize 
the momentum of different parts of it, and hence reduce the energy ’.* 
How does the mean circulation regain this energy ? Evidently, it can 
do so only if the poleward forces (namely, pressure gradients) which 
produce it are allowed to push the air a bit.t Thus, if dissipation is 
taken into account, the mean winds in the upper troposphere can no 
longer be taken simply along the isobars, but must be imagined to have, 
at least sometimes, a component down the pressure gradient. 

Now, these velocity components down the pressure gradient must tend 
to generate a poleward flow of west-east angular momentum. For the 
work done by pressure gradient is (0p/dy)v per unit volume, where v is 
the velocity down it.{ The flow of angular momentum is proportional 
to (pw)v per unit volume, provided that the isobars in the upper troposphere 
lie roughly in the east-west direction. Since by eqn. (1) for the geos- 
trophic wind, pu is proportional to dp/dy, any north-south movements 
that put energy into the general circulation must also give a poleward 
flow of angular momentum. In the trade-wind zone, for example, 
the upper air moves polewards and the lower air moves equatorwards, 
and this transfers angular momentum because the upper air has the 
stronger westerly circulation. 

This is the only place in the theory where any secondary north-south 
circulation (‘due to hot air rismg at the equator’, a circulation which 
played such an important part in the older accounts of the subject) is 
mentioned. Note however its completely secondary function, which 
is merely to maintain a slight preponderance of west-east angular | 
momentum in higher latitudes, so that the prevailing winds near the 
ground may have westerly components in higher latitudes and easterly 
components in lower latitudes. Near the equator, the process can 
reasonably be described as ‘ hot air rising’, for there Coriolis forces are at 
their weakest, so that the tendency for air to move down a pressure 


* Jeffreys (1926) made calculations which led to the rejection of a mechanism 
similar to that of Prandtl. The reason, however, why the Jeffreys mechanism 
gave insufficient flow of angular momentum, was that only the direct action 
of turbulence (in the frictional layer) on the mean circulation was taken into 
account as a means of reducing its energy. Now, although turbulence is the 
final translator of kinetic energy into heat, the large-scale disturbances form 
a most important intermediate stage in the degradation process. This was 
recognized by Jeffreys in a later paper (1933), which gives a theory more 
closely similar to that put forward here. : 

+ They are enabled to do so, of course, by the fact that, if the geostrophic 
wind is weakened from any cause, the air can flow down the pressure gradient, 
which then exceeds the opposing Coriolis force. 

tie. in the negative y-direction. 
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gradient must be at its strongest. The reason why even equatorial 
air, at the great heights to which it rises (expanding up a saturated 
adiabatic), is deflected into a westerly direction, even though no Coriolis 
force is available for the purpose, must be that exchange of westerly 
momentum with neighbouring latitudes is continually occurring in the 
upper troposphere. 

However, the secondary circulation certainly does not penetrate to 
30°N, where Starr and White (1951, 1952) found negligible mean north— 
south wind components. There, it is the upper tropospheric air with 
strongest westerly components which moves North, while the upper 
tropospheric air with weakest westerly components moves South ; 
however, this process equally transfers angular momentum polewards 
and puts energy into the general circulation. 

It is possible to compute the mean rate of kinetic energy loss by the 
general circulation which would be required, according to Prandtl’s 
theory, to produce the angular momentum flow which Priestley (1951) 
inferred from the east-west components of wind velocity observed near 
the ground. The average loss required is 5x 10% ergs/sec per square 
centimetre of the earth’s surface. This agrees with a direct estimate 
of atmospheric energy dissipation due to Brunt (1926). Both estimates 
have only order-of-magnitude significance. They imply that the kinetic 
energy of the atmosphere would decay by a factor e~1 in around one day 
if the temperature distribution which maintains the motion ceased to be 
maintained (in turn) by the difference in geographical distribution 
between the incoming solar radiation and the outgoing terrestrial 
radiation. The estimates give an efficiency of around 2% for the 
atmosphere as a heat engine, taking the heat supply to be that portion 
of solar radiation which is not instantly reflected. An engine describing 
a Carnot cycle with a temperature difference of one-fifth of the absolute 
temperature would have an efficiency of about 20%. The atmosphere 
can be congratulated if it has really achieved an efficiency of one-tenth 
of ideal efficiency, under distinctly unfavourable conditions. 

The Prandtl mechanism undoubtedly pays too little attention to angular 
momentum transport by the disturbances to the mean flow—whereas 
the importance of all kinds of disturbances in reducing the kinetic 
energy of the mean circulation is explicitly emphasized in the theory. 
This omission is particularly important in the cyclone belt—both because 
the disturbances there are particularly large, and because, in a region 
where the angular momentum per unit mass is decreasing polewards, 
any disturbances are likely, in smoothing the distribution, to cause 
further transport polewards. Again, the hypothesis of the Prandtl 
theory that on the whole the isobars in the upper troposphere lie in an 
east-west direction is violated more constantly in the cyclone belt than 
anywhere else. ae 

In the lower latitudes, however, including those around 30° (where 
the greatest flow of angular momentum is required), disturbances are 
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probably less important. For any irregular system of disturbances 
would make angular momentum flow towards the region of least angular 
momentum—i.e. towards the tropics. On the other hand, specially 
regular systems of disturbances, which have been postulated by some 
authors, do not seem very likely. 

One kind of disturbance that is doubtless present is a fluctuation in 
the northerly components of both wind and pressure gradient. The 
Prandtl theory is sufficiently elastic, however, to admit fluctuations 
of this kind. At any instant the work being done by pressure gradient 
on unit volume of fluid is its product with the velocity component down 
it. Provided that the isobars in the upper troposphere remain approxi- 
mately in the east-west direction, this product is proportional to the 
poleward flow of angular momentum at that instant. The annual 
means of both energy input and flow of angular momentum at the 
latitude in question would then be approximately in this fixed, known, 
proportion to one another—as was assumed above in obtaining the 
figure of 5 x 10% for the mean dissipation rate. 

Starr and White (1951, 1952) found that in fact there are no significant 
mean north-south motions in the upper troposphere at 30°N. But the 
strong northerly components in that region occur on the average where 
strong westerly components are also present, while strong southerly 
components coincide with weak westerly components. Evidently, the 
mean effect of this is both to replenish the energy of the general circulation 
and to produce poleward flow of angular momentum. 

To sum up, the required flow of angular momentum in the tropics and 
sub-tropics can largely be achieved by means of such northerly winds 
down the pressure gradients in the upper troposphere as are required to 
maintain the energy of the westerly circulation in that region. It is 
probable, however, that in the cyclone belt an important amount of 
transport is achieved by the cyclones themselves. 
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ABSTRACT 


The specific heat of liquid helium has been measured at several densities 
over the temperature range 1° to 2°k. The values obtained under the 
saturated vapour pressure are about 10° higher than those observed 
by earlier workers; the results are briefly discussed. 


$1. INTRODUCTION 


Previous workers have shown that the application of pressure to liquid 
helium produces large changes in the position of the A-point, the thermal 
conductivity, the viscosity and the velocities of first and second sound. 
Apart from some measurements in the immediate region of the A-point 
there are, however, hardly any values for the specific heat under pressure. 
The present measurements were therefore undertaken to obtain more 
detailed information, particularly in the range 1° to 2°K where most of 
the other experiments have been performed. Besides the measurements 
at higher pressures, we have also obtained values along the liquid—vapour 
equilibrium curve ; these are about 10% higher than those given recently 
by Kramers, Wasscher and Gorter (1952). | Although we cannot suggest 
any reason for this discrepancy, we are publishing the present results as a 
first step in clarifying the position. 


§ 2. THe EXPERIMENT 


A schematic diagram of the apparatus is shown in fig. 1. The vessel A 
held liquid helium under atmospheric pressure and served as a thermal 
shield which allowed the cryostat B to be cooled to about 1-0°K without 
undue loss of liquid helium. This cryostat was so constructed that 
liquid helium totally surrounded the thick walled calorimeter C which 
contained the mole or so of liquid helium on which the measurements were 
made, The calorimeter also contained a long resistance heater which 
was non-inductively wound, and a carbon film thermometer deposited 
on a perspex former. The necessary electrical leads were brought out 
of the calorimeter through seals and were then brought into good thermal 
contact with the cryostat before passing up to higher temperatures. 
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Although standard methods of low temperature calorimetry were 
employed, the column of liquid helium, through which the pressure was 
applied to the calorimeter, necessitated some precautions in the design of 
the apparatus. This capillary, D, was made as narrow as was convenient, 
in order toreduce the heat flowalong it ; thesections between the hydrogen 
bath E and vessel A and between vessels A and B were each about 20 cm 
long and of 0-3 mm I.D., while that between the cryostat and the calori- 
meter was 10cm long and of 0:15mm I.D. Within the two baths 
A and B, the german silver capillary was replaced by copper tubing of 
greater diameter (2mm I.D.) which passed right through the liquid 


Fig. 1 
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Schematic diagram of the calorimeter and cryostat. 


helium in the baths. This arrangement was designed to reduce the 
temperature difference between the liquid helium in the tubes and that 
in the surrounding baths. A calculation, based on the results of Keesom, 
Saris and Meyer (1940) for the heat flow along a capillary of liquid helium 
and on the results of Kapitza (1941) for the thermal resistance of a 
helium—copper interface, shows that in the case of the bath B the tem- 
perature difference should be negligible ; this conclusion was confirmed 
by a subsidiary experiment with and without helium in the capillary. 
Were it not so, the temperature of the specimen and the thermometer 
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would not have been the same as that of the bath, and the thermometer 
calibration would have been inaccurate. Good thermal contact between 
the helium in the capillary and the bath B also ensured that the liquid 
in the capillary above B was thermally isolated from that in the 
calorimeter. 

Particular attention was paid to obtaining as steady temperature 
conditions as possible, and the following experimental procedure was 
adopted. The temperature of the cryostat was maintained at some 
convenient value for about 30 minutes to allow the calorimeter to attain 
thermal equilibrium. The resistance of the carbon thermometer was 
then measured and the vapour pressure of the helium in the cryostat 
noted so as to obtain a calibration point for the thermometer. (Below 
1-4°, the temperature of the cryostat was measured more accurately by a 
magnetic susceptibility thermometer calibrated against the vapour pres- 
sure at the higher temperatures.) Heat was then supplied for a period 
of about 20 seconds by a stop-watch controlled circuit and observations 
made of the new temperature of the calorimeter. As the capillary of 
liquid helium acted as a comparatively large thermal link, the calorimeter 
drifted back towards the cryostat temperature and this drift rate was 
observed in order to make the usual extrapolation. In all cases save for 
one or two points near 1-0° the extrapolation was linear and introduced 
only a very small error. 


$3. THe Speciric Heat 


While making a measurement of the specific heat the helium in the 
calorimeter was isolated by closing off the top of the capillary; thus, 
after a small correction for the dead volume of this capillary, our results 
gave the specific heat at constant volume, C,. As a matter of con- 
venience we chose to maintain the pressure over the helium at some 
constant value during each run, only isolating the calorimeter during an 
actual measurement ; hence each set of experimental points consists of 
values of C,, at densities corresponding to a given pressure. Our results 
obtained under the saturated vapour pressure and under one atmosphere 
pressure are shown in fig. 2; as is to be expected the difference between 
the two sets of values of C, is less than the experimental error. The 
earlier results of Keesom and Clusius (1932) and Keesom and Keesom 
(1932 and 1935) for the specific heat under the saturated vapour pressure 
are also shown, corrected to the 1949 temperature scale ; in the region 
with which we are concerned their points are rather scattered. The 
measurements of Hull, Wilkinson and Wilks (1951) were primarily con- 
cerned with the temperature region below 1°K and were carried out using a 
closed capsule technique not very suitable for measurements at higher 
temperatures ; the experiments of Keesom and Westmijze (1941) have 
never been reported in any detail. Thus the most reliable of the earlier 
measurements in this region would appear to be those of Kramers, 
Wasscher and Gorter (1952), and it is a matter of concern that our values 
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are about 10° higher than theirs. We became aware of this discrepancy 
at an early stage in the measurements and checked over our work in 
some detail but could find no reason to modify our values. We estimate 
our accuracy at +2°% ; Kramers, Wasscher and Gorter do not estimate 
theirs but the scatter on their results suggests that it is about the same 
or a little less. Thus, 5°% discrepancy between the results would not call 
for special comment, but 10°% is more difficult to account for. Figure 3 
shows the only relevant results obtained previously under higher pressures — 
—a few scattered points due to Keesom and Clusius (1932) at ‘ about 
19 atm’; two of these are actually higher than our values for 19:8 
atm, which are also shown. 


Fig. 2 
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Values of C’, may also be obtained, far less directly, from the velocity 
of second sound together with observations of either the torque on a 
Rayleigh Disc due to waves of second sound (Pellam and Hanson 1952) 
or the motion of a pile of plates oscillating in the liquid (Andronikashvili 
1946). Here again the evidence is inconclusive quite apart from any 
questions concerning the validity of the two fluid model which is used 
in the calculations. Figure 4 shows values of the specific heat under the 
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_ vapour pressure deduced by Pellam and Hanson (1952) from Rayleigh 
Disc experiments. These points agree better with our results than with 
the Leyden ones, but the scatter on them indicates that the experiment 


Fig. 3 
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The specific heat of liquid helium (C,). 
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Fig. 4 
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The specific heat of liquid helium under the saturated vapour pressure. 
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is a delicate and difficult one. The oscillating disc experiment is also 
difficult (Andronikashvili 1946, 1948, Hollis-Hallet 1952), and the absolute 
values of the density of the normal fluid, p,, obtained from it are accurate 
to only about 10% because of uncertain geometrical corrections. 

Using our results, together with the isopycnal data of Keesom and 
Keesom (1933), we have constructed table 1 which SIE smoothed 
’ interpolated values of C’, for five different densities. (These densities 
correspond roughly with the pressures at which sets of measurements 


Table 1. The Specific Heat of Liquid Helium, C, 


C, (cals g—} deg.—*) 


p (g cm-%) 0-1465 0-1525 0-1630 0-1675 0-1715 
GPT 
1-10 0-052 0-070 
1-15 0-068 0-090 
1:20 0-089 0-099 0-11, 0-13, 
1-25 0-11, 0-12, 0-14, 0-17, 
1:30 0-13, 0-15, 0-18, 0-21, 
1-35 C1, 0-18, 0-22, 0-26, 
1-40 0:20, 0-21, 0-23, 0:27, 0-32, 
1:45 0-254 0-25, 0-27, 0-32, 0:38, 
1-50 0-29. 0:30; 0-34, 0-39, 0-47, 
1-55 0-355 0-365 0-40, 0-46, 0-545 
1-60 0-41, 0-46, 0-48, 0-54. 0-65, 
1-65 0-48, 0-50, 0-57, 0-63; 0-505 
1:70 0-569 0-58, 0-68, 0-75, 0-87; 
1-75 0-64, 0-67, 0-79. 0-88, 1-03 
1-8¢ 0-745 0-78, 0-93; 1-04 
1-85 0:86, 0-90; 1-08 
1-90 0-99, 1-04 1-21 
1-95 1:14 1-19 
2-00 1:30 
2-05 1-48 


were made.) The increase of specific heat with density may of course 
be deduced from Keesom and Keesom’s isopyenal data. This is not a 
very satisfactory procedure as it involves taking the second differential 
of a smoothed curve, but to within the possible accuracy the calculated 
pressure dependence is the same as we observe. This isopycnal data has 
also been used to obtain approximate values of the difference between 
C, and C', for the five densities mentioned in table 1; as may be seen in 
table 2, the difference is small except in the region near the A-point. 
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Table 2. C,—C, expressed as a percentage of C,, 


Pressure of Helium 


eK 1 Atm 5 Atm 15 Atm 20 Atm 25 Atm 
1-5 0-01 0-2 1:0 2-1 3:3 
1-6 0:03 “0:3 1:8 3-0 5-9 
1-7 0-05 0-5 2-9 2 9-3 
1:8 0-09 0:7 4-4 8-0 

1-9 0-14 1-0 10-1 

2-0 0-17 


§ 4. THe ENTROPY 


It is, of course, not possible to obtain absolute values of the entropy 
from our measurements because they do not extend to the lowest tem- 
peratures. In the case of helium under its saturated vapour pressure 
one could take the entropy at about 1-2°K to be that given by Kramers, 
Wasscher and Gorter (1952) but in view of the discrepancy between the 
specific heat results this hardly seems justifiable. However, by plotting 
our values of the specific heat logarithmically it may be seen that below 
1-3° the temperature dependence follows a 7'* law, as do the Leyden values. 
We have therefore assumed that this relation may be extrapolated down to 
0-6° where the entropy has become very small (Hull, Wilkinson and Wilks 
1951). In this way we obtain values of the entropy which are about 
10% higher than those of Kramers, Wasscher and Gorter over the entire 
temperature range. A similar extrapolation has been made at the 
higher pressures and table 3 shows values of the entropy at five different 
densities obtained by a numerical integration of C,,. 

These new values of the entropy would explain two sets of experimental 
observations which have not hitherto agreed with theory. Kapitza 
(1941) measured the fountain effect in liquid helium and observed greater 
pressure differences than those predicted by the relation given by H. 
London (1939) 4p/47’=pS ; that is, the values for the entropy which he 
deduced from his experiments were greater than those calculated from 
the specific heat measurements. Kapitza also measured the heat of 
transport by observing the rate of flow through a superleak into a vessel 
to which heat was supplied. Assuming the relation from the London 
theory that this quantity of heat @ should be equal to 7'S, he deduced 
further values for the entropy in agreement with his previous ones. 
Recently Chandrasekhar (1952) has repeated the latter experiment 
replacing the superleak by a link of helium film ; his results are in agree- 
ment with Kapitza’s. Although there has been much discussion of the 
validity of the London relations, the fact that both types of experiment 
(on fountain pressure and heat of transport) lead to the same value for 
the entropy suggests that we are dealing with reversible processes for 
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which they should apply. In any case it is not at all clear how both 
experiments could give rise to too great a value for the entropy. The 
discrepancies are however removed by taking our new values for the 
entropy as is shown in fig. 5. It has however to be noted that the results 
of Meyer and Mellink (1947) and van den Meij denberg, Taconis, Beenakker 
and Wansink (1954) for the fountain pressure lead to somewhat lower 
values for the entropy; unfortunately the scatter on some of their 
points is of the order of 10% and it is hard to draw any firm conclusions. 


Table 3. The Entropy of Liquid Helium 


S (cal g- deg.-) 


p (g em?) 0-1465 0-1525 0-1630 0-1675 0-1715 
TSK 
1-20 0-014 0-016 0-019 0-022 
1-25 0-018 0-020 0-025 0-028 
1:30 0-023 0-026 0-031 0-036 
1-35 0-029 0-033 0-039 0-045 
1-40 0-036 0-036 0-040 0-048 0-056 
1-45 0-044 0-044. 0-049 0-058 0-068 
1:50 0-053 0-054. 0-059 0-071 0-083 
1-55 0-064 0-064 0-072 0-084 0-099 
1-60 0-075 0-077 0-086 0-105 O-ll, 
1-65 0-088 0-091 0-10, 0-11, 0-14, 
1-70 0-10, 0-10, 0-12, 0-13, 0-16, 
1-75 0-12, 0-12, 0-14, 0-16, 0-19, 
1-80 0-14, 0-14, 0-16, 0-195 
1-85 0-16, 0-16, 0-19, 
1-90 0-18, 0-19; 0-22, 
1-95 0-21, 0-22, 
2-00 0-24, 
2-05 0:28, 


One further consequence of accepting higher values for the entropy 
would be that some recalculation of the absolute temperature scale would 
be called for but as the changes will be very small we see no point in 
discussing this in detail until there is a generally accepted value for the 
specific heat, 

§5. Lanpavu’s THEORY 


Finally we may note that these results offer, at least in principle, an 
interesting check on Landau’s theory of liquid helium (Landau 1941, 
1947). Such confirmation as this theory has so far received has come 
from observations on the liquid under its saturated vapour pressure, 
whereas the theory should presumably apply equally well to liquid helium 
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II of any density. As many properties of the liquid are very sensitive 
to pressure, it would therefore be instructive to see how far the theory 
can account for their behaviour at higher densities.* However, the only 
relevant results previously obtained under pressure have been the velo- 
cities of second sound (Maurer and Herlin 1950) and these by themselves 
were insufficient to check the theory. At first sight it appears that by 
combining these results with our specific heat measurements one could 
both check the theory and evaluate, for each density, the three parameters 
(A p p,) in Landau’s equation for the roton energy spectrum. 

Landau considers helium at some fixed density and gives expressions 
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for the specific heat and the density of the normal component, p,, in 
terms of the three parameters 4, 4 and py. It so happens that the tem- 
perature dependence of both the specific heat and Py are determined 
only by the parameter 4. We may thus determine 4 by two indepen- 
dent methods, either directly from the temperature dependence of the 
specific heat or by deriving pp, from the observed values of the specific 


*TIn this connection it should be noted that as the difference between Cp 
and C, is not negligible at the higher densities, it is important to specify 
carefully the conditions under which measurements are made. 
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heat and the velocity of second sound and then determining its tempera- 
ture dependence. The other two parameters, 1 and py»), may then be 
found from the absolute magnitudes of the specific heat and of p,. 
However, if we attempt to determine 4, and pp for helium at a given 
density we find that the values of 4 deduced from the temperature depen- 
dence of the specific heat and from p, are not in complete agreement. 
It is true that this discrepancy is only about 5%, but as 4 enters into the 
expression for the specific heat as a high exponential power, a 5% varia- 
tion corresponds to about a 50% variation in specific heat. Thus, to 
decide how good an approximation to the truth Landau’s roton formula 
is, we require values for the specific heat, entropy and velocity of second 
sound, which are accurate to about 0-1%. Any further investigation of 
this point—and also of the previously discussed question as to whether 
the London relations for the fountain pressure and allied phenomena are 
really exact—would require much more accurate measurements of all 
the quantities concerned than are yet available. 
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ABSTRACT 


The surface energy of sodium, calculated by the use of the free-electron 
model, is modified by considering the increased electrostatic energy in the 
surface if the positive charge is assumed to be localized in the planes of the 
atomic nuclei instead of distributed uniformly. While the theoretical 
surface energy is increased significantly, the theory does not yet account 


_ for the energy observed experimentally. 


SEVERAL attempts to calculate the surface energy of metals by the use of 
the free-electron model have been made (Brager and Schuchowitzky 1946, 
Huang and Wyllie 1949, Huntington 1951). Such calculated energies 
should most nearly agree with experimental observation for the alkali 
metals, for which the Fermi energy corresponds closely to that calculated 
for free electrons. However, the most recent of the surface energy treat- 
ments, that of Huntington, yields a value of surface energy for sodium 
treated as a free-electron metal which is only about one-fourth of the 
observed value. 

Huntington used the surface barrier computed by Bardeen (1936), and. 
evaluated the phase shifts for electrons of various energies up to the Fermi 
limit which occur upon cutting a new surface in the metal. He found the 
surface enthalpy to be 88 erg/em?; at the melting point the surface free 
energy was found, by the use of the surface entropy calculated by Huang 
and Wyllie (1949), to be 54 erg/cm?. The charge distribution in the 
surface was that calculated by Bardeen (1936), and is shown in fig. 1 (a). 
The positive charge is assumed to be distributed uniformly right up to the 
geometrical surface, beyond which its density is zero, while the electronic 
charge distribution, calculated from the free-electron theory, decreases 
more gradually and extends well beyond the surface. Charge is conserved, 
that is the areas under the positive and negative distribution curves are 
equal. The density of electronic charge in the interior is — po. 

It may be seen that a part of the surface energy is electrostatic and arises 
from the charge distribution in the surface as compared to the interior. 
Let this part be H,. The remainder, #’,, is the result of phase shifts and 
enna ey ee eee SF 


* Communicated by the Author. 
+ Now at Department of Chemistry, Oregon State College, Corvallis, Oregon, 


UD S.A. 


1174 Allen B. Scott on the 


changes in exchange and correlation potentials for the electrons caused by 
the introduction of the surface. ; 

Even though the positive charge is not distributed uniformly in an 
actual metal, we may, owing to the insensitivity of the electronic wave 
functions to the periodic field of the ions, take Bardeen’s negative charge 
distribution and H,, calculated from Huntington’s result, as essentially 
correct ; we are then free to adopt a model of the surface more nearly 
representative of an actual metal to compute the electrostatic part of the 
surface energy. The model used in the calculation is illustrated in fig. 
1 (0). 

Fig. 1 


-@a O a 
x 


Charge density and potential near the metal surface. 


We suppose the charge density to be a function only of x, the direction 
perpendicular to the surface. The negative charge is distributed 
according to Bardeen’s calculation. Instead of taking the positive charge 
to be uniform, as in Bardeen’s model, we consider it to be uniformly 
distributed in the planes passing through the nuclei of the metal atoms. 
For the charge to be conserved, the density of positive charge per unit 
area in the plane must be ap), where a is the interplanar spacing. For the 
(110) surface in sodium, a is 5-73 atomic units ; Po 18 3°80 x 10-3 atomic 
units. 

The negative charge density curve from Bardeen’s figure was tested 
by determining the moment of the double layer, when the positive charge 
was distributed as in fig. 1 (2). This was done by integrating Poisson’s 
equation numerically from a point outside the metal into the interior. 
Minor adjustments were made in the curve until the moment was the same 
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as that obtained by Bardeen (0-4v). The negative density was brought 
to —py at x=—ja. To maintain the same moment, the first plane of 
positive charge in the non-uniform distribution used here must be located 
3@ units from the geometrical surface of Bardeen. 

We need first H, for the case of uniform positive charge distribution so 
that #,, can be found from Huntington’s result. Then we need #, for 
the case in which the positive charge is distributed in planes. This, plus 
E,,, will give the total surface energy for the second case. 

The potential, 4, was found in each case by numerical integration of 
Poisson’s equation. Since the negative and positive charges are equal, ¢ 
and d¢é/dx may be taken as zero at large values of x, and the integration 
carried toward the interior. 

In the case of the uniform positive charge distribution the electrostatic 
energy per unit area of the surface layer is 


B=1| pod de, 


where p, is the actual charge density, that is, the algebraic sum of the 
positive and negative densities. The factor } enters to avoid counting 
interactions twice. Since p,=0 for x<—4a, the electrostatic energy in 
the interior is zero and the integral above gives simply the electrostatic 
part of the surface energy. JH, for this case was found to 4:7 10~° 
atomic units, or 37 erg/cm?. H, is the difference between this value and 
88 erg/em?, or 51 erg/cm?. 

For the second case we expect that there will be a large electrostatic 
energy in the interior as well as in the surface ; H, is the difference between 
that in the surface and that in the interior. Figure 1 (c) shows the varia- 
tion with x of the potential obtained from the integration of Poisson’s 
equation. The potential at the plane of positive charge nearest the surface 
is ¢,; at a plane of positive charge in the interior the potential is ¢,. 
The energies per unit area of positive charge in the surface layer and in an 
interior layer are apy ¢, and apy ¢; respectively. 

The energy of the negative charge, per unit cross section, in an interior 
atom is 


and in a surface atom is 
oe) 
= | pp dx, 
—ta 


where —p is the variable electronic charge density in the surface atom. 
The excess energy per unit area of surface is then 


B.=3| ap, (6.8) — | pd det [ed ar] 


where again the factor } enters so that interactions are counted only 
once. ¢, was found to be 0-444 and 4, to be 0-424. The energy of the 
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negative charge was found to be —3-76 x 10~° per unit cross section for a 
surface atom and —3-53 x 10-3 for an interior atom. All values are given 
in atomic units. Z, for this case is 1-0 x 10-4 atomic units, or 77 erg/em?. — 

This quantity, added to #,, gives a total surface energy of 128 erg/em?. 
The product of the absolute temperature at the melting point times the 
surface entropy calculated by Huang and Wyllie is 34 erg/cm?, so that the 
surface free energy at the melting point is the difference, or 94 erg/cm?. 

While this value is larger by 70°% than that calculated by Huntington, 
it is still well below the experimental value given by Poindexter and 
Kernaghan (1929), who found the surface energy to be 206 erg/cm? near 
the melting point and to be almost independent of temperature. The 
small temperature coefficient may mean that the calculated surface entropy 
is too large. 

It is evident that even when we consider the effect of the non-uniformity 
of positive charge on the electrostatic energy at the surface, the free- 
electron theory does not yet quantitatively account for the observed 
surface energy of sodium. It is illustrated here, however, that the electro- 
static energy arising from the interaction of the surface ions with the non- 
uniform electronic charge distribution is a very important part of the 
surface energy. | 
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SUMMARY 


The analysis contained in an earlier paper on the inhibition of thermal 
convection by a magnetic field is extended to include the case when the 
impressed magnetic field (H) acts in a direction different from that of 
gravity (g). A variational method is described for determining the critical 
Rayleigh number for the onset of convection of any preassigned pattern ; 
and the result is established that when H and g act in different directions 
the convection which appears at marginal stability is in the form of rolls 
extended in directions parallel to the plane containing H and g. 


$1. IyTRODUCTION 


In an earlier paper (Chandrasekhar 1952: this paper will be referred to 
hereafter as I; see also Thompson 1951) the instability of a horizontal 
layer of fluid heated below has been examined in the case when the fluid 
considered is an electrical conductor and a uniform external magnetic — 
field is impressed on the fluid. In developing the general theory in I no 
restriction was made regarding the direction of the impressed magnetic 
field (H) relative to that of gravity (g). But in deriving the principal 
results on the critical Rayleigh numbers for the onset of instability, the 
assumption was made that the directions of H and g coincide. The case 
when these directions do not coincide was considered only very briefly : 
it was, in fact, dismissed with the remark that if one restricted oneself to 
the onset of convection in the form of infinitely extended rolls (in directions 
parallel to the plane containing H and g) then the critical Rayleigh 
numbers derived for the case when H and g are parallel can be applied to 
the more general case (when H and g are not parallel) if in the relevant 
formulae we interpret H (=| H | in the cases to which the analysis strictly 
applies) to mean the component of H in the direction of g. The question 
whether lower Rayleigh numbers can be reached by considering more 
complicated patterns of motion than rolls was left open with the conjecture 
that this is unlikely. In this paper we return to the question which was 
left open and prove what was conjectured by developing a method which 
enables one to determine the critical Rayleigh number for the onset of 
convection of any preassigned pattern. 


* Communicated by the Author. 
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§ 2. THe REDUCTION TO A CHARACTERISTIC VALUE PROBLEM 


In the notation of I, §7 (eqns. (117) and (118)) the basic equations are 


020 826\— p2Ho a\2 
vivn—2 (554+ 5a) +" aa “(cos # 5 srineg) oy vo ghee 


Bly (O2w  2w\  p?H?o 0 5 
and V®w= By ax? + aye ai a cos sf ;, sin Ee on “y OP 5 
where it may be recalled that the orientation of the co-ordinate axes has 


been so chosen that H lies in the wz-plane and 9 denotes the inclination of 


H to the vertical. 
If we measure all linear distances in units of the depth of the layer, d, 


and let 


pa lb ly ga pes Tes 


and od? ,. "Se 
KV Pov 
eqns. (1) and (2) take the forms 
7-0" 070 0 \2 
ayes 2 
VG. “d (Sa + Fa) +2 (5. +tan da, Jw - ») oe 
w w 0 : ; 

and Vew=R ( — =) +@ @ +tan > =) V202 ee (5) 


Corresponding to an analysis of the disturbance in the horizontal plane in a 
double Fourier integral, we seek solutions of eqns. (4) and (5) which are of 
the forms 


w= W(2) exp [i(a,v+agy)] | (6) 


and 6=0 (2) exp [i(a,2-++-ayy)]. 
where a, and a, are the wave numbers (in the unit 1/d) of the disturbance 
in the x- and the y-directions and W and @ are functions of the vertical 
co-ordinate, z, only. It should be particularly noted that W (z) and @ (z) 
will in general be complex. 

With w and 6 given by (6), eqns. (4) and (5) reduce to 


[(D*—| a P—Q(D-+iw)"| WX a | alt @> 4.6 cane 


and (D?—| a P)[(D?—| a P?})?}—Q(D+iw)*]W=—R|aPWw, . (8) 
d 
where Da aq? JaP=a’+a2 and w=a,tan% . . (9) 


The boundary conditions with respect to which eqns. (7) and (8) must be 
solved are (cf. I, eqns. (53) and (54)) : 


| W |=| 6 |=0 for z=0~ and'71 1G} 
and |DW|=0 ona rigid surface, 7 
|D°W|=0 ona free surface. ; a 
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In view of eqn. (7) the boundary conditions (10) are equivalent to 
| W [=| {(D2—| a [?)?—Q(D+ia)}W |=0 forz=0 and 1. (12) 


Since W is complex, eqn. (8) together with the boundary conditions (11) 
and (12) constitute a characteristic value problem in an equation, 
genuinely, of order twelve. However, eqn. (8) being linear and one with 
constant. coefficients, the problem presents no difficulty of principle. But 
as the solution of the physical problem requires the determination of the 
Minimum value of R as a function of the two variables |a| and 
ow (<| a|tan %) for various assigned values of Q, it is clear that the 
practical problem is a very formidable one if the characteristic value 
problem has to be solved in any direct manner. Fortunately, it will 
appear that this is not necessary, since, as in I, a more convenient method 
which is practicable can be devised. This method which is based on a 
variational principle is described in the following section. 


§ 3. THE VARIATIONAL METHOD 


Letting 
F=[(D*|aPP—QD+io)W, . . . . (13) 
we first rewrite the equation governing W as 
(D?—|aP)F=—R|aPw. ey el Le) 
The boundary conditions (12) now require that 
ay Omeetorz—O0)-andel. os 4. | (15) 


Now multiply eqn. (14) by /*, the complex conjugate of F’, and integrate 
over the range of z (i.e. from 0 to 1). The left-hand side of the equation 
gives 

1 1 
| F*(D2—|a P)F de=[F*DF]— { (|DFP+apl FP) de. (16) 
0 0 
The integrated part vanishes on account of (15) and we are left with 
1 1 
I F*(D*—|aP)F de=— | (| DFP+aP]FpP)de . (17) 
0 0 
Considering next the right-hand side of eqn. (14) we have (cf. eqn. (13)) 


{ * WE* de= | ; W[(D?—| a [2)?—Q(D—iw)?|W* dz. —. (18) 
0 0 


After two successive integrations by parts, the first of the two integrals on 
the right-hand side of the foregoing equation becomes (cf. I, eqns. (65) and 


(66)) F 5 
[ wwrla pew d= | | (Dla PW Pde emake 
0 “0 


4AH2 
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while the second of the two integrals can be transformed in the manner 
| * W(D—iw)W* dz 
a ii W(D2W*—w2W*—2iwD W*) de 
0 
=—[ (DW P+0%| WP) de—2iee | WDW* dz 


—— | (| DW P-+-a| W Pie WDW*+iwW*DW) dz 
0 


=~ | |DW+iwW Pde, pO Ae 
Thus, we have 

ir WF* de= i; | (D—laP)WP+Q| DW+ieW Pd: . (21) 
The result of nate aerer eqn. (14) by F* and integrating over the range of 


is therefore 1 
4 | (DF P+laPl ep) @ 
lifes “ (22) 


See Ss . 
la Py (D2—|a2)W P+Q| DW+ioW -] dz 


This formula expresses R as the ratio of two positive definite integrals ; 
and as in I, eqns. (70)-(76), it can be shown that eqn. (22) similarly provides 
the basis for a variational method of solving the underlying characteristic 
value problem. And the method as applied to eqns. (8), (11) and (12) is 
the following : 

Assume for F/ an expression involving one or more parameters A ,, and 
which vanishes at z=0 and 1. With the chosen form for F determine W 
as the solution of the equation 


[(D?—] a ?)?—Q(D+iw)]JW=F, . . . . . (23) 
which satisfies the remaining four boundary conditions on W. Then 


evaluate RF according to formula (22) and minimize it with respect to the 
parameters A ,. 


§ 4. THe DETERMINATION OF R WHEN H AND g ACT IN DIFFERENT 
DIRECTIONS AND FOR THE CasEs (i) WHEN BotH BounpING 
SURFACES ARE Riqtp, (ii) WHEN Botu BounpInG SURFACES ARE 


FREE AND (iii) WHEN ONE Bounpin@ SuRFACE Is RIGID AND THE 
OTHER IS FREE 


(i) Both bounding surfaces rigid. In this case it is convenient to trans- 
late the z-axis so that the bounding surfaces are at z=}. 


As the simplest trial function for F which vanishes at z=-+} we shall 


assume E008 nap ven UN eee an an ae 


This is a trial function with no variational parameters. Nevertheless, 
from our earlier experience with problems of this kind (ef. I, for example) 
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we may be confident that already this trial function will lead to values of R 
which are accurate to within 1%. 
With F given by (24) the equation to be solved for W is 
[(D2—| a ,?)?—Q(D+iw)?|W=cosmz, . . . . (25) 
and the general solution of this equation is 


W=TI, exp (imz)+T, exp (—imz)+ z A; exp (q;2), . (26) 
where Py=3{(a?+| a [?)?+Q( (rt epi | (27) 
P,=3 (77+ a PP+Q(9—w)? 1, 
the q;’s (J=1, ..., 4) are the roots of the quartic equation 
(| aPRaOGtia?) 2. 1. (28) 
and the A,’s (j=1,... , 4) are constants of integration to be determined by 


the boundary conditions on W. 
It can be readily verified that the roots of eqn. (28) are given by 


qy=a-+ie, do=—a-+ic, qg=b—c, Qa=—b—ic, . (29) 


where a=4(+/Q+«,), b=4(1/Q—«,), C=405, 3 « (30) 
=[4{(Q+4| a P)?+1602Q}”?+4(9+4] a 2) (31) 
md a= [H(Ot+-4] a PY L1eog)*=_-HO+4l a pr. 
The boundary conditions to be satisfied by W are 
W=0 and DW=0 [Ole ee, te (32) 


and these conditions require : 


4 4 
2 A; exp (+29))=— 2 A, exp (—3q,)=—1(,—L9) 
j 


and es Bis . (33) 
a Ajq; exp (+3¢)=— aos exp (—39,)=7(1,+14). 
j=1 i 


It is found that the solution of these equations can be expressed in the 


iis =A" = e708 dee 1 a) Kaeo 1's) Y 3} . (34) 
a 
and AeA a es ee ee Ween Deer Y<1,- . (36) 
where 
A=2[(a.?+@Q) cosh «1+(«12—Q) cos «.—(«42+a,") cosh +/Q], . (36) 


X 1=G4({3—s2) [XP (Y3+Y2)—eXP (Ya) ]+-42(4a—Js) 

x [exp (da +43)—eXP (72) ]+-43(G2—Ya) [EXP (Y2+G4)—EXP (Ys)], (37) 
Y .=(42—4s)[€XP (V2 +93)—eXP (Ya) ]+(U3—@a)lEXP (Ys-+4)—EXP (Ze) ] 

+ (q,—2)[exP (a t92)—exP Gs)];  - - - + + - + + + (88) 
and X, and Y, are given by expressions which can be obtained by replacing 


qo and q; in eqns. (37) and (38) by g, and qp», respectively. 
With F and W given by eqns. (24) and (26) we readily find that 


[ (DF P+la PL? P) de=Moe+]aP), ea 59) 
0 
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exp (+34;) exp (—24)) 

qj +n 
Substituting for the q;'s and the A,’s in accordance with eqns. (29), (34) 
and (35) we find after some very lengthy but straightforward reductions 
that the integral on the left-hand side of eqn. (40) can be reduced to the 
form 


Fil 4 
and | WE* de=3(0,+T2)-+7 2 A, . (40) 
0 I= 


1 
| WF* dz=4(, 41.) 
0 


Qa 


ve ((,—L2)(91 Fi th 61-92 F.—fe 62) 


w 


a SPE nh hE) 1—91 Jitfhe Bo—92 52), + - - (41) 
where 
f1=84L(7? +0? —c?)(1+-e% cos c)+2 e* ac sin ¢], 
91=8,[(7?+a?2—c?) e* sin c— 2ac(1+e* cos c)], 
FF ,=2ab(cos c—e~*) cosh b—2(b?+-2c?)(cos c+-e~*) sinh b 
+ 2ab(cos 2c—e-* cos c)— 2bc(sin 2c-+-e~“ sin ¢) 
—2c(a sinh b+-6 cosh b) sin ¢, 


@,=— 2ac(cos c+e~*) sinh b—2bc(cos c—e~*) cosh b 


—2ab cosh b sin c+ 2(b?+ 2c?) sinh b sin ¢ (42) 
—2ab(sin 2c+e-* sin c)—2bc(cos 2c—e~* cos ¢c), 

= — 2a(cos c+e~*) sinh b+ 2b(cos c—e~*) cosh b 
+4c sinh 6 sin c+ 2b(cos 2c—e~‘* cos c), 

¥,=4c(cos c+e~*) sinh b+ 2(a sinh b—6 cosh b) sin ¢ 
—2b(sin 2c+e~? sin c), 

8, =[(a?-+-a2—c?)24 4a2e2}-1, 

F .=2ab(cos c—e~”) cosh a— 2(a?+ 2c?)(cos c+e~°) sinh a 
+ 2ab(cos 2c-—e~? cos c)—2ac(sin 2c+e~? sin c) 
—2c(a cosh a+-b sinh a) sin c, 

6 = — 2be(cos c+-e~*) sinh a—2ac(cos c—e~®) cosh a 
—2ab cosh a sin c+-2(a?+ 2c?) sinh a sin ¢ (43) 


— 2ab(sin 2c-+-+e~® sin c)—2ac(cos 2c—e-® cos c), 
3) .—= — 2b(cos c+e~?) sinh a+ 2a(cos c—e~°) cosh a 

+ 4c sinh a sin c++ 2a(cos 2c—e~? cos c), 
3,=4c(cos c+e~?) sinh a+2(b sinh a—a cosh a) sin ¢ 

— 2a(sin 2c+-e~° sin c), 


8 9=[(m2-+-b®—c2)24 4b%02]-1, 


) 
So=8.[(7?-+b?—c?)(1+e° cos c)+2 e? be sin c], 
J2=5.[(7?-+b?—c?) e” sin c—2bc(1+e? cos c)], 
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Finally, the required value of R (in the chosen approximation) is given 
(apart from the factor | a |?) by the ratio of the expressions on the right- 
hand sides of eqns. (39) and (41). 

Values of R obtained with the aid of the foregoing formulae for various 
values of @, | a| anda are listed in table 1. We return to a discussion of 
these values in § 5, 

(ii) Both bounding surfaces free. With the same trial function for F as 
before (eqn. (24)) the solution for W will again be given by (26). But the 
boundary conditions which determine the constants A j are now different ; 
they are 

Ww=0 and 172 W=-0 for 2=+4. . . «7(44) 


These conditions require : 


4 4 
2 A, exp (+39,)=— 2 A, exp (—3q,)=—t(1,-Iy) 
j=1 j=l , 
: : (45) 
and 2 A,g? exp (+3%))=— moe exp (— 39) =i17?(I;—L). 
j= 


p= 


And the solutions of these equations are found to be expressible in the 
forms 


A,=A,*=1(P\—I5) 


ak 
Sh eo a CO) 


exp (+295) 


Bnd A,=4*=i(P,-9,) 


(Xe7tY.). ©. 1.. (47) 
where 
A=2[(a1?+«,42)Q cosh +/Q—«1?(a.?+@) cosh «1-+%9?(«12—Q) cos x9], 
| (48) 
X 1 =G (73? —12?)[€XP (Ts +2) — EXP (Fa) ]+927(Ya?—45") 
x [exp (da-+43)—eXP (Ya) ]+437(¢2?—Ga”) [EXP (Yo+44)—eXP (3)], 
(49) 
Y,=(q2?—¢3")[€XP (¢2+43) — XP (44) ]+(Ys?—4?)LeXP (43 +44) —eXP (2) | 
+ (¢4?—72")[exp (Va +72) exp (%3)] 5 eee ee wee |. (50) 


and X, and Y, are given by expressions which can be obtained by 
replacing q, and q; in eqns. (49) and (50) by g, and q,, respectively. ; 

The result of substituting the foregoing expressions for the A,’s in 
eqn. (40) (which continues to be valid) is to give 


1 
| WF* dz=3(', +1.) 
0 


ai 3 ((,—T'2){9.(F 1— 7B) + f/1(61—-775) 


—9(F .—7°B.)—fo(G2—7F 2) }; Ee ag eh ED) 
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Table 1. Rayleigh Numbers for Various Assigned Values of @, | a | anda 
for the Case when both Bounding Surfaces are Rigid. (In all 
cases the entry for 7=0 is for a value of | a| for which it is known 
(cf. I, table 2) that the Rayleigh Number attains its Minimum 
Value for Convection setting in as Rolls) 


Q w |a| Rh Q co Ja] R 
sore fare 0 1954 
0-5 1958 
10 1-0 3-25 1972 
0-0 | 40 | 3768 2a) 2026 
39 «| 3797 | 4-0 2007 
0-5 40 | 3793 u 2811 
4-1 3795 = 1-0 he 2876 
3-9 3875 2-0 3072 
1-0 4-0 3878 4:0 3876 
4-1 | 3871 0 10120 
oe 3.9 | 4190 || 500 2-0 5:16 11280 
2-0 4-0 PTW 4 ee pee ee cE ____|__ 15240 
4:1 4178 0 17120 
3.9 5523 1000 3-0 5-80 21580 
4-0) 4-() 5494 6-0 39880 
4:1 5474. Hens 0 si 54710 
2-0 60440 
| 10000 : 8-66 Biaial 


| 4-0 182700 
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where, now 


fi=8,[ (7? +0?—c?)(1+e? cos c)+ 2e% ac sin c], 
91=9,[ (77+ a?—c?) e* sin c— 2ac(1-+-e? cos c)], 
=2BD{(cos 2c—e~* cos c)+ (cos c—e~*) cosh b} 
. +2AD{(sin 2c+e-* sin c)+ cosh 6 sin c} 
+2BC sinh 6 sin c+ 2(C?+ _D?— AC)(cos c+e~*) sinh b, 
6,—=2AD{(cos 2c—e~? cos c)-+(cos c—e~*) cosh 5} 
—2BD{(sin 2c+-e~* sin c)+cosh 6 sin c} ee (D2) 
+2BC(cos c+e~*) sinh b—2(C?+ D?— AC) sinh 6 sin c, 
%, = — 2D {(sin 2c+-e~‘ sin c)+cosh 6 sin c}— 2B sinh 6 sin ¢ 
+2(A—C) (cos c+e-*) sinh b, 
J ,——2D{(cos 2c—e~* cos c)-++ (cos c—e~*) cosh b} 
—2B(cos c+e~*) sinh b—2(A—C) sinh 6 sinc, 
34 =[(n?-pa?—c2)?-4 dare? ]-1, 
| 


fo=5 ol (7?-+b2—c?)(1+-e? cos c)+ 2e” be sin c], 
Jo=8 o[ (7?-+b2—c?) e? sin c—2be(1+e? cos ¢)], 
Ff .—2BD{(cos 2c—e~? cos c)+(cos c—e-®) cosh a} | 
+2BC {(sin 2c+e? sin c)+cosh a sin c} 
+2AD sinh asin c+2(A?+B?— AC)(cos c+e7°) sinha, 
@,—2BC {(cos 2c—e~® cos c)+(cos c—e~°) cosh a} 
—2BD{(sin 2c+e~® sin c)+ cosh @ sin ch | 
+24 D(cos c-++e~°) sinh a—2(A?+ B?—AC) sinh a sin ¢, | 
= —2B{(sin 2c+e~? sin c)+ cosh a sin c}—2D sinh a sin ¢ | 
—2(A—C)(cos c+e-®) sinh a, | 
| 


(53) 


9= —2B{(cos 2c—e~® cos c)+ (cos c—e- >) cosh a} 
—2D(cos c+e~*) sinh a+2(A—C) sinha sine, 
5=[(7?-+b?—c?)?+ 46707]. 


In eqns. (52) and (53) a, 6 and c have the same meanings as in eqns. (42) 
and (43); in addition we have let 


A=a—e, C=b?—c?, B=2ac and D=2be. . (54) 


And again the required value of # (in the chosen approximation) is given 
(apart from the factor | a |?) by the ratio of the expressions on the right- 
hand sides of eqns. (39) and (51). 

Values of R obtained with the aid of the aes formulae for various 
values of Q, | a | and are listed in table 2 
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Table 2. Rayleigh Numbers for Various Assigned Values of Q, | a] anda 
for the Case when both Bounding Surfaces are Free. (In all cases 
the Entry for 7=0 is for a Value of | a | for which it is known (ef. I, 
table 1) that the Rayleigh Number attains its Minimum Value for 
Convection Setting in as Rolls) 


= ja | R Q 7 |a| R 


6-0 3:7 6843 6-0 37820 


3-8 3 0 51520 
yk 4000 7-3] 


2-0 56330 


0 119800 


177200 
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~ (iii) One bounding surface free and the other rigid.t In this case it is 
convenient to have the bounding surfaces at z=0 and 1 and choose for F 


the trial function : 
ieesthas © 4. APCS eee (515) 


The corresponding solution for W is 
. * . D x 
W=—1I, exp (imz)-+iI, exp (—imz)+ 2 A; exp (q;2) ~. (56) 
j=l 


where I,, I’, and the q,'s have the same meanings as hitherto. The 
present boundary conditions on W are 


W=DW=0 | for z—0, Wea D*We-0e tor7— pore) 
These boundary conditions require that 
4 4 
Fee OP et Es): 
= j= 
4 4 6 (58) 
= Jeu 


The solution of these equations are found to be expressible in the forms 


A= GEi-lyXjt 54+ Py¥;, Gabe. 69) 
where 
A=2[ (a 2+a92)1/Q sinh //Q—(Q+a«>”)a, sinh %4—(Q—a«47)a, sin al, 
60 
X 5=9(4a?—{3") CXP (4s+44)—Ya(Ya" —@2") EXP (V2+a) s 
+-¢4(¢s?—42") EXP (dats) — (2? +77) (Za—Ws) EXP (Va) 
+ (qs2-+7)(Ga—Wa) EXP (Vs) — (Ga? +7)(Ya—Wa) EXP (Ya), ++ (81) 
Y =(q1?—43") exp (Ys +44) — (442 —@2”) EXP (Yo) 
+ (322?) EXP (Yas), Ree ee Seen rid 2. 1,905, (62) 


and the expressions for the remaining X’s and Y’s can be obtained by 
suitable permutations of the subscripts of the q’s in eqns. (61) and (62). 
With F and W given by eqns. (55) and (56) we now find that 


1 
[ (DF P+ a PL FP) ded? H1 a) Wei. 4 (63) 
& 4, itexp (g) 

and I, WE de UP4+T)t0 2 Aaa (64) 


J eS ee 

+ In obtaining the solution for this case (when H and g are parallel) in ii 
it was pointed out that on account of symmetry the ‘ odd solutions ’ for the 
case when both bounding surfaces are rigid provide solutions for the case 
when one bounding surface is rigid and the other is free ; but that these latter 
solutions are applicable to a cell depth d' =1d and a Rayleigh number which 
is sixteen times smaller. But it should also have been pointed out that the 
corresponding value of the parameter @ is four times smaller. Hence, the 
values of Q given in the first column in table 3 should be divided by 4 to 
make this table comparable to tables 1 and 2 ; similarly the values of a given 
in the second column should be divided by 2. 
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The reduction of the right-hand side of eqn. (64) with the A,’s given by 
eqns. (59)-(62) is particularly troublesome: it 1s very long and very 
complicated, and requires among other * tricks ’ the repeated use of the 
identities | 
4 i 4 ; 
ee, IT (i7—4;) and ‘am IT (ia-9)) 
j=l 


Ri 2 j=1 


But eventually one finds that 
8a 

{ WE* de=ACAT 2) + Ts 
0 


x [(2(7,-T,) [B(C’-+D?) + D(a?+B?)] +1701 }(Q+a57)a, cosh a, 

— (PI) B(C2-+D*)— D(@P+ BY) +7 Qos? +222) VQ cosh /Q 

am {2([,—I,)[A(€2+ D*®)— €(A? + B*)]+ a9 }(ay?—Q) xg COS %» 

—8(I',—I,) {a[® sin c sinh 6+ cos c cosh 6] 

-Lb[#F sin c sinh a+6 cos c cosh a]}]]. + eee ah tap Oe 

where 

F=(a?—B?)(2,C—B?) + 2BA(€+2)), ! 
6— B{24(2,C — B*)—(a?—B?)(€+2,)}, | (67) 
h=(C— D*)(X,a+ D*?)—2D°€(A—2'), 
J=DLC(EA+D)1(C—DYA-Z,)}, | 


A=n?4 A=n+at—c?; 2y=4(a/Q—a22), | (68) 
J 


1 
2 
C= 4+CaP+P—P 3 2=3(%1V0+a,"), 


anda, b,c, A, B, C and D have the same meanings as in cases (i) and (ii) 
above (eqns. (29)-(31) and (54)). 

Finally, the required value of R (in the chosen approximation) is given 
(again, apart from a factor | a|?) by the ratio of the expressions on the 
right-hand sides of eqns. (63) and (66). 

Values of & obtained with the aid of the foregoing formulae for various 
values of Q, | a | anda are listed in table 3. 


§ 5. THE ONSET OF ConvecTION As LoNGITUDINAL ROLLS WHEN H ANp g 
ACT IN DIFFERENT DIRECTIONS 

We now return to an examination of the results summarized in tables 1, 
2 and 3. 

First, it should be noted that in all cases the entry for o=0 is for a 
value of | a | for which it is known from the results in I that the minimum 
Rayleigh number occurs when the pattern of convection which appears at 
marginal instability is in the form of longitudinal rolls (i.e. rolls extended 
in directions parallel to the plane containing H and g). Thus fora layer of 
liquid confined between two rigid planes, the lowest Rayleigh number is 


ae 
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Table 3. Rayleigh Numbers for Various Assigned Values of Q,|a| anda 
for the Case when One Bounding Surface is Free and the Other 
Bounding Surface is Rigid. (In all cases the Entry for 70 is for a 
Value of | a | for which it is known (cf. I, table 3) that the Rayleigh 
Number attains its Minimum Value for Convection Setting in as 


Rolls) 
Q@ ap |a| R @ wo |a| R 
0 1713 
1-0 1754 
25 3°17 
2-0 1877 
0 5-1 9319 
4-0 2378 
5-0 9389 
0 2231 
0-5 5-1 9387 
1-0 2297 
5-2 9393 50 3°45 
2-0 2496 
5:0 9596 
4-0 3250 
1-0 5-1 9593 — 
0 5628 
5:2 9599 
250 0:5 4-50 5671 
5:0 10456 
500 1-0 5803 
2-0 5-1 10450 — — 
0 16130 
5-2 10453 
a 0-5 16240 
5-0 14396 1000 5°75 
1-0 16580 
4-0) Sal 14371 
4-0 24540 
5:2 14358 || — 
ee eer 0 35060 
5:0 22684 
0-5 35290 
6-0 5-1 22598 2500 6-75 
1-0 35980 
52 22528 
4-0 52490 
0 122200 
10000 8-65 
1-0 125200 
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obtained for | a |=4-00 when @= 100 and convection sets in as longitudinal 
rolls. From the results for Q=100 given in table 1 it follows that for a 
slightly different from zero the minimum Rayleigh number still occurs for 
| a|=4-00; also, that the new minimum slightly exceeds that for a=0. 
For example, for 7=0-:5 and Q=100 the minimum Rayleigh number is 
3793: this should be contrasted with R,=3768 for a=0. For larger 
values of a the minimum appears to shift to somewhat higher values of 
|a|; but the Rayleigh numbers now exceed fk, at o=0 by a substantial 
margin. Thus in the illustrative example considered the minimum 
Rayleigh number for w=1-0 is 3871 and this occurs for |a|=4-1. The 
results for other values of Q@ and other boundary conditions exhibit the 
same general behaviour and confirm the conjecture made in I that when 
H and g act in different directions convection when it first sets in appears 
as longitudinal rolls. 

In some ways the conclusion arrived at in the last paragraph is a 
paradoxical one. For, when H and g are parallel the Rayleigh number 
depends only on | a ?}=a,?+-a,? (cf. eqn. (9)) and we expect convection at 
marginal instability to have a cellular—and, indeed, a hexagonal—pattern. 
How then, one may ask, does this situation change discontinuously when 
H is inclined only very slightly to the vertical and when, it is asserted, 
convection at marginal instability must appear as longitudinal rolls ? 
The resolution of this paradox is as follows : 

When H and g act in different directions, differing Rayleigh numbers 
are required for the marginal appearance of convection with a preassigned 
pattern. In particular, the most ‘ difficult ’ pattern to excite (other things 
being equal) is the system of the transverse rolls. Once the Rayleigh 
number is high enough to excite these transverse rolls, we may expect 
a proper cellular pattern of convection to emerge. The difference in the 
minimum Rayleigh numbers required to excite the longitudinal and the 
transverse rolls is, therefore, a measure of the extent to which a proper 
cellular pattern of convection is suppressed at marginal instability. But, 
and this is the main point of the argument, this difference in the two 
Rayleigh numbers tends to zero as w (=a, tan $<|a| tan 9) tends to 
zero. In other words, when H is only very slightly inclined to the 
direction of gravity, the extent to which transverse rolls are suppressed at 
marginal instability (when longitudinal rolls appear) is also only very 
slight. Finally, when H is exactly parallel to g this suppression of the 
cellular pattern ceases and longitudinal and transverse rolls appear 
simultaneously at marginal instability. 

The picture which we have thus arrived at as to what really ‘ happens ” 
when H and g act in different directions constitutes an essential simpli- 
fication for the general theory of thermal instability. For, it effectively 
ensures that conclusions reached after an examination of special 
geometrical situations (such as H and g are parallel) are of wider scope 
and generality than the geometrical restrictions under which they were 
arrived at would appear to warrant. Thus in the problem of the inhibition 
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of thermal convection by rotation (Chandrasekhar 1953) the principal 
results were derived on the assumption that the direction of the angular 
velocity, Q, is parallel to that of gravity, g. One may now feel confident 
that the conjecture made in that context namely that if Q and g are not 
parallel, convection at marginal instability will be in the form of longitu- 
dinal rolls is a correct one. Similarly, in the discussion of the condition for 
overstability when H and g are acting (I, §8) or when @ and g are acting 
(Chandrasekhar 1953, §8) or when H, ® and g are all acting 
(Chandrasekhar 1954) the assumption that all these vectors are parallel 
is not likely to be a restrictive one. 


In concluding this paper I should like to acknowledge my indebtedness 
to Miss Donna Elbert who carried out all the very laborious numerical 
calculations involved in the preparation of this paper. 

The research reported in this paper has in part been supported by the 
Geophysics Research Directorate of the Air Force Cambridge Research 
Center, Air Research and Development Command, under Contract 
AF 19(604)—299 with the University of Chicago. 


REFERENCES 
CHANDRASEKHAR, S., 1952, Phil. Mag., 48, 501; 1953, Proc. Roy. Soc. A, 
217, 306; 1954, Ibid., 225, 173. 
Tuompson, W. B., 1951, Phil. Mag., 42, 1417. 


fratio2 ] 


CXXXII. Nuclear Orientation of °*Mn 


By M. A. Grace, C. E. Jounson, N. Kurti, H. R. LEMMER 
and F. N. H. RosBinson 
Clarendon Laboratory, Oxford* 


[Received July 3, 1954] 


ABSTRACT 
Orientation of =4Mn at very low temperatures has been attained. 
The anisotropy of the 0-835 Mev y-radiation emitted in its decay shows 
that the first excited state of ®4Cr has spin 2. Detailed interpretation of 
the results is discussed. 


CoNTINUING our experiments on nuclear orientation in paramagnetic ions 
(Bleaney, Daniels, Grace, Halban, Kurti, Robinson and Simon 1954, 
Ambler, Grace, Halban, Kurti, Durand, Johnson and Lemmer 1953) 
we have investigated the 310 day isotope °*Mn which decays by K-capture 
followed by a single y-ray. This particular nucleus was chosen because 
the magnetic properties of the Mn** ion in various paramagnetic crystals 
had been established by paramagnetic resonance experiments and, also, 
because orientation experiments were expected to yield further information 
on the spin assignments and parity changes in the decay of ®4Mn. 


T,21,2,3 Mn 54 
310 DAYS 


% 835 KEV 


1, =0* TI UITHATIT Gr eS 4 


The decay scheme of *4Mn. 


The decay scheme for *4Mn has been determined by Deutsch and 
Elliott (1944) and is shown in fig. 1. The energy available for the 
K electron capture 0-54-+-0-02 mev has been deduced from the @ value of 
—2-162+0-005 mev for the ®4Cr(p, n)54Mn reaction measured by Loving- 
ton, McCue and Preston (quoted by Collins, Nier and Johnson 1952) 


* Communicated by the Authors. 
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From the half-life of 310 days and using the nomogram of Moszkowski 
(1951) a value for log ft of 6-2 is obtained. Comparison of this figure 
with those tabulated by Nordheim (1951) suggests that the transition is 
allowed or first forbidden, involving a spin change AJ of 0 or 1. From 
the existence of coincidences between K x-rays and y-rays observed 
by Deutsch and Elliott (1944) the lifetime of the 0-835 Mev state of °4Cr 
can be assumed to be 10° sec or less. Since the empirical lifetime 
given by Goldhaber and Sunyar (1951) for octupole transitions is of the 
order 10-3 sec and for quadrupole transitions of the order 10-8 sec or less 
it can be concluded that the transition is quadrupole or dipole. 

Since °4Cr is an even-even nucleus J, must be 0 with even parity : 
I, must, therefore, be 1 or 2 and this restricts J; to 3, 2 or 1. Parity 
assignments cannot be made to either of these levels. 

The three Mn** salts that had been investigated by paramagnetic 
resonance are manganese fluosilicate, manganese ammonium sulphate 
and manganese bismuth nitrate (Bleaney and Ingram 1951, Trenam 
1953). The resonance spectrum could be interpreted in all cases by a spin 
Hamiltonian 

A —gfH .S+DS2+AS.1 
with S=5/2 and J=5/2 (for the stable **Mn). Here DS,? represents the 
effect of the crystalline electric field which splits the 6-fold degenerate 
electronic ground state into 3 doublets while AS . I describes the magnetic 


_ hyperfine structure coupling. It so happens that in all three crystals the 


overall Stark splitting is of the order of 0-3°x. This means that below 
about 0-05°K, even in the absence of an external magnetic field, only the 
lowest electronic level will be appreciably populated. This lowest state is 
S,=+5/2 for the fluosilicate and the double nitrate (D<0) and pre- 
dominantly S,=-+1/2 for the Tutton salt (D>0). The magnetic hyper- 
fine coupling is isotropic and since it is of the same order as the Stark 
splitting, alignment of the electron spins is accompanied by alignment of 
the nuclear spins. 

Of these three possible salts we chose the double nitrate because it was 
known from previous work (Daniels and Robinson 1953), that, by incor- 
porating cerium ions as cooling agent, temperatures of the order of 
0-003°K could easily be reached by demagnetization from comparatively 
small fields (10 kilo gauss) at 1°. Our specimen consisted of a single 
crystal of Ce,(NOs), . Mg3(NO3), . 24H,O, with some of the Mg replaced by 
20 micro Curie of >4Mn which had been prepared through the **Fe(d «)°4Mn 
reaction. After demagnetization the anisotropy « of the y-radiation 
was observed as the sample warmed up (e=1—ZJ(0)/I(7/2) where I(@) 
is the counting rate at an angle 0 to the crystal axis). 

Figure 2 shows the observed variation of « with 1/7’*. The temper- 
ature 7 calculated from Curie’s law is known to be identical with the 
absolute temperature down to 0-006°K and deviates little down to 0-003°K. 
The results differ both qualitatively and quantitatively from what one 
would expect. The anisotropy, instead of reaching ‘a temperature 
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independent saturation value at low temperatures, passes pian a 
maximum of 28% at 1/7*=100 and then decreases to 21% at 1/T* =320. 
The dotted line shows the anisotropy calculated theoretically using the 
method of Simon, Rose and Jauch (1951) assuming a nuclear gyro- 
magnetic ratio equal to that of the stable isotope. The calculation has 
been carried to terms in 1/74. At 1/7’=20 the last term contributes 
nearly 25% of the anisotropy. The agreement is poor, and by no choice 
of the gyromagnetic ratio is it possible to fit the experimental and theore- 
tical curves between 1/7’*=1 and 20. In view of the behaviour of the 
anisotropy at even lower temperatures we have not thought it worth 
while to pursue this calculation further. 


Fig. 2 
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Variation of anisotropy ¢ with 1/7* in zero external field. 
The theoretical curve is shown as a dashed line. 


The general behaviour, especially the existence of a maximum in the 
anisotropy, is similar to that found previously in our experiments on ®Co 
in cobalt cerium nitrate (Ambler et al. 1954). It is thought to be caused 
by the internal field due to the cerium ions which in these crystals are very 
close to the divalent ions. 

We decided therefore to carry out experiments using an external 
magnetic field of 1000 gauss parallel to the crystal axis. Such a field does 
not cause a prohibitively large temperature rise since the g value of the 
cerium ion is very anisotropic (g,=0-25, g,=1-84) (Cooke, Duffus and 
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Wolf 1953). If H,, 7, are the starting field and temperature respectively, 
H, the field applied parallel to the crystal axis after demagnetization and 
7’; the corresponding temperature, we have approximately 

9 H/T =9,H,/T;, 
e.g., with H;=30 kgauss, 7';=1°K, H,=1 kgauss, we have 7'=0-0045°K. 
At the same time a field of 1000 gauss is large enough to swamp the effects 
of both the internal field and the crystalline field. 


Fig. 3 


° [OmmeecOn 630. 400s SOMmCOnN 70° (805290. = 100 110°  i20- 130 
/_® 
Variation of anisotropy € with 1/7* in an external field of 1000 gauss. 
A correction for the finite aperture of the counters has been made. 


The general form of the dependence of ¢ on 1/7’ found in these magnetic 
h.f.s. polarization experiments is shown in fig. 3 where it is seen that 
anisotropies up to +90% were recorded ; this corresponds to a polari- 
zation with 95°, of the nuclei in the lowest state. The sign and magnitude 
of this effect show that the y-ray cannot be dipole: therefore it must be 
quadrupole with J,=2. Furthermore, since the anisotropy is so large the 
K capture transition must carry away the minimum angular momentum 
compatible with the spins of the initial and final states : 

ig—1,—I,. 
Thus if [,=3,%,=1 or if [,=2, 42=0. No other transitions can give 
anisotropies greater than 33%. If we assume that the anisotropy 
extrapolates to 100%, we can conclude that °*Mn has either spin 2, 
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decaying by an allowed Fermi transition or spin 3, decaying by an 
allowed Gamow-—Teller transition or a Ist forbidden Fermi transition. 
Although the experimental points of fig. 3 follow qualitatively the 
expected trend we could not obtain quantitative agreement with any 
existing theoretical curve. For this reason it is not possible to distin- 
guish between the transitions mentioned above or to deduce a value for 
the magnetic moment of °4Mn. 


We are greatly indebted to Dr. J. M. Daniels for his generous advice and 
help during the early phases of this work, to Dr. H. Halban for his valuable 
cooperation and to Professor Sir Francis Simon for his interest and 
encouragement. Our grateful thanks are due to the Nuffield Cyclotron 
Group, Birmingham, for the irradiation of the iron target and to 
Dr. P. F. D. Shaw of this Laboratory for his great assistance in effecting 
a chemical separation of the activity. 
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J. PEREZ 
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[Received July 3, 1954] 


ABSTRACT 


The plane polarization of the y-ray emitted in the decay of 54Mn to *4Cr 
has been measured using nuclei oriented by magnetic hyperfine structure 
at very low temperatures. It is found that the excited state of ®4Cr decays 
by an E2 transition, and it is furthermore concluded that this state has 
even parity. 


In an earlier publication (Bishop, Daniels, Goldschmidt, Halban, Kurti 
and Robinson 1952) we have described the technique used in this labor- 
atory to measure the plane polarization of the y-rays emitted by oriented 
radioactive nuclei. The measurement of the electric or magnetic char- 
acter of the radiation gives the relative parity of the nuclear states 
involved. We have recently measured the plane polarization of the 
y-rays emitted in the decay of >4Mn which we oriented by the two methods 
described in the preceding communication (I) (Grace, Johnson, Kurti, 
Lemmer and Robinson 1954). 

The polarimeter was placed to receive y-rays emitted in a direction at 
right angles to the single axis of orientation of the nuclei. For this 
direction a simple relation (Steenberg 1953) exists between the anisotropy 
« defined in the preceding paper and the degree of plane polarization 
p=(1,—1,)/(1,+1,). Here J, is the intensity of o-radiation, 1.e., radiation 
polarized with its electric vector perpendicular to the axis of nuclear 
orientation and J, is the intensity of 7 radiation, i.e., the radiation polar- 
ized with its electric vector parallel to the axis. The measurements 
reported in (I) had established that the radiation was quadrupole and 
for this case p—ae, where a=+1 for electric and —1 for magnetic 
radiation. The ratio of the intensities of the o and 7 radiation is therefore 

I, 1+ae 
a E me Ge 
and the ratio of the coincidence counting rates parallel and perpendicular 
to the axis is 


Ny /N,=(1 +Rn)/(R-+n), 
ate ee oe ee ee 
* Now at the University of British Columbia, Vancouver, B.C. 
+ Communicated by the Authors. 
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where R is the analyser efficiency ratio which may be calculated from the 
differential cross section for Compton Scattering and the geometry of the 
polarimeter. ; 
In fig. 1 is plotted the observed ratio Ny /N,; normalized to the ‘ warm 
rates, against 1/7*, the reciprocal of the magnetic temperature, for 
nuclei aligned by magnetic h.f.s. in the absence of an external field. 
Since « is positive and NV ,/N, is greater than unity we see that a is positive 
and the radiation is electric in character. 


Fig. | 


The ratio of counting rates Ny/N, against 1/7'* for aligned nuclei. 


In order to obtain a larger effect polarization of the ®4Mn nuclei was then 
produced by the application of a homogeneous magnetic field of 600 gauss. 
Owing to the restricted space between the magnetizing coils and the 
necessity of shielding the photo-multiplier tubes from the fringing magnetic 
field, the three scintillation phosphors constituting the polarimeter were 
placed at the end of perspex light guides. Curve A of fig. 2 shows 
N,/N, against 1/7*. Curve B shows the anisotropy « for a magnetic 
field of 600 gauss. This had to be determined separately since the 
measurements mentioned in (I) were done with a field of 1000 gauss. 

These results confirm those obtained with the aligned nuclei and prove 
conclusively that the excited state of ®4Cr decays to the ground state 
(even parity, spin 0) by the emission of electric quadrupole radiation. 
It follows that the excited state has even parity and spin 2. 


We wish to thank Dr. H. Halban for his continued interest and 
encouragement. 


Polarized y-Rays from Oriented ®4Mn Nuclei 1199 


Fig. 2 


50 US 


N,,/N, (curve A) and ¢ (curve B) against 1/7'* for nuclei polarized by magnetic 
h.f.s. in a field of 600 gauss. The points refer to Vy/N,. 
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CXXXIV. Interaction of 19 Mev Deuterons with Carbon 
By R. G. Freemantie,* W. M. Greson t and J. Rorsxat } 
[Received September 3, 1954] 


Tun 19 ev deuteron beam from the Birmingham cyclotron was used in 
a study of the nuclear reactions taking place at the bombardment of 


Fig. 1 
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carbon. The scattering camera described by Burrows, Powell and Rotblat 
(1951) was used, after slight modification. The target was acetylene gas 
at a pressure of 29-3cm; the charged products of the reactions were 
recorded in photographic plates (Ilford C2, 400 microns thick). 

Range experiments were carried out at a number of angles of emission 
of the particles. Histograms plotted for each angle showed the presence 
of a number of groups of particles, most of which could be identified from 
the measurement of the energy and its variation with angle. Apart from 


Fig. 2 
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the group of protons due to the scattering by hydrogen, several other 
proton groups were identified as originating from the C(d, p) reaction 
and corresponding to different energy states of 8C. Three deuteron 
groups were identified : the elastic group, and two inelastic groups corre- 
sponding to the energy levels of #2C at 4:43 and 9-61 mev. No trace was 
found of the level reported at 7-5 or 7-68 Mev, which has been observed in 
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several reactions, i.e. *Be(«, n) (Guier et al. 1952) and ™*N (d, «) (Holloway 
and Moore 1940, Dunbar e¢ al. 1953, Prowse 1954). The absence of this 
level in d, d scattering cannot be due to a selection rule for the isotopic 
spin, as the lowest level of 2C with 71 must lie much higher ; the fact 
that no sign of this level could be found must therefore be due to a different 


cause. 
Fig. 3 
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The angular distributions in the centre of mass system for the three 
groups of deuterons are shown in figs. 1 and 2; the ordinates give the 
differential cross section in millibarns per sterad. The errors indicated 
are statistical only, other relative errors are quite small. The absolute 
cross section may be too low by about 5%, because of a possible systematic 
error in the beam integrating system. Figure 1 contains also the curve for 
Rutherford scattering. No analysis of these distributions has so far been 
carried out. 
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Figures 3 and 4 give the differential cross sections for the protons from 
the C(d, p) reaction corresponding to #C being formed in the ground 
state, and in the 3-08 and 3-89 Mev states respectively. The latter group 
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contains some protons from the 3-68 Mev state which was not resolved. 
These angular distributions are qualitatively similar to those observed 
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with an 8 Mev deuteron beam (Rotblat 1951). An analysis of these curves, 
according to Butler’s stripping theory and assuming a value 7»=4-5 x 10~-¥# 
em, was carried out by Mr. W. M. Fairbairn. The dashed curves in figs. 3 
and 4 are the normalized curves calculated for J,=1, 0 and 2 for the three 
states respectively. As is seen, the agreement is very good for the excited 
states, confirming the assignment of even parities and spin values of 1/2 
to the 3-08 Mev state and 3/2 or 5/2 to the 3-89 mev state. In the case of 
the ground state, however, the agreement is very poor, and experimental 
points at smaller angles would be necessary to confirm from this experi- 
ment the assignment of J, =1 for the transition to this state. 
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CXXXV. Angular Distribution of Neutrons from the Reaction “B(d, n)??2.C 
at 0-85 Mev Bombarding Energy 


By ARNFINN GRAUE 
Fysisk Institutt, Norway * 


[Revised MS. received August 19, 1954] 


SUMMARY 


Angular distributions and energy spectra of neutrons from ™B(d, n)?C 
have been measured at 0:85 Mev bombarding energy, using the photo- 
graphic plate technique. Energy levels in 1°C are found at 4-4, 7-7, 9-6 
and 12-7 mey. The stripping process is suggested to be responsible for a 
large part of the neutron yield. The orbital angular momentum of the 
captured proton has been determined as 1,—1 and 2 for the 4-4 Mev and 
9-6 Mev level respectively. 


§ 1. INTRODUCTION 


SEVERAL experiments have confirmed the existence of levels in °C at 
4-4, 7-5 and 9-6mey. The 14B(d, n)¥C reaction has been studied by 
Johnson (1952) and Gibson (1949). Johnson failed to observe the 
7-5 mev level whereas Gibson from his data concluded that this level could 
be confirmed less certainly than the 4-4 and 9-6 Mev levels. 

Spin and parity assignments have been made only for the ground state 
and the 4:4 Mev level. Gibson (1953) has measured the angular distri- 
butions of neutrons leaving 12C in its ground state and in its first excited 
state at a deuteron energy E,—8-1 Mev. In both cases he found agree- 
ment with Butler’s stripping theory (Butler 1951) for /,=1. 

Investigations regarding the competition between stripping and 
compound nucleus formation can be done at smaller bombarding energies. 
Burke et al. (1954) have measured the angular distribution of neutrons 
emitted in the 14B(d, n)!2C reaction at bombarding energies 0-71, 1-06 
and 1:43 Mev. The distribution curves had peaks in the forward direction. 
As, however, neutron groups of all energies were counted, further experi- 
ments with energy discrimination are desired. For the present measure- 
ments of the angular distributions of neutrons from 11B(d, n)!#C the 
photographic plate technique has been used. At the moment this method 
seems to be the best one available for investigations of the more 
complicated neutron spectra. 


§ 2. EXPERIMENTAL CONDITIONS 


The deuterons were accelerated in the Van de Graaff generator at the 
University of Bergen. The target used was a separated ™B target 
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obtained from Dr. J. Koch, Institute of Theoretical Physics, Copenhagen. 
The target-thickness was approximately 25 kev for 0-85 Mev deuterons. 
Ilford C2 Nuclear Research Plates with emulsion thickness 200 were 
used. The plates received the neutrons at grazing incidence. They 
were arranged in such a way that they covered the range of the angle 
6 from 0°-160°, 6 being the angle between deuteron and neutron directions. 
The distance between target and the scanned areas was 125-135 mm, the 
same for all angles 6. The total deuteron charge on the target was 
100 wA-hours, which gave a suitable density of recoil-proton tracks in the 
emulsions. The methods for measurement and calculation of the neutron 
energy and the angular distribution have been described by many workers 
in previous papers (e.g. Trumpy ef al. 1952 and 1953). 


Fig. 1 
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0° neutron spectrum from the !B(d, n)!2C reaction. 


§ 3. Toe Neutron Enrerey Spectra 


Neutron energy spectra have been measured at the following angles @ : 
0°, 15°, 30°, 45°, 75°, 105°, 120°, 135° and 155°. A total of about 6000 
tracks have been measured. Figure 1 shows one of these spectra. This 
spectrum is corrected for the variation with energy of the neutron—proton 
scattering cross section and for the escape of protons from the emulsion. 
The errors indicated are statistical errors only. Because of the lower 
intensity of the higher energy neutron groups, a greater area was scanned 
for measurements of these neutrons. This is the reason for the dis- 
continuity in the statistical errors. 

All the peaks shown in the spectrum correspond to levels in !C except 
those marked otherwise. In the lower part of the observed spectrum 
homogenous groups from 1B(d, n)C reaction are superimposed upon the 
continuous spectrum from 14B(d, n)3« reaction. 

All the neutron spectra have a high energy group corresponding to 120 
in its ground state. In accordance with previous measurements they also 
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demonstrate the existence of energy levels in ?C at 4-43, 9-67 and 
12:67 Mev. The neutron spectra also indicate a level in C at 7-67 Mev. 
This level was found at all angles. The existence of this neutron group 
has been less certain in previous measurements (Gibson 1949) on the 
11B(d, n)!2C reaction. It is, however, clearly resolved in fig. 1. The 
group can not be due to any contamination of the target because no 
conceivable contaminant gives rise to a neutron group of this energy. 
The only contamination groups observed in the spectra are from the 
D(d, n)*He and the *N(d, n)45O reactions. Both D and 'N have 
previously been observed to be picked up from the vacuum system during 
exposure (Hnge 1952). 


§ 4. ANGULAR DISTRIBUTION OF THE NEUTRONS 


No theory for the angular distribution of the neutrons at low deuteron 
bombarding energies has so far been published. For energies higher than 
the Coulomb-barrier there are theories based on the stripping process. 
These theories, which predict a maximum in a forward direction, have also 
been successfully applied for deuteron energies lower than the Coulomb- 
barrier. However, in this region it is expected that compound nucleus 
formation also contributes appreciably to the reaction yield. 
Wolfenstein’s theory (Wolfenstein 1951) for compound nucleus formation 
predicts a distribution symmetrical about 90° in the centre of mass system, 
and generally with maxima in the forward and backward directions. 

Angular distributions of the neutron groups corresponding,.to the ground 
state, 4:4 Mev and 9-6 Mev levels are shown in fig. 2. Considering first 
the reaction forming 12C in its ground state (fig. 2 (@)), we find a maximum 
in the forward direction. Such a distribution is characteristic for the 
stripping process, and theoretical curves were calculated, using the 
stripping theory in the form given by Bhatia et al. (1952). Two curves are 
plotted corresponding to [-values of the captured proton of 1 and 2. The 
nuclear radius used in the calculations was k=6-4 x 10-18 cm computed 
from the expression R=(3-7+ 1-22 x A183) x 10-18 em used by Bhatia et al. 
There is no good agreement between the experimental distribution and any 
of the theoretical curves. The /-value of the captured proton, however, is 
known from previous measurements at 8-1 Mev deuteron energy to be 
1,=1 (Gibson 1953). Butler has reported some preliminary results 
about the influences of Coulomb forces on the deuteron stripping cross 
section (Butler 1954). He then finds that the maximum in the angular 
distribution curve is somewhat displaced towards larger angles. ‘This 
would may be give a better agreement between the experimental results 
and a theoretical /,—1 curve. 

The reaction forming 1°C in its 4-4 Mev excited state also shows a distri- 
bution with a maximum in the forward direction (fig. 2 (b)). The curve 
with /,,=1 gives the best fit which is in accordance with Gibson’s results at 
E,=81mev. Since the ground state of 1B has spin 3/2, and odd parity, 
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the 4-4 Mev excited state of 12C must have even parity and spin 0, 1, 2 or 3. 
The spin of this state is known to be 2 (Ajzenberg and Lauritsen 1952). 

No parity assignment of the 9-6 Mev excited level of 12C has so far been 
given. The distribution in fig. 2 (c) shows a minimum at 0° and a maxi- 
mum at 80°. If the stripping process is responsible for this distribution 
the calculated curves show that protons with /,—2 are captured. This 
leads to the results that the 9-6 Mev excited level of °C has odd parity 
and spin 0, 1, 2, 3 or 4. 

Fig. 2 
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Angular distributions of the neutrons from the UB(d,n)#?C reaction and 
theoretical curves calculated from the stripping theory in the form given 
by Bhatia et al. (1952) for angular momentum transfers of 1 and 2. 
(a) °C formed in its ground state, (b) 12C formed in the 4:4 Mev excited 
state, (c) °C formed in the 9-6 Mev excited state. 


Because of the low intensity of the neutron group corresponding to the 
7 cu Mev level of C, the statistical errors in the experimental angular 
distribution for these neutrons are large. It can be stated, however, that 


the deviation from spherical symmetry in the centre of mass system is not 
more than 25%. 
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The agreement between the experimental curves and those calculated 
from the stripping theory is quite good for the 4:4 and 9-6 mev levels. 
This is surprising because the theory takes no account of the Coulomb 
effects, and since contributions from the compound nucleus formation are 
expected. At large angles one should not expect the experimental and 
theoretical curves to agree. 

The problem regarding the relative importance of the stripping process 
and the compound nucleus formation can only be solved on the basis 
of further experimental investigations and a theory for the stripping 
process valid also for low deuteron energies. 
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CXXXVI. CORRESPONDENCE 


On the Acceleration of Cosmic Ray Particles 


By W. B. THomMPson 
Atomic Energy Research Establishment, Harwell* 


[Received July 21, 1954] 


A MECHANISM which seems capable of accelerating particles to cosmic ray 
energies has been proposed by Fermi (1949), and has recently been 
extended by Fermi, and by Morrison, Olbert and Rossi. Energetic 
charged particles are elastically scattered by magnetized clouds of 
partially ionized gas in interstellar space, and since head-on collisions 
are slightly more frequent than those in which the particle overtakes the 
cloud, a particle experiences a gradual increase in energy at a rate 


1 dH v\? 
edt Cy 


where v is the mean velocity of the gas clouds, ¢ the velocity of light, 
E the energy of the charged particle, assumed relativistic, and r the rate 
at which the particle is scattered. 

If particles have life-times which are distributed exponentially about 
a mean, this theory predicts the correct power law for the cosmic ray 
energy spectrum, the exponent in this fixing a relation between the 
mean life and the rate of energy increase. The mean life is determined 
by the presence of heavy nuclei in the primary cosmic rays, which implies 
that the time available for acceleration cannot be much greater than the 
average time between nuclear collisions involving cosmic ray particles, 
which is of order 10’ years. Since the velocity of magnetized clouds v 
can be associated with the observed turbulent velocities of interstellar 
matter, about 30 km/sec, the theory determines a relation between the 
maximum energy for which the power law spectrum holds and the 
frequency of collisions between particles and clouds. Morrison, Olbert 
and Rossi conclude that the extension of the power law spectrum to 
101° ev implies scattering centres about a light year apart, an uncomfort- 
ably small separation. Recent air shower experiments (e.g. Cranshaw 
and Galbraith) have shown that the power law spectrum holds for 
energies of 10! ev ; hence some more efficient accelerating mechanism 
must be found. 

Fortunately such a mechanism is concealed in the Fermi theory. In 
the irregular magnetic field one expects to find fluctuations varying in 
diameter from several light years to small fractions of a light year ; and, 
while the small scale fluctuations by themselves do not seem sufficient 
to accelerate particles, they enable the large scale clouds to operate much 
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more effectively. The scattering by these small irregularities changes the 
paths of the energetic particles from regular spirals to irregular diffusive 
tracks, and a particle undergoing such a motion can gain energy from a 
large scale fluctuation without being reflected. If, as a model of a large 
disturbance, we consider a region in which the field changes rather 
rapidly by an amount dH, remains constant for some space, then returns 
to its original value, one can approximate the solutions to the equations 
of motion by expanding in powers of 5H/H. The first term in this. 
expansion is enough to describe the reflections proposed by Fermi ; 
while the second has associated with it an energy increase at a rate 
1 dH dH\2 
Edi (=r) 

where the constant K depends on the details of the model, and for a 
uniform distribution of particle velocities, and linearly varying fields 
is 9/4. Since 5H/H, while assumed small, can easily be 10%, this 
acceleration rate may easily be 100 times that given by Fermi. This 
mechanism gives importance to relatively small variations in field 
strength, and more recent pictures of the galactic magnetic field 
(Chandrasekhar and Fermi 1953 a, b) suggest that the field is somewhat 
more uniform than had previously been supposed. If the galactic field 
has the regularity ascribed to it, and if the acceleration process described 
here is operative, it is possible to extend the theoretical range of the 
power law spectrum to about 10!% ev, but it is difficult to see how it 
can hold for much greater energies. 


Vv 
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~ Kr 
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On the Evaluation of Characteristic Temperatures for Cubic Crystals 


By A. B. Batra and G. E. TAUBER 


Division of Physics, National Research Council, Ottawa* 
[Received June 21, 1954] 


Tur calculation of the Debye characteristic temperature @» involves the 
evaluation of the integral 


T= | (2006, Or] eee oe en) 
Sayan ee ee ee See ere 
* Communicated by the National Research Council, Ottawa. 
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where v, (0, 6) (i=1, 2, 3) are the velocities of the three elastic waves for a 
given direction (0, ¢) of propagation and integration is over all directions ot 
In terms of elastic constants the v; are given by the well- 


propagation. 
known Christoffel’s equations of elasticity theory. Op and J are related 
by the equation h /9s\ 1/8 

0,= (3) igi Sk 


in which J is the cell-volume, s the number of atoms per cell, and other 
symbols have their usual meaning. 

For elastically anisotropic crystals [j= 2c44(c1;—¢12) *~1] J may either 
be evaluated numerically or by the Hopf—Lechner (1914) method ; both 
involve lengthy calculations. Recently Blackman (1951) has obtained a 
semi-theoretical formula for J valid for 7 sufficiently greater than unity. 
One of the constants occurring in this formula was fixed by fitting the 
value of @p (obtained numerically) for Li. 

The purpose of this note is to obtain an approximate formula for I 
in the following manner : since the integrand in (1) is invariant under the 
operations of the cubic symmetry group, we may expand it in terms of the 
appropriate Kubic Harmonics K,,, introduced by Von der Lage and Bethe 
(1947).t Thus ite 
jp | ['SaeK |dQ—4ree | en 

m=0 


The first three lowest order Kubic Harmonics, apart from normalization 
factors, are (7?+-y?+-2?=1) 


Ky=1, K,=(«*+y*+24—3), 


and K,=xy'2+ 2 K,— 4. . 
Now if one assumes that the first three terms in the integrand of (3) give a 
reasonable approximation for (1), one may determine «, «%, and a, 
analytically since for three directions in a cubic crystal, viz. (100), (110) 
and (111), the v; can be easily obtained in terms of the elastic constants. 
This procedure gives for J the expression 


T= 3 [10e,1-*?-+16(c,,—}e)-9-+ 16(Cqy-+ $e) 9? 
+9(¢y4—3e)-9/2-+ 18(6 44+ Fe)-9/24 36c4,-9/2], . . (5) 


with e=c,;—Cy9—2¢ 44. 

With the help of (5) and (2) we have evaluated @, for several crystals. 
The results are given in table 1. @p values calculated by other methods 
are also included. Since according to Blackman (1951) the latter estimates 
may be in error by several per cent, the @-values calculated here seem 
reasonable. 


+ It ought to be mentioned that the present method has been previously 
employed by W. V. Houston (1948, Rev. Mod. Phys., 20, 162), to obtain an 
approximate form for the frequency spectrum of crystals. 
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It should be noticed that (5) would be nearer to the correct value of (1), 
the closer 7 is to unity. A glance at the table shows that the formulae 


Li Na K Cu Ag |) Au Al KCl | NaCl | 


@,(P)* _| 324 LS es 329 | 220 |165 | 405 | 239-5 | 311 
©7p(O) 354 144 cae 342 4) 212-. |. 158 ~— | 394-227 305 
©,(B) 351 147 82 | 355 | 203 | 160 | 400 oa — 


2c : 
= 7-8 8-2} 8-4 3:2 2-9 2-9 1-2 0-52 0-74 


* @>(P)=present calculation; ©p(B)=Blackman’s values (loc. cit.) ; 
©p(O)=other calculations. 


given here are approximately valid over a wide range of values of 7. The 
accuracy in any particular case can of course only be estimated by calcu- 
lating the influence of the next higher harmonic on a». 


The authors are grateful to Mr. D. Patterson for many helpful 
discussions and to Miss C. de Chantigny for help with the numerical work. 
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The Elastic Scattering of Protons by Fluorine 
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[Received August 20, 1954] 


Tur differential elastic scattering cross section of fluorine for protons 
has been measured at four angles, with bombarding energies from 0-5 to 
9-1 Mev. Anomalies were observed corresponding to the resonances 
in 19F--p at 669, 831, 874, 935, 1355, 1381, 1431, 1690, 1940, and 2030 kev. 

A beam of protons, monoenergetic to 1 kev and defined to 2mm by 
3mm was provided by the Cavendish electrostatic generator. The 
target was 3 kev thick lithium fluoride evaporated on to a film of carbon 
a few hundred angstroms thick, made by the technique due to Bradley 
(1954). The scattered protons were analysed by a rotatable magnet 
analyser, and detected with a proportional counter ; protons scattered 
from 19F, 2C, 7Li and °Li in the target were well resolved. The scattermg 
angles chosen were those for which the respective Legendre polynomials 
of the angle in the centre-of-mass system are zero, i.e. 90°, 125-3°, 140°8°, 
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159-8°, since an anomaly due to resonance scattering of protons of angular 
momentum / disappears at an angle for which P",=0. A Geiger counter 
detected gamma-rays from the target, and provided a voltage calibration 
of the generator. Monitoring was by integration of the beam current, 
after D.C. amplification. 


Fig. 1 
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The 1°F(p, p) elastic cross section, at 140-8° in O.M. system, and gamma-ray 
excitation function. 


Figure 1 shows the elastic scattering cross section curve obtained at 
141°, together with the gamma-ray excitation function. The relative 
accuracy of the points on the curve is around 2%; the absolute cross 


section, obtained by comparison with 7Li (Wart F 
accurate to 10%. (Warters et al. 1953) is probably 
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Detailed analysis of the results is in progress, but certain features of 
the levels in ?°Ne are apparent. Thus, in the resonances for the reaction 
19K (p, xy)?°Ne, the alpha-particle yield (Chao 1950) gives a value for 
the ratio of proton to alpha partial widths, but does not indicate which 
is the larger. The marked scattering anomaly at 669 kev has been 
fitted, taking a value of proton partial width, I’,, equal to the total 
resonance width I’ which is 7-5 kev. The gamma-ray width is negligible. 
The alpha particle yield shows that [.,=0-98F or 0-021, and so the 
present evidence leads to the conclusion that the proton width is larger, 
and some fifty times the alpha width. However, at the 935 kev resonance, 
which is also known to be (1, +) (Seed and French 1952), the scattering 
curve cannot be fitted by [,=I, the observed anomaly being quite 
small. For this resonance, [.,=0-9I or 0-1’, and it therefore appears 
that the proton width takes the smaller value. The 669 kev level also 
differs from the other (1, +) resonances, at 340, 935 kev, in that some 
2°, of the transitions from it occur by gamma emission to the first 
excited level of 2°Ne (Devons and Hereward 1948, Rae et al. 1950, 
Carver and Wilkinson 1951). 

Sinclair (1954) has observed three other resonances for the reaction 
19F(p, y)2°Ne at 1092, 1324, and 1431 kev. The first corresponds to 
a resonance in 1°F(p, «y)!®O, whereas the last two do not. No appreciable 
variations in the proton scattering cross section occur at the first two 
of these levels, but at the third there is a large isotropic anomaly. The 
1431 kev level is therefore formed by s-wave protons, allowing assignment 
(0, +) or (1, +), but not (0, —) as suggested by Sinclair. No alpha 
emission to levels in 1*O with even spin and parity are observed from 
this level (Chao 1950), so that it is probably (1, +). 

For the levels formed at 874, 1355 and 1381 kev, the anomalies disappear 
at 90° and become more marked as the scattering angle increases, showing 
formation by p-wave protons. This agrees with the assignments, 
(2, —) for the 874 kev level (Seed and French 1952) and 1381 kev level 
(Sanders 1953). The spins and parities of the 1355 and 1381 kev levels 
deduced tentatively by Chao (1950) from the gamma-ray distribution 
and quoted by Ajzenberg and Lauritsen (1952) are incorrect. As (p, ay) 
transitions to levels in 10 with odd spin and even parity occur from the 
1355 kev level, the assignment must be (2, —). The small size of the 
874 kev scattering anomaly in relation to the resonance cross section 
and width supports the argument of Wilkinson (1935) based on the 
radiative width, that the proton width is small. Similarly, at 598 kev, 
no appreciable anomaly occurs, though the resonance cross section is 
high, again indicating that I’, takes the smaller of the two possible values. 

It remains unexplained why, of the four (1, +) levels at 340, 699, 
935 and 1431 kev, two decay by alpha-emission, one by gamma-emission 
alone, and one by both processes. 

I am indebted to Dr. A. P. French for suggesting this problem, to 
Mr. E. 8. Shire for permitting the use of the electrostatic generator, 
and to the Department of Scientific and Industrial Research for a grant, 
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On the Interaction of a Dislocation and a Charged Solute Ion 


By Routr LanpAvUER 


Research Department, International Business Machines Corporation, 
Poughkeepsie, New York 


[Received August 30, 1954] 


CoTTRELL, HUNTER AND NABARRO (1953) have pointed out that the 
electric field surrounding a dislocation in copper may exert an electrical 
force on a polyvalent substitutional solute atom. According to Friedel 
(1952) it is not likely that the polyvalent solute atom contributes more 
than two electrons to the conduction band of the copper matrix. The 
electrons which are not contributed to the conduction band, have an 
energy below the bottom of the conduction band, and are bound to the 
solute atom. This still leaves the possibility of having, for example, 
a Zn** ion in solution, screened only by conduction electrons. Now 
if we consider motion of the Zn++ ion, the screening cloud must follow 
the ion, and the motion does not correspond to any net charge transfer. 
Therefore, there are no purely electrostatic energy changes, and con- 
sequently no net electrostatic forces. Cottrell, Hunter and Nabarro 
have only taken that part of the screening cloud which lies within the 
cell of the Zn** ion as being effectively tied to the Zn++ ion in its motion, 
but it is not clear why this arbitrary cut-off can be introduced. 

It should perhaps be pointed out that there is a great deal of evidence 
(Schwarz 1940) that an electric field in a conductor can move solute ions, 
despite the presence of screening electrons. However, this evidence 
deals with non-equilibrium situations in which there is an electronic 
current. One cannot talk in any simple way about the energy change 
of the screening cloud, as it is moved about, in the presence of the current. 
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Even if one discounts the direct electrostatic force on the Zn+* ion, 
it may still seem possible that a charged solute atom interacts with 
dilations, in a manner which is not satisfactorily represented by the 
usual picture (Cottrell 1953) of a misfitting atom, which produces elastic 
strains. After all, one does not convert a copper atom to a zinc atom, 
be inserting an extra layer of copper, as the usual elastic model presumes, 
but rather by adding an extra nuclear charge to a copper ion and sim- 
ultaneously an extra electron to the conduction band. If this is to be 
represented by a classical elastic model, it seems more reasonable to 
consider the extra charges that have been added to the original copper 
matrix as providing localized body forces acting on the pure copper that 
was initially present. 

We shall assume that these volume forces f(r), have radial symmetry 
about the solute ion, that the forces are directed radially, and exist in an 
isotropic medium with bulk modulus K. Then the change in volume of 
a finite sample, taking into account Eshelby’s (1954) image terms can 
be shown to be 


5V= ag |flr)- ra, Wer sa) 


where the range of integration can be limited to the small volume in 
which f(r) is appreciable. If we subject the copper matrix containing 
one set of these volume forces, to an externally applied dilation, jabs 
then the work we do against the volume forces is 


U=— SE ft(r) dV =—4,K (8). a sae tdci (2) 


In the more usual argument (Cottrell 1953) one pictures an inclusion 
which is too large for its cavity of radius r, by a factor (1+), where 


1 da 
== eT ° e e e ° ° e 3 
<~ a de (3) 
In (3), a is the lattice constant of a dilute solution and ¢ is the atomic 
concentration of solute atoms. The energy of interaction of this over- 
sized atom with a hydrostatic pressure is then taken to be either 


U;=—4rer®K A, ° e ° ° e e (4) 


or else 3(1—v)/(1-+-v) times that value where v is Poisson’s ratio. The 
latter factor serves to include the strain energy in the solute atom. (3) 
combined with (4) yields 


cen) (Ha) (Ha) 0 


where N is the total number of atoms in the specimen under consideration. 
The first factor in (5) is the volume under consideration. The second. 
factor is the fractional change in volume that results from the addition 
of one solute atom. We therefore have U;=—KA, (5V) in agreement 
with eqn. (2) which we found on the basis of the distributed force model. 


1218 Correspondence 


The equation, U;=—KA, (8V), simply expresses the fact that if the 
energy of solution depends on the dilation of the lattice, then this energy 
dependence provides a force tending to change the crystal volume away 
from the value it has in the pure solvent. Consider a volume of copper 
into which the extra electron and nuclear charge are introduced to con- 
vert a Cu atom to Zn. Let this energy of conversion or solution be 
E,, where H, assumes that all copper nuclei are held in place. If we 
then permit the nuclear positions to relax, we can describe the relaxation 
by a set of coordinates £;, in addition to a dilation 4,. The coordinates, 
B;, are chosen so that in the pure copper matrix a dilation, 4, does not 
simultaneously cause a change in the value of any of the 8;. Furthermore 
the 8; are taken to be zero in pure unstrained copper. The energy of 
a unit volume, containing n solute atoms can then be written 
B=By +n +n St A, + Sn SB ts KAP+ DaghB. . . (6) 
i OB," 2 ii 
Here EH, is the energy of the initial state, before the conversion. The 
last term gives the elastic energy associated with the B;. If we permit 
relaxation we must have 


04 =n a7 +Kk4,=0, Pele 5 3 ee (7) 
or ok, Ka, 
aA = = KO). o& Wil eisai 


Here, as before, 5V is the volume change produced by a single conversion 
of copper to zinc. The interaction with an externally enforced dilation, 
os 


OE, n: 
U1= 4g 5 = — 4aK (8V) See oe 
in agreement with eqn. (2). 
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Further Evidence for the B-Decay of K-Mesons 


By C. DananayakeE,* P. E. Francois, Y. Fustmoro,;, P. [REDALE, 
C. J. WADDINGTON and M. YAsIN 
H. H. Wills Physical Laboratory, University of Bristol 


[Received September 15, 1954] 


AN event has recently been reported from this laboratory which appears 
to be due to the decay of a K-meson into an electron of ~ 90 Mev 
(Friedlander et al. 1954). In this letter we report a similar event in which 
the secondary electron has an energy of only 50 Mey. ‘This event was 
found in the course of an investigation, currently proceeding in this 
laboratory, in which K-mesons emitted from stars are detected in stacks 
of stripped emulsions (Dahanayake ¢¢ al. 1954, see also Report of Padua 
Conference 1954). 

The primary particle, K,,, was emitted from a star of type 19+ 6p, 
and traversed five plates with a total range of 11-79 mm before coming 
to rest. All other tracks associated with the parent star, with grain- 
densities above plateau, have been examined and appear to be due to 
stable particles. The mass of the K-meson, measured by ionization- 
range is 970+100m,; by scattering-range (Menon and Rochat 1951, 
Menon and O’Ceallaigh 1954) 690-+160m,; and by scattering-range 
using the constant sagitta method (Dilworth et al. 1954, Fay et al. 1954), 
760+110 m,. , 

The secondary particle, which has an ionization indistinguishable 
from the plateau value, traverses 23 plates with a total range of 25:6 mm 
before leaving the stack. Initially the length per plate of the secondary 
track is about 2 mm, but it becomes steeper, and just before leaving the 
stack the length per plate is only 0-8 mm. 

For convenience the track was divided into 5 segments and blob 
density-scattering measurements were made on each. These results 
are shown in the table. 


Scattering and Blob Density Measurement 


Multiple scattering pp 
Segment Length/plate Normalized angle per 100 u in Mev/c 
mn blob density oc) a (3) from a) 


1 6-35/4 plates 1-0040-04 0-5240:04 0-4640-04 4914 
2 5-36/4 plates 1-0340-08 0-6440-05 0-7140-06 31438 
3 5-68/5 plates 1-06-£0-09 0-82-4006 0-7040-05 3142 
4 3-96/5 plates rod 0-82-40-08 0:7040-:07 3143 
5 3-76/5 plates = ~1 0-82-40-:09 1:2640:15 1742 


x2) ig obtained from the second differences of coordinates, and a3) from the 
third differences to eliminate distortion. The value of pf is obtained from a). 
ee 


* On leave of absence from the University of Ceylon. 
+ On leave of absence from the University of Kyoto. 
+ Communicated by Professor C. F. Powell, F.R.S. 
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Values of Z‘2) were obtained by noise-elimination, using the procedure 
of Daniel and Perkins. Values of ~‘) obtained by the method of third 
differences are also quoted, as an indication of the comparative lack of 
distortion in the plates. A rigorous 4D-‘ cut off’ was employed. The 
errors quoted are purely statistical. 

Extensive measurements in the present stack of plates (Fowler and 
Perkins 1954) indicate that a j.-meson with an ionization near the plateau 
value must have a value of pf greater than about 120 Mev/c. It therefore 
appears from the results obtained in the first segment alone, which 
provides the most favourable conditions for measurements, that the 
secondary particle is most probably an electron. This identification 
is supported by the results obtained from the other segments. 

Since the observed track-length of 2-56cm is ~ 0-9 cascade units, 
we may expect an appreciable decrease in the value of p8 due to brems- 
strahlung if the track is indeed due to an electron. The values of pf 
obtained in the five segments of track are plotted in fig. 1 and may be 


Fig. 1 
ENERGY of SECONDARY os o FUNCTION of DISTANCE TRAVERSED 
60 


sO 


40 


w 
ce) 


Energy in Mev. 


Distance Traversed in Cascade Units. 


Plot of the energy of the secondary as a function of the distance traversed in 
cascade units. (a) shows the expected variation due to ionization loss 
alone, and (b) that due to ionization loss together with average 
bremsstrahlung losses. Points marked x show values calculated from 
2nd differences, and those marked © values from 3rd differences. 


compared with the expected variations due to (a) ionization loss alone, 
and (b) to ionization loss together with average bremsstrahlung losses. 


* appears that the particle has in fact suffered some bremsstrahlung 
Oss. 
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The present event increases to four the number of heavy mesons which 
apparently undergo B-decay (Friedlander et al. 1954, O’Ceallaigh 1954). 
The spread of the pf values in these events suggests that the B-decay of 
the K-meson is probably a three-body process. Since a secondary 
particle can be identified as an electron only if’the initial value of pB 
is very low, or if there is considerable bremsstrahlung loss, the results 
suggest that K-mesons not infrequently decay in this mode. 
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Paramagnetic Resonance at Radio Frequencies 
in Low Temperature Carbons 


By D. J. E. Ineram and J. G. TAPLEY 
Department of Electronics, University of Southampton 


[Received September 21, 1954] 


Tx order to confirm (Ingram and Bennett 1954, Castle 1954 a, Bennett 
et al. 1954) that two distinct types of paramagnetic resonance absorption 
occur in different kinds of carbon we have now made measurements in 
the radio frequency region of 25-50 Me/s. The apparatus used is a 
modification of Pound and Knight’s nuclear spectrometer (1950), the 
carbon ‘samples filling the inductive coil, which has a volume of about 
0-5 cm, The magnetic field is obtained from a Helmholtz coil-pair 
fed with current at 50 c/s, no steady d.c. field being produced. Since 
the magnetic field strength required for resonance at these frequencies 
is only of the order of 15 gauss, the modulating field, of some +20 gauss 
amplitude, will go through resonance at each side of zero, and hence 
two absorption lines are obtained on the oscilloscope screen, as seen In 
the figure. The frequency of the absorbed radiation (v) is measured by 
a crystal-controlled wavemeter. This determination also acts as a 
direct calibration of the magnetic field, since the spectroscopic splitting 
factor of the carbon samples is the same as that of the hydrazyl marker 
to within -—-0-001, and hence the separation between the centres of the 
+wo resonant lines is given by 2hv/gB (=7-178 x 10-*v gauss). 
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The photographs shown in the figure are the resonances obtained 
from (i) diphenyl trinitrophenyl hydrazyl, (ii) eek activated, 
decolourizing charcoal, (iii) coal carbonized at 550°c, this being the 
carbonizing temperature which produces the highest free-radical con- 
centration (Ingram et al. 1954). The different widths of all three samples 
can be clearly seen, and the measured values for full width at half-power 
are (i) 2-5-+-0-5 gauss, (ii) 4-5+2-0 gauss and (iii) 842-0 gauss. The fact 
that the widths of the resonances from the two different carbon samples 
are still the same as those observed at microwave frequencies (Bennett 
et al. 1954) is in striking contrast to the results obtained by Henning et al. 
(1954) on graphites; where the width changed from the 8 to 40 gauss, 


Paramagnetic Resonance Absorptions at 40 Mc/s. 


Top: Diphenyl trinitrophenyl hydrazyl 
Centre: Activated, decolourizing charcoal 
Bottom : Coal, carbonized at 550°c 


observed by Castle (1953, 1954 b) at 9000 Me/s, to 0-1 gauss at radio 
frequencies. Any increase in width due to random size, or orientation, 
of the particles is ruled out by the fact that both carbon samples have 
the same range of particle size, and previous measurements on graded 
sizes have shown no change in width (Bennett et al. 1954). 

These low frequency measurements therefore confirm the suggestion 
(Bennett e¢ al. 1954) that there are two different mechanisms responsible 
for paramagnetic absorption in carbon. (i) In carbons which have not 
been heated above 650°c there is a resonance due to free radicals, or 
‘broken bonds ’, in the carbon structure. this normally having a width 
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of 8 gauss, independent of frequency of absorption ; and (ii) in graphites, 
formed above 1200°c, there is a resonance from the conduction electrons 
in the two-dimensional plane, the width of this absorption varying 
markedly with frequency. This difference is also supported by the 
recent measurements of signal strength against carbonizing temperature 
(Ingram et al. 1954). 

The reason for the reduced line width in the decolourizing charcoal 
is still not clear, but is thought to be due to the presence of adsorbed 
oxygen, and further measurements are in progress to clarify this point. 
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CXXXVII. Notices of New Books and Periodicals recewed 


The Concise History of Mathematics. By Dirk J. Srrurk. (G. Bell and Sons 

Ltd.) [Pp. xix+299.] Price 14s. 

Tuts book gives a.short account of the history of mathematics from the earliest 
times to the year 1900. The work is based on sound scholarship and is likely 
to become the standard book on the subject. The author has managed to 
give a balanced account of all the important periods in the development of 
mathematical thought. Owing to limitations of space the author is not in 
a position to discuss controversial questions at great length, or to point out 
the large number of unsolved problems which may offer scope for successful 
research. 

Though the author displays sound judgement in most places, the reviewer 
is not convinced by his arguments on the relatively undeveloped state of formal 
algebra in classical Greece. In particular the sentence “ If the classical authors 
had been interested in algebra they would have created the appropriate 
symbolism with which Diophantos actually made a beginning ” (p. 80) appears 
to be highly dogmatic. The author tries to point out relations between social 
changes and changes in the trend of mathematical research in various 
communities. As he has not the space at his disposal to develop these difficult 
problems in greater detail, these passages do not enlighten the reader to a — 
considerable extent. 

At the end of each chapter the author gives a good bibliography. The value 
of the book is further enhanced by a carefully compiled index. 

BAB. 


Optical Aberration Coefficients. By H. A. Bucupann. (Oxford: University 
Press.) [Pp. 336.] Price 50s. 


Royal Society Mathematical Tables. 3. Binomial Coefficients. (Cambridge : 
University Press.) [Pp. 162.] Price 35s. 


Royal Society Shorter Mathematical Tables. 3. Bessel Functions of Integer Orders 
and Large Arguments. By L. Fox. (Cambridge: University Press.) Price 6s. 6d. 


Elements of Statistical Mechanics. By D. ter Haar. (Rinehart & Company 
Inc.) [Pp. 468.] Price $8.50. 


Atomphysik III and IV. By K. Brcuerr. (Berlin: Walter de Gruyter.) 
Price D.M. 4.80 each. 


Statistical Analysis in Chemistry and the Chemical Industry. By Caru A. 
caer and Norman L. Frankury. (Wiley Publications in Statistics.) 
rice 58s. 


Sound. By R. C. Brown. (Longmans, Green and Co.) 


The Method of Trigonometrical Sums in the Theory of Numbers. By I. M. 
VinoeRaDoy. (Interscience Publishers, Ltd.) Price 33s. 


Gas Dynamics of Thin Bodies. By F. I. Franxu and E. A. Karpovicu. 
(Interscience Publishers, Ltd.) Price 38s. 


A Selection of Graphs for Use in Calculations of Compressible Airflow by the 
Compressible Flow Tables Panel. (Oxford : University Press.) Price 84s. 
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